// shape8.java  plot 3d data  Utah Graphics  to shape8.dat
//                  javac -cp . shape8.java
//                  java  -cp . shape8 > shape8_java.out
//                  java  -cp . dat_to_stl shape8.dat shape8.stl
//                      
//  leafs    /\
//           \/  this shape, right side, defined by points
//           |                           then reverse x to -x
//        <>-+-<>
//           |      then rotate 90, 180, 270
//           /\
//           \/

import java.io.*;

public class shape8
{
  boolean debug = true;
  int pts, pt, n2, n3, n4, n5, n6;
  int new_points = 0; // for each x,y,z point kept
  int new_polys = 0; // changed k/4
  int k = 0;         // index for newpts
  double xyz[][] = new double[10000][3]; // x, y, z for output
  int newpts[] = new int[40000]; // 3, xyz1, xyz2, xyz3 for output
  String fileout = "shape8.dat";
  String fileouta = "shape8a.dat";
  String fileoutc = "shape8c.dat";
  double Pi = Math.PI;
  double r;
  double ang;
  int np; // first new_points
  int n = 40;
  double x, y, z, dz;
  double sxy[][] = new double[40][2];
  int ntri = 39;
  int tri[][] = {{ 2, 4, 3},{ 4, 6, 5},{ 2, 6, 4},
                { 8,10, 9},{10,12,11},{ 7,10, 8},
                { 7,12,10},{12,14,13},{16,18,17},
                {18,20,19},{20,22,21},{22,24,23},
                {20,25,22},{25,24,22},{18,25,20},
                {25,27,26},{27,29,28},{18,27,25},
                {16,27,18},{16,29,27},{16,30,29},
                {31,33,32},{33,36,34},{34,36,35},
                {31,37,33},{37,36,33},{31,37,36},
                {31,39,38},{15,30,16},{15,31,30},
                {15,39,31},{39,31,15},{ 7,14,12},
		{39,15,14},{39,14, 7},{31,38,37},
                {39, 7, 6},{39, 6, 2},{39, 2, 1}}; // 39
  int tri2[][] = {{45,47,46},{48,50,49},{45,48,47},
                  {45,50,48},{57,59,58},{56,59,57},
                  {56,60,59},{54,56,55},{52,54,53},
                  {56,61,60},{54,61,56},{54,62,61},
                  {54,63,62},{52,63,54},{63,65,64},
                  {52,65,63},{51,65,52},{51,66,65},
                  {50,66,51},{71,73,72},{69,71,70},
                  {67,69,68},{69,73,71},{69,74,73},
                  {75,77,76},{77,79,78},{75,79,77},
                  {42,80,79},{42,79,75},
                  {42,75,74},{42,74,69},{42,69,67},
                  {42,67,66},{42,66,50},{42,50,45},
                  {42,45,44},{42,44,43}};
  int ntri2=37;

