org.simantics.g3d/src/org/simantics/g3d/shape/Sphere.java

   1 /*******************************************************************************\r
   2  * Copyright (c) 2012, 2013 Association for Decentralized Information Management in\r
   3  * Industry THTH ry.\r
   4  * All rights reserved. This program and the accompanying materials\r
   5  * are made available under the terms of the Eclipse Public License v1.0\r
   6  * which accompanies this distribution, and is available at\r
   7  * http://www.eclipse.org/legal/epl-v10.html\r
   8  *\r
   9  * Contributors:\r
  10  *     VTT Technical Research Centre of Finland - initial API and implementation\r
  11  *******************************************************************************/\r
  12 package org.simantics.g3d.shape;\r
  13 \r
  14 import java.util.ArrayList;\r
  15 import java.util.List;\r
  16 \r
  17 import javax.vecmath.AxisAngle4d;\r
  18 import javax.vecmath.Vector3d;\r
  19 \r
  20 import org.simantics.g3d.math.MathTools;\r
  21 \r
  22 public class Sphere {\r
  23 \r
  24         public static Mesh create(double radius, int s, int p) {\r
  25                 if (s < 3 || p < 3 || radius < MathTools.NEAR_ZERO)\r
  26                         throw new IllegalArgumentException();\r
  27                 List<Vector3d> vertices = new ArrayList<Vector3d>((s-2)*p + 2);\r
  28                 List<Vector3d> normals = new ArrayList<Vector3d>(vertices.size());\r
  29                 List<Integer> indices = new ArrayList<Integer>(((s-3)*p*2 + p * 2)*3);\r
  30                 Vector3d v = new Vector3d(0.0,-radius,0.0);\r
  31                 Vector3d vp = new Vector3d();\r
  32                 for (int ip = 0; ip < p; ip++) {\r
  33                         if (ip == 0) {\r
  34                                 vertices.add(new Vector3d(0.0,-radius,0.0));\r
  35                         } else if (ip == p - 1) {\r
  36                                 vertices.add(new Vector3d(0.0, radius,0.0));\r
  37                                 int off = 1 + (ip-2)*s;\r
  38                                 for (int is = 0; is < s; is++) {\r
  39                                         indices.add(vertices.size() - 1);\r
  40                                         indices.add(is+off);\r
  41                                         if (is < s -1)\r
  42                                                 indices.add(is+off+1);\r
  43                                         else\r
  44                                                 indices.add(off);\r
  45                                         \r
  46                                         \r
  47                                 }\r
  48                         } else {\r
  49                                 AxisAngle4d aa = new AxisAngle4d(1, 0, 0, ((double)ip/(double)(p-1))*Math.PI);\r
  50                                 MathTools.rotate(MathTools.getQuat(aa), v, vp);\r
  51                                 for (int is = 0; is < s; is++) {\r
  52                                         aa = new AxisAngle4d(0, 1, 0, ((double)is/(double)s)*Math.PI*2);\r
  53                                         Vector3d vs = new Vector3d();\r
  54                                         MathTools.rotate(MathTools.getQuat(aa), vp, vs);\r
  55                                         vertices.add(vs);\r
  56                                 }\r
  57                                 if (ip == 1) {\r
  58                                         for (int is = 0; is < s; is++) {\r
  59                                                 indices.add(0);\r
  60                                                 if (is < s -1)\r
  61                                                         indices.add(is+2);\r
  62                                                 else\r
  63                                                         indices.add(1);\r
  64                                                 indices.add(is+1);\r
  65                                         }\r
  66                                 } else {\r
  67                                         int off = 1 + (ip-1)*s;\r
  68                                         for (int is = 0; is < s-1; is++) {\r
  69                                                 indices.add(off + is - s);\r
  70                                                 indices.add(off + is+1);\r
  71                                                 indices.add(off + is);\r
  72                                                 \r
  73                                                 \r
  74                                                 indices.add(off + is - s);\r
  75                                                 indices.add(off + is + 1 - s);\r
  76                                                 indices.add(off + is + 1);\r
  77                                                 \r
  78                                         }\r
  79                                         indices.add(off - 1);\r
  80                                         indices.add(off);\r
  81                                         indices.add(off + s - 1);\r
  82                                         \r
  83                                         indices.add(off -1);\r
  84                                         indices.add(off - s);\r
  85                                         indices.add(off);\r
  86                                         \r
  87                                 }\r
  88                         }\r
  89                 }\r
  90                 for (int i = 0; i < vertices.size(); i++) {\r
  91                         Vector3d n = new Vector3d(vertices.get(i));\r
  92                         n.normalize();\r
  93                         normals.add(n);\r
  94                 }\r
  95                 \r
  96                 return new Mesh(vertices,normals,indices);\r
  97                 \r
  98         }\r
  99         \r
 100         public static void main(String arg[]) {\r
 101                 Mesh s1 = create(1.0, 3, 3);\r
 102                 Mesh s2 = create(1.0, 4, 4);\r
 103                 System.out.println("debug " + s1 + " " + s2);\r
 104         }\r
 105 }\r