// shape7.java  plot 3d data  Utah Graphics  to shape7.dat shape7.stl
//                 javac -cp . shape7.java
//             _   java  -cp . shape7 > shape7_java.out
//
//
//
//

import java.io.*;

public class shape7  // only 2d lines first cut, z=0
{
  boolean debug = true;
  int pts, pt, np, j, kk;
  int k = 0; // index into newpts k++
  datread DR; // = new datread(filename); // defines data structures
  int new_points = 0; // for each x,y,z point kept
  int new_polys = 0; // k/4 for datwrite
  double xyz[][] = new double[10000][3]; // x, y, z for output
  int newpts[] = new int[40000]; // 2, xyz1, xyz2 for output
  String fileout = "shape7.dat";
  int new_points2 = 0; // for each x,y,z 2point kept
  int new_polys2 = 0; // k/3 for datwrite
  double xyz2[][] = new double[10000][3]; // x, y, z for output
  int newpts2[] = new int[40000]; // 2, xyz1, xyz2 for output
  String fileout2 = "shape7a.dat";
  int newpts3[] = new int[40000]; // 2, xyz1, xyz2 for output
  int new_polys3 = 0; // k/4 for datwrite
  String fileout3 = "shape7b.dat";
  int new_points4 = 0; // for each x,y,z 2point kept
  int new_polys4 = 0; // k/3 for datwrite
  double xyz4[][] = new double[10000][3]; // x, y, z for output
  int newpts4[] = new int[40000]; // 2, xyz1, xyz2 for output
  String fileout4 = "shape7c.dat";
  double Pi =  Math.PI;
  double ang = 0.0;
  double dang = 0.0;
  double r1 = 10.0; // center at 0.0
  double r2 = 20.0;
  double r3 = 30.0;
  int n=32;
  double x[] = new double[8];  // for inside
  double y[] = new double[8];
                               // 0 to 7, 1 to 2, 3 to 4, 5 to 6 
  double xx[] = new double[8]; // for outside x,y to xx,yy
  double yy[] = new double[8];
  double ra, rb;
  int mov[] = {67,68,69,70,71,100,99,101,102,
               75,76,77,78,79,104,103,105,106,
               83,84,85,86,87,108,107,109,110,
	       91,92,93,94,95,112,111,97,98};
  int tri[][] = {{ 1,51, 2},{ 2,51, 3},{ 3,51,23},
                 { 3,23, 4},{ 4,23, 5},{ 5,23,24},
		 { 5,24, 6},{ 6,24, 7},{ 7,24,32},
		 { 7,32, 8},{ 8,32, 9},{ 9,32,33},
		 { 9,33,10},{10,33,11},{11,33,41},
		 {11,41,12},{12,41,13},{13,41,42},
		 {13,42,14},{14,42,15},{15,42,50},
		 {15,50,16},{16,50, 1},{ 1,50,51},
		 {51,52,17},{51,17,18},{51,18,19},
		 {51,19,23},{23,19,20},{23,20,21},
		 {23,21,22},{24,25,26},{24,26,27},
		 {24,27,28},{24,28,32},{32,28,29},
		 {32,29,30},{32,30,31},{33,34,35},
		 {33,35,36},{33,36,37},{33,37,41},
		 {41,37,38},{41,38,39},{41,39,40},
		 {42,43,44},{42,44,45},{42,45,46},
		 {42,46,50},{50,46,47},{50,47,48},
		 {50,48,49}};
		 
