// gear12.java       compile    javac -cp . gear12.java  2D to3D  new z=10.0
//                   execute    java  -cp . gear12   writes gear12d3.java
//                   reads gear12.dat
//                   comes with center hole 12 points, make triangles
// there are 16 points in each tooth, 192,   1,   2, ...  16
//                                          17,  18, ...  32 
//                                          33,  34, ...  48
//                                 ...     177, 178, ... 192  bot z=0
//       217 ...408                                           top z=10
// hole 193, 194, ... 204  --- to 1,  16,  17, 204  3 triangles bot
//      409, 410, ... 420       217, 232, 233, 408              top 
// tooth center 205, 206, ... 216                               bot
//              421, 422, ... 432                               top
public class gear12 
{
  int numteeth = 12;
  int ptteeth = 12*16; // plus numteeth tooth centers and numteeth hole
  double xyz[][] = new double[2*(216+8)][3]; // 216+8 layer plus bottom 
  double xyzt[][] = new double[2*(216+8)][3]; // gear12tt.dat top 
  int polyst[] = new int[4*4*(216+8)];
  int num_points = 204; // coming in from gear12.dat 1..192 outline
                        // 193..204 hole
  int num_polys = 204;  // coming in from gear12.dat 
  int num_tpoints = 0;  // for gear12t.dat + 12 tooth centers
  int num_tpoly = 0;    // computed
  int new_points = 0;   // computed
  int new_polys = 0;    // computed
  int newt_polys = 0;
  int newt_points = 0;
  int oldpolys[] = new int[3*204];  // just 2 1 2
  int polys[] = new int[4*4*244];   // 3 1 2 205, double times 2 for outside
  String filename = "gear12.dat";   // datread
  String filetri  = "gear12t.dat";  // datwrite with triangles
  String filetrit  = "gear12tt.dat";  // datwrite top z=10
  String fileout  = "gear12d3.dat"; // datwrite
  datread DR;
  double g12[][] = new double[12][3]; // for new tooth center points
  double r = 0.0;
  double rc = 22.73; // from gear12.dat, tooth 2
  int pc = 205;      // first of numteeth tooth centers
  double x = 0.0;
  double y = 0.0;
  double z = 10.0;
  double Pi = Math.PI;
  double ang = 0.0;
  double dang = -Pi/6.0; // numteeth/2 of tooth centers
  int k = 0;
  int kn = 0;
  int pt = 0;
  int pt2 = 0;
  int pt3 = 0;
  int pt4 = 0;
  int pt5 = 0;
  int pts = 0;
  int ptc = 204; // first added point
  int ptz = 204; // tooth center 
  double xmin, xmax, ymin, ymax, zmin, zmax;
    
  public gear12()
  {
    System.out.println("gear12,java running");

    DR = new datread(filename); // defines data structures
    num_points = DR.num_points;
    num_polys = DR.num_polys;
    System.out.println("num_points="+DR.num_points+
                       ", num_polys="+DR.num_polys);
    System.out.println("print vertices");
    xmin = DR.data_points[0].x;
    xmax = DR.data_points[0].x;
    ymin = DR.data_points[0].y;
    ymax = DR.data_points[0].y;
    zmin = DR.data_points[0].z;
    zmax = DR.data_points[0].z;
    for(int i=0; i<DR.num_points; i++) // vertices
    {
      System.out.println(+(i+1)+"  "+DR.data_points[i].x+
                         "  "+DR.data_points[i].y+"  "+DR.data_points[i].z);
      xyz[i][0] = DR.data_points[i].x;
      xyz[i][1] = DR.data_points[i].y;
      xyz[i][2] = DR.data_points[i].z;
      xmin = Math.min(xmin, DR.data_points[i].x);
      xmax = Math.max(xmax, DR.data_points[i].x);
      ymin = Math.min(ymin, DR.data_points[i].y);
      ymax = Math.max(ymax, DR.data_points[i].y);
      zmin = Math.min(zmin, DR.data_points[i].z);
      zmax = Math.max(zmax, DR.data_points[i].z);
    } // end print vertices
    System.out.println("xmin="+xmin+"   xmax="+xmax);
    System.out.println("ymin="+ymin+"   ymax="+ymax);
    System.out.println("zmin="+zmin+"   zmax="+zmax);
    
    System.out.println("print polygons");
    k = 0; // polygons
    for(int i=0; i<DR.num_polys; i++)
    {
      pts = DR.data_polys[k++];
      oldpolys[k-1] = pts;
      System.out.print(+pts); // no cr
      for(int j=0; j<pts; j++)
      {
        pt = DR.data_polys[k++];
        oldpolys[k-1] = pt;
        System.out.print(" "+ pt);
      }
      System.out.println(" ");
    } // end print polygons


    System.out.println("make triangles for gear12t1.dat");
    k = 0; // 3 point triangles to hole 193 to 204
    polys[k++] = 3;    // 193
    polys[k++] = 193;
    polys[k++] = 1;
    polys[k++] = 16;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 193;
    polys[k++] = 16;
    polys[k++] = 17;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 193;
    polys[k++] = 17;
    polys[k++] = 194;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      

    pt = 194;
    pt2 = 17;
    pt3 = 32;
    for(int i=0; i<10; i++)
    {  
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt2;
      polys[k++] = pt3;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt3;
      polys[k++] = pt3+1;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt3+1;
      polys[k++] = pt+1;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      pt = pt+1;
      pt2 = pt2+16;
      pt3 = pt3+16;
    } // end i
    polys[k++] = 3;
    polys[k++] = 204;
    polys[k++] = 177;
    polys[k++] = 192;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 204;
    polys[k++] = 192;
    polys[k++] = 1;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 204;
    polys[k++] = 1;
    polys[k++] = 193;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    new_polys = k/4;
    System.out.println("finished triangles to hole, new_polys="+new_polys);
    System.out.println("Make tooth centers");
    pt = num_points;
    pc = num_points;
    ang = 0.0;
    for(int i=0; i<12; i++)
    {
      x = rc*Math.cos(ang);
      y = rc*Math.sin(ang);
      z = 0.0;
      System.out.println("tooth center "+(i+1)+" pt="+pt+" x="+x+"  y="+y);
      xyz[pt][0] = x;
      xyz[pt][1] = y;
      xyz[pt][2] = z;
      ang = ang + dang; 
      pt += 1;
    }
    new_points = num_points + numteeth;
    System.out.println("12 tooth center points added, new_points="+new_points);
    pt = 205; // to 216
    pt2 = 1;
    pt3 = 1;
    pt4 = 16;
    for(int i=0; i<numteeth; i++) // tooth center
    {
      for(int j=0; j<15; j++) // points in tooth 
      {
        polys[k++] = 3;
        polys[k++] = pt;
        polys[k++] = pt2;
        polys[k++] = pt2+1;
        System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);
	pt2 = pt2+1;
      } // end j around tooth
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt4;
      polys[k++] = pt3;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      pt = pt+1;
      pt2 = pt2+1;
      pt3 = pt3+16;
      pt4 = pt4+16;
    } // end i tooth center