  public shape8()
  {
    System.out.println("shape8.java running");
    new_polys = 0;
    k = 0;
    new_points = 0;

    System.out.println("sphere8.java running ");
    System.out.println(" ");
    System.out.println("generate top side_right points ");
    sxy[0][0]=0.0;     // side_r at z = 0.0
    sxy[0][1]=28.0;
    sxy[1][0]=0.5;
    sxy[1][1]=26.5;
    sxy[2][0]=1.0;
    sxy[2][1]=26.5;
    sxy[3][0]=1.5;
    sxy[3][1]=26.0;
    sxy[4][0]=2.0;
    sxy[4][1]=26.5;
    sxy[5][0]=1.5;
    sxy[5][1]=24.0;
    sxy[6][0]=2.0;
    sxy[6][1]=22.5;
    sxy[7][0]=3.0;
    sxy[7][1]=23.0;
    sxy[8][0]=4.0;
    sxy[8][1]=23.5;
    sxy[9][0]=3.5;
    sxy[9][1]=22.5;
    sxy[10][0]=4.0;
    sxy[10][1]=22.0;
    sxy[11][0]=3.0;
    sxy[11][1]=20.5;
    sxy[12][0]=3.5;
    sxy[12][1]=20.0;
    sxy[13][0]=3.0;
    sxy[13][1]=18.0;
    sxy[14][0]=3.0;
    sxy[14][1]=17.0;
    sxy[15][0]=4.0;
    sxy[15][1]=18.0;
    sxy[16][0]=5.0;
    sxy[16][1]=20.0;
    sxy[17][0]=5.5;
    sxy[17][1]=19.5;
    sxy[18][0]=6.5;
    sxy[18][1]=21.0;
    sxy[19][0]=7.0;
    sxy[19][1]=21.0;
    sxy[20][0]=7.5;
    sxy[20][1]=21.5;
    sxy[21][0]=8.0;
    sxy[21][1]=21.5;
    sxy[22][0]=9.0;
    sxy[22][1]=22.5;
    sxy[23][0]=8.5;
    sxy[23][1]=21.0;
    sxy[24][0]=8.0;
    sxy[24][1]=20.5;
    sxy[25][0]=8.5;
    sxy[25][1]=20.0;
    sxy[26][0]=8.0;
    sxy[26][1]=19.0;
    sxy[27][0]=9.0;
    sxy[27][1]=18.0;
    sxy[28][0]=8.0;
    sxy[28][1]=16.5;
    sxy[29][0]=6.5;
    sxy[29][1]=15.0;
    sxy[30][0]=6.0;
    sxy[30][1]=14.0;
    sxy[31][0]=8.5;
    sxy[31][1]=14.0;
    sxy[32][0]=8.5;
    sxy[32][1]=13.0;
    sxy[33][0]=9.0;
    sxy[33][1]=12.5;
    sxy[34][0]=11.0;
    sxy[34][1]=11.5;
    sxy[35][0]=8.0;
    sxy[35][1]=11.5;
    sxy[36][0]=6.5;
    sxy[36][1]=10.5;
    sxy[37][0]=6.0;
    sxy[37][1]=10.0;
    sxy[38][0]=0.5;
    sxy[38][1]=10.5;
    sxy[39][0]=0.5;
    sxy[39][1]=0.5;
    n = 40;

    System.out.println("top_r ");
    z = 0.0;
    for(int i=0; i<n; i++) // top_r  1..n
    {
      x = sxy[i][0];
      y = sxy[i][1];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("top_r xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end 
    System.out.println("top_r points="+new_points);
    System.out.println(" ");

    System.out.println("shift 90 side_r top");
    n2 = n;
    for(int i=1; i<=n; i++) // shift 90 side_r top  1..n
    {
      x = sxy[n2-i][1];
      y = sxy[n2-i][0];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("side_r top xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("side_r top points="+new_points);
    System.out.println(" ");

    System.out.println("top_r triangle test"); // 1 to 39
    for(int i=0; i<ntri; i++)
    {
      newpts[k++] = 3;
      newpts[k++] = tri[i][0];
      newpts[k++] = tri[i][1];
      newpts[k++] = tri[i][2];
      System.out.println("tri="+newpts[k-3]+" "+
               	         newpts[k-2]+" "+newpts[k-1]);
    } // end i

    System.out.println("side_r top triangle test"); // 42 to 80 
    for(int i=0; i<ntri2; i++)
    {
      newpts[k++] = 3;
      newpts[k++] = tri2[i][0];
      newpts[k++] = tri2[i][1];
      newpts[k++] = tri2[i][2];
      System.out.println("tri="+newpts[k-3]+" "+
               	         newpts[k-2]+" "+newpts[k-1]);
    } // end i
    
    new_polys = k/4;
    System.out.println("new_points="+new_points+"  new_polys="+new_polys);

    new datwrite(fileoutc, xyz, new_points, newpts, new_polys);
    System.out.println("shape8c.dat written");
    System.out.println(" ");
    k = 0; // end test triangle


    System.out.println("shift 90 side_r bot");
    for(int i=0; i<n; i++)
    {
      x = sxy[i][1];
      y = -sxy[i][0];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("side_r bot xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("side_r bot points="+new_points);
    System.out.println(" ");


    System.out.println("shift 180 bot_r");
    n2 = n;
    for(int i=1; i<=n; i++) // bot_r  1..n
    {
      x = sxy[n2-i][0];
      y = -sxy[n2-i][1];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("bot_r bot xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("bot_r points="+new_points);
    System.out.println(" ");
    
    System.out.println("shift 180 bot_l");
    for(int i=0; i<n; i++)
    {
      x = -sxy[i][0];
      y = -sxy[i][1];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("bot_l bot xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("bot_l points="+new_points);
    System.out.println(" ");

    System.out.println("shift 270 side_l bot");
    n2 = n;
    for(int i=1; i<=n; i++) // shift 270 side_l bot
    {
      x = -sxy[n2-i][1];
      y = -sxy[n2-i][0];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("side_l bot xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("side_l top points="+new_points);
    System.out.println(" ");

    System.out.println("shift 270 side_l top");
    for(int i=0; i<n; i++)
    {
      x = -sxy[i][1];
      y = sxy[i][0];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("side_l top xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("side_l bot points="+new_points);
    System.out.println(" ");


    n2 = n;
    System.out.println("top_l ");
    for(int i=1; i<=n; i++) // top_l  1..n
    {
      x = -sxy[n2-i][0];
      y = sxy[n2-i][1];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("top_l  xyz="+new_points+" "+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
    } // end i
    System.out.println("top_l points="+new_points);
    System.out.println(" ");


    for(int i=0; i<8*n-1; i++)
    {
      newpts[k++] = 3;
      newpts[k++] = i+1;
      newpts[k++] = i+2;
      newpts[k++] = i+2; // just lines
    } // end i
    

    new_polys = k/4;
    System.out.println("new_points="+new_points+"  new_polys="+new_polys);

    new datwrite(fileouta, xyz, new_points, newpts, new_polys);
    System.out.println("shape8a.dat written");
    System.out.println(" ");

    np = new_points;
    // top z=5
    z = 5.0;
    k = 0;
    System.out.println("copy top ");
    // copy points, newpts, add triangles to sides
    for(int i=0; i<np; i++)
    {
      xyz[new_points][0] = xyz[i][0];
      xyz[new_points][1] = xyz[i][1]; 
      xyz[new_points][2] = z;
      new_points = new_points +1;
    
      newpts[k++] = 3;
      newpts[k++] = i+1;
      newpts[k++] = i+2;
      newpts[k++] = i+1+np;
      System.out.println("i="+i+" newpts= "+newpts[k-4]+" "+newpts[k-3]+
                         " "+newpts[k-2]+" "+newpts[k-1]);
      newpts[k++] = 3;
      newpts[k++] = i+1+np;
      newpts[k++] = i+2;
      newpts[k++] = i+2+np;
      System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
                         " "+newpts[k-2]+" "+newpts[k-1]);
    } // end i
    newpts[k-5] = 1;
    newpts[k-3] = 1;
    newpts[k-1] = 2;
    System.out.println("k="+k+" newpts= "+newpts[k-4]+" "+newpts[k-3]+
                         " "+newpts[k-2]+" "+newpts[k-1]);


    System.out.println("triangle top_bot"); // first of 2 each doing 4
    for(int j=0; j<4; j++)
    {
      for(int i=0; i<ntri; i++)
      {
        newpts[k++] = 3;
        newpts[k++] = tri[i][0]+80*j;
        newpts[k++] = tri[i][1]+80*j;
        newpts[k++] = tri[i][2]+80*j;
        newpts[k++] = 3;
        newpts[k++] = tri[i][0]+np+80*j; // other z
        newpts[k++] = tri[i][1]+np+80*j;
        newpts[k++] = tri[i][2]+np+80*j;
      } // end i
    } // end j

    System.out.println("triangle top_bot"); // second of 2 each doing 4
    // n2 = ntri;
    for(int j=0; j<4; j++)
    {
      for(int i=0; i<ntri2; i++)
      {
        newpts[k++] = 3;
        newpts[k++] = tri2[i][0]+80*j;
        newpts[k++] = tri2[i][1]+80*j;
        newpts[k++] = tri2[i][2]+80*j;
        newpts[k++] = 3;
        newpts[k++] = tri2[i][0]+np+80*j; // other z
        newpts[k++] = tri2[i][1]+np+80*j;
        newpts[k++] = tri2[i][2]+np+80*j;
      } // end i
    } // end j



    System.out.println("center vertical"); // triangles center vertical
    newpts[k++] = 3;
    newpts[k++] = 39;
    newpts[k++] = 282;
    newpts[k++] = 199;
    newpts[k++] = 3;
    newpts[k++] = 39+np; // other z
    newpts[k++] = 282+np;
    newpts[k++] = 199+np;
    newpts[k++] = 3;
    newpts[k++] = 199;
    newpts[k++] = 122;
    newpts[k++] = 39;
    newpts[k++] = 3;
    newpts[k++] = 199+np; // other z
    newpts[k++] = 122+np;
    newpts[k++] = 39+np;
    System.out.println("center horizontal"); // triangles center horizontal
    newpts[k++] = 3;
    newpts[k++] = 42;
    newpts[k++] = 279;
    newpts[k++] = 202;
    newpts[k++] = 3;
    newpts[k++] = 42+np; // other z
    newpts[k++] = 279+np;
    newpts[k++] = 202+np;
    newpts[k++] = 3;
    newpts[k++] = 202;
    newpts[k++] = 119;
    newpts[k++] = 42;
    newpts[k++] = 3;
    newpts[k++] = 202+np; // other z
    newpts[k++] = 119+np;
    newpts[k++] = 42+np;

    System.out.println("center top"); // triangles
    newpts[k++] = 3;
    newpts[k++] = 1;
    newpts[k++] = 282;
    newpts[k++] = 39;
    newpts[k++] = 3; // other z
    newpts[k++] = 1+np;
    newpts[k++] = 282+np;
    newpts[k++] = 39+np;
    System.out.println("center right"); // triangles 
    newpts[k++] = 3;
    newpts[k++] = 80;
    newpts[k++] = 42;
    newpts[k++] = 119;
    newpts[k++] = 3; // other z
    newpts[k++] = 80+np;
    newpts[k++] = 42+np;
    newpts[k++] = 119+np;
    System.out.println("center bottom"); // triangles
    newpts[k++] = 3;
    newpts[k++] = 160;
    newpts[k++] = 122;
    newpts[k++] = 199;
    newpts[k++] = 3; // other z
    newpts[k++] = 160+np;
    newpts[k++] = 122+np;
    newpts[k++] = 199+np;
    System.out.println("center left"); // triangles 
    newpts[k++] = 3;
    newpts[k++] = 240;
    newpts[k++] = 202;
    newpts[k++] = 279;
    newpts[k++] = 3;
    newpts[k++] = 240+np;
    newpts[k++] = 202+np;
    newpts[k++] = 279+np;

    
    new_polys = k/4;
    System.out.println("new_points="+new_points+"  new_polys="+new_polys);

    new datwrite(fileout, xyz, new_points, newpts, new_polys);
    System.out.println("shape8.dat written");
    System.out.println(" ");
    System.out.println("shape8.java finished");

  } // shape8 constructor

  public static void main(String[] args)
  {
    new shape8(); // writes shape8.dat
  }
} // end class  shape8.java