  public shape7()
  {
    double xc, yc, r, ang1, ang2, rxc, ryc;
    double pt1x, pt1y, pt2x, pt2y, dx, dy;
    int ca0, ca1, ca2, ca3, ca4, ca5, ca6, ca7; // inner arcs
    System.out.println("shape7.java running");
    new_polys = 0;
    new_points = 0;
    ra = Math.sqrt(r2*r2-r1*r1);
    rb = Math.sqrt(r3*r3-r1*r1);
    System.out.println("r1="+r1+"  r2="+r2+"  r3="+r3);
    System.out.println("ra="+ra+"  rb="+rb);
    
    System.out.println(" ");
    System.out.println("generate n points of pattern, inside");
    dang = 2.0*Pi/(double)n;
    ang = 0.0;
    for(int i=0; i<n; i++)
    {
      xyz[i][0] = r1*Math.cos(ang);
      xyz[i][1] = r1*Math.sin(ang);
      xyz[i][2] = r1/2.0;
      new_points = new_points+1;
      ang = ang + dang;
    } // end i inside
    System.out.println("new_points="+new_points+" generated, inside");
    System.out.println(" ");

    System.out.println("connect inside lines ");
    for(int i=0; i<n; i++)
    {
      newpts[k++] = 2;
      newpts[k++] = i+1;
      newpts[k++] = i+2;
    } // end i
    newpts[k-1] = 1; // close inside
    System.out.println("k="+k+" inside lines connected ");
    System.out.println(" ");

    System.out.println("generate middle points ");
    ang = 0.0;
    pt = new_points+1;
    for(int i=0; i<n; i++)
    {
      xyz[new_points][0] = r2*Math.cos(ang);
      xyz[new_points][1] = r2*Math.sin(ang);
      xyz[new_points][2] = r1/2.0;
      new_points = new_points+1;
      ang = ang + dang;
    } // end i
    System.out.println("new_points="+new_points+" generated, middle");
    System.out.println(" ");
    
    System.out.println("generate middle lines ");
    for(int i=0; i<n; i++)
    {
      newpts[k++] = 2;
      newpts[k++] = pt+i;
      newpts[k++] = pt+i+1;
    } // end i
    newpts[k-1] = pt; // close middle
    System.out.println("k="+k+" middle lines generated ");
    System.out.println(" ");


    System.out.println("generate outside points ");
    ang = 0.0;
    pt = new_points+1;
    for(int i=0; i<n; i++)
    {
      xyz[new_points][0] = r3*Math.cos(ang);
      xyz[new_points][1] = r3*Math.sin(ang);
      xyz[new_points][2] = r1/2.0;
      new_points = new_points+1;
      ang = ang + dang;
    } // end i
    System.out.println("new_points="+new_points+" generated, outside");
    System.out.println(" ");
    
    System.out.println("generate outside lines ");
    for(int i=0; i<n; i++)
    {
      newpts[k++] = 2;
      newpts[k++] = pt+i;
      newpts[k++] = pt+i+1;
    } // end i
    newpts[k-1] = pt; // close outside
    System.out.println("k="+k+" outside lines generated ");
    System.out.println(" ");

    System.out.println("generate lines ");

    x[0] = ra;
    y[0] = r1;
    xx[0] = rb;
    yy[0] = r1;

    x[1] = r1;
    y[1] = ra;
    xx[1] = r1;
    yy[1] = rb;

    x[2] = -r1;
    y[2] = ra;
    xx[2] = -r1;
    yy[2] = rb;

    x[3] = -ra;
    y[3] = r1;
    xx[3] = -rb;
    yy[3] = r1;

    x[4] = -ra;
    y[4] = -r1;
    xx[4] = -rb;
    yy[4] = -r1;

    x[5] = -r1;
    y[5] = -ra;
    xx[5] = -r1;
    yy[5] = -rb;

    x[6] = r1;
    y[6] = -ra;
    xx[6] = r1;
    yy[6] = -rb;

    x[7] = ra;
    y[7] = -r1;
    xx[7] = rb;
    yy[7] = -r1;

    pt = new_points+1;
    for(int i=0; i<8; i++)
    {
      xyz[new_points][0] = x[i];   // pt
      xyz[new_points][1] = y[i];
      xyz[new_points][2] = r1/2.0;
      new_points = new_points+1;
      xyz[new_points][0] = xx[i];  // pt+1
      xyz[new_points][1] = yy[i];
      xyz[new_points][2] = r1/2.0;
      new_points = new_points+1;
    } // end i
    System.out.println("new_points="+new_points+" generated, lines");
    System.out.println(" ");

    System.out.println("generate lines ");
    for(int i=0; i<8; i++)
    {
      newpts[k++] = 2;
      newpts[k++] = pt+2*i;
      newpts[k++] = pt+2*i+1;
    } // end i
    System.out.println("k="+k+" lines generated ");
    System.out.println(" ");
    newpts[k++] = 2;
    newpts[k++] = pt+2;
    newpts[k++] = pt+4;
    newpts[k++] = 2;
    newpts[k++] = pt+6;
    newpts[k++] = pt+8;
    newpts[k++] = 2;
    newpts[k++] = pt+10;
    newpts[k++] = pt+12;
    newpts[k++] = 2;
    newpts[k++] = pt+14;
    newpts[k++] = pt;
    System.out.println("k="+k+" iner lines generated ");
    System.out.println(" ");

    new_polys = k/3;
    System.out.println("new_polys="+new_polys);
    new datwrite(fileout, xyz, new_points, newpts, new_polys);
    System.out.println(fileout+" written");
    System.out.println(" ");

    System.out.println("generate inside lines2 ");
    k = 0;
    for(int i=0; i<16; i++)
    {
      xyz2[new_points2][0] = xyz[2*i][0];
      xyz2[new_points2][1] = xyz[2*i][1];
      xyz2[new_points2][2] = xyz[2*i][2];
      new_points2 = new_points2+1;
    }
    System.out.println("new_points2="+new_points2+" generated, lines2");
    System.out.println(" ");

    System.out.println("connect inside lines2 ");
    for(int i=0; i<16; i++)
    {
      newpts2[k++] = 2;
      newpts2[k++] = i+1;
      newpts2[k++] = i+2;
    } // end i
    newpts2[k-1] = 1; // close inside
    System.out.println("k="+k+" inside lines2 connected ");
    System.out.println(" ");


    System.out.println("generate outside points2 ");
    pt = new_points2;
    for(int i=0; i<36; i++)
    {
      xyz2[new_points2][0] = xyz[mov[i]-1][0];
      xyz2[new_points2][1] = xyz[mov[i]-1][1];
      xyz2[new_points2][2] = xyz[mov[i]-1][2];
      new_points2 = new_points2+1;
    }
    System.out.println("new_points2="+new_points2+" generated, outside lines2");
    System.out.println(" ");

    System.out.println("connect outsid2 lines ");
    for(int i=0; i<36; i++)
    {
      newpts2[k++] = 2;
      newpts2[k++] = pt+i+1;
      newpts2[k++] = pt+i+2;
    } // end i
    newpts2[k-1] = pt+1; // close outside
    System.out.println("k="+k+" inside lines connected ");
    System.out.println(" ");

    new_polys2 = k/3;
    System.out.println("new_polys2="+new_polys2);
    new datwrite(fileout2, xyz2, new_points2, newpts2, new_polys2);
    System.out.println(fileout2+" written");
    System.out.println(" ");

    System.out.println("generate triangles to cover shape ");
    System.out.println(" ");
    k = 0; // newpts3
    for(int i=0; i<52; i++) // triangle cover
    {
      newpts3[k++] = 3;
      newpts3[k++] = tri[i][0];
      newpts3[k++] = tri[i][1];
      newpts3[k++] = tri[i][2];
    }
    
    new_polys3 = k/4;
    System.out.println("new_polys2="+new_polys3);
    new datwrite(fileout3, xyz2, new_points2, newpts3, new_polys3);
    System.out.println(fileout3+" written");
    System.out.println(" ");

   
    System.out.println("copy for z=0");
    System.out.println(" ");
    k = 0;
    for(int i=0; i<new_points2; i++) // from previous
    {
      xyz4[new_points4][0] = xyz2[i][0];
      xyz4[new_points4][1] = xyz2[i][1];
      xyz4[new_points4][2] = xyz2[i][2];
      new_points4 = new_points4+1;
    }
    for(int i=0; i<new_points2; i++) // from previous z=0
    {
      xyz4[new_points4][0] = xyz2[i][0];
      xyz4[new_points4][1] = xyz2[i][1];
      xyz4[new_points4][2] = 0.0;
      new_points4 = new_points4+1;
    }
    System.out.println("new_points4="+new_points4+" generated, z=0");
    System.out.println(" ");

    System.out.println("generate triangles for z=0 ");
    System.out.println(" ");
    for(int i=0; i<new_polys3; i++) // triangles z=0 copy
    {
      newpts4[k++] = newpts3[4*i];
      newpts4[k++] = newpts3[4*i+1];
      newpts4[k++] = newpts3[4*i+2];
      newpts4[k++] = newpts3[4*i+3];
    }
 
    for(int i=0; i<new_polys3; i++) // triangles z=0 copy
    {
      newpts4[k++] = newpts3[4*i];
      newpts4[k++] = newpts3[4*i+3]+new_polys3; // reverse order
      newpts4[k++] = newpts3[4*i+2]+new_polys3;
      newpts4[k++] = newpts3[4*i+1]+new_polys3;
    }
    np = new_points2; // 52
    System.out.println("connect z=0 z=z np="+np);
    System.out.println(" ");
    for(int i=0; i<16; i++) // triangles z=0 z=z inside
    {
      newpts4[k++] = 3; 
      newpts4[k++] = i+1; 
      newpts4[k++] = i+np+1;
      newpts4[k++] = i+np+2;  // k-5 at close
      newpts4[k++] = 3; 
      newpts4[k++] = i+np+2;  // k-3 at close
      newpts4[k++] = i+2;     // k-2 at close
      newpts4[k++] = i+1;     // k-1 at close
    }
    newpts4[k-5] = np+1; // close
    newpts4[k-3] = np+1;
    newpts4[k-2] = 2;
    newpts4[k-1] = 1;
    for(int i=16; i<52; i++) // triangles z=0 z=z outside
    {
      newpts4[k++] = 3; 
      newpts4[k++] = i+1; 
      newpts4[k++] = i+np+1;
      newpts4[k++] = i+np+2;
      newpts4[k++] = 3; 
      newpts4[k++] = i+np+2;
      newpts4[k++] = i+2;
      newpts4[k++] = i+1;     // k-1 at close
    }
    newpts4[k-5] = 17+np; // close
    newpts4[k-3] = 17+np;
    newpts4[k-2] = 18;
    newpts4[k-1] = 17;

    new_polys4 = k/4;
    System.out.println("new_polys4="+new_polys4);
    new datwrite(fileout4, xyz4, new_points4, newpts4, new_polys4);
    System.out.println(fileout4+" written");
    System.out.println(" ");

    System.out.println("shape7.java finished");

  } // shape7 constructor

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