    new_polys = k/4;

    System.out.println("datwrite "+filetri+" new_points="+new_points+
                       " new_polys="+new_polys);
    new datwrite(filetri, xyz, new_points, polys, new_polys);
    System.out.println("datwrite wrote "+filetri);

    for(int i=0; i<new_points; i++)  // copy points  z=10.0
    {
      xyz[i+new_points][0] = xyz[i][0];
      xyz[i+new_points][1] = xyz[i][1];
      xyz[i+new_points][2] = 10.0;
      xyzt[i][0] = 0.0; // fill to get correct point numbers for tt
      xyzt[i][1] = 0.0;
      xyzt[i][2] = 0.0;
      xyzt[i+new_points][0] = xyz[i][0];
      xyzt[i+new_points][1] = xyz[i][1];
      xyzt[i+new_points][2] = 10.0;
    } // end i 

    System.out.println("new_polys="+new_polys+" k="+k);
    for(int i=0; i<k; i++) // copy polys z=10.0
    {
      if(i%4==0)
      {
        polys[i+k] = polys[i]; // 3
        polyst[i] = polys[i]; // 3
      }
      else
      { 
	polys[i+k] = polys[i]+new_points;
	polyst[i] = polys[i]+new_points;
      }
    } // end i

    newt_points = 2*(num_points+numteeth); // includes dill
    newt_polys = num_points+numteeth+12;
    System.out.println("datwrite "+filetrit+" newt_points="+newt_points+
		       " newt_polys="+newt_polys);
    new datwrite(filetrit, xyzt, newt_points, polyst, newt_polys);
    System.out.println("datwrite wrote "+filetrit);

    new_points = 2*new_points;
    System.out.println("new_points="+new_points);
    k = 2*k;
    new_polys = k/4;
    System.out.println("new_polys="+new_polys+"  k="+k);
    System.out.println("outside tooth 1-128 145-272");
    // outside tooth
    pt = 1;
    pt2 = 145;
    for(int i=0; i<127; i++) // one less, not tooth hole
    {
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt+1;
      polys[k++] = pt2;
      if(i>3 && i<125) System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      polys[k++] = 3;
      polys[k++] = pt+1;
      polys[k++] = pt2+1;
      polys[k++] = pt2;
      if(i>3 && i<125) System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      pt = pt+1;
      pt2 = pt2+1;
    } // end i  then 128 1 145  1 145 272
    polys[k++] = 3;
    polys[k++] = 128;
    polys[k++] = 1;
    polys[k++] = 145;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 1;
    polys[k++] = 145;
    polys[k++] = 272;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    System.out.println("outside hole 129-136 273-280");
    // outside hole 129-136  273-280
    pt = 129;
    pt2 = 273;
    for(int i=0; i<7; i++) // one less
    {
      polys[k++] = 3;
      polys[k++] = pt;
      polys[k++] = pt+1;
      polys[k++] = pt2;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      polys[k++] = 3;
      polys[k++] = pt+1;
      polys[k++] = pt2+1;
      polys[k++] = pt2;
      System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
      pt = pt+1;
      pt2 = pt2+1;
    } // end i  then 136 129 280  129 272 280
    polys[k++] = 3;
    polys[k++] = 136;
    polys[k++] = 129;
    polys[k++] = 280;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    polys[k++] = 3;
    polys[k++] = 129;
    polys[k++] = 272;
    polys[k++] = 280;
    System.out.println("polys["+(k-4)+"]="+polys[k-4]+" "+polys[k-3]+" "+polys[k-2]+" "+polys[k-1]);      
    
    new_polys = k/4;
    System.out.println("new_polys="+new_polys+"  k="+k);

    System.out.println("datwrite "+fileout+" new_points="+new_points+
                       " new_polys="+new_polys);
    new datwrite(fileout, xyz, new_points, polys, new_polys);
    System.out.println("datwrite wrote "+fileout);
    System.out.println("gear12.java ends");
  }

  public static void main (String[] args)
  {
    new gear12(); // filenames built in
  }
} // end gear12.java