// shape3.java  plot 3d data  Utah Graphics  to shape3.dat
//                  javac -cp . shape3.java
//           ___    java  -cp . shape3 > shape3_java.out
// |\___/|  |   |
// | ___ |  |   |
// |/   \|  |___|

import java.io.*;

public class shape3
{
  boolean debug = true;
  int pts, pt;
  datread DR; // = new datread(filename); // defines data structures
  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 = "shape3.dat";
  double edge, da;
  int n = 32;
  int nz = 19;
  int np = 0;
  double x, y, z, dz;
  double sxy[] = new double[19];     

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

    sxy[0]=2.0;     // sides at z = 0.0
    sxy[1]=1.7369;
    sxy[2]=1.4517;
    sxy[3]=1.2022;
    sxy[4]=0.9894;
    sxy[5]=0.8141;
    sxy[6]=0.6771;
    sxy[7]=0.5788;
    sxy[8]=0.5197;
    sxy[9]=0.5;     // sides at z = 2.5
    sxy[10]=0.5197;
    sxy[11]=0.5788;
    sxy[12]=0.6771;
    sxy[13]=0.8141;
    sxy[14]=0.9894;
    sxy[15]=1.2022;
    sxy[16]=1.4517;
    sxy[17]=1.7369;
    sxy[18]=2.0;    // sides at z = 5.0

    System.out.println(" ");
    System.out.println("generate new points ");

    dz = 5.0/(double)(nz-1);
    da = 4.0/(double)(n-1);
    System.out.println("edge1 ");
    z = 0.0;
    for(int i=0; i<nz; i++) // edge1  1..19 nz
    {
      x = -sxy[i];
      y = -sxy[i];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("edge1 xyz="+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
      if(i<nz-1)
      {
        newpts[k++] = 3;
        newpts[k++] = i+1;
        newpts[k++] = i+2;
        newpts[k++] = i+1+nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
                           " "+newpts[k-2]+" "+newpts[k-1]);
        newpts[k++] = 3;
        newpts[k++] = i+1+nz;
        newpts[k++] = i+2;
        newpts[k++] = i+2+nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
       	                     " "+newpts[k-2]+" "+newpts[k-1]);
      }
      z = z + dz;
    } // end i

    System.out.println("edge2 ");
    z = 0.0;
    for(int i=0; i<nz; i++) // edge2  20..38  nz+1 .. 2*nz
    {
      x = -sxy[i];
      y = sxy[i];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("edge2 xyz="+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
      if(i<nz-1)
      {
        newpts[k++] = 3;
        newpts[k++] = i+1+nz;
        newpts[k++] = i+2+nz;
        newpts[k++] = i+1+2*nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
                           " "+newpts[k-2]+" "+newpts[k-1]);
        newpts[k++] = 3;
        newpts[k++] = i+1+2*nz;
        newpts[k++] = i+2+nz;
        newpts[k++] = i+2+2*nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
       	                     " "+newpts[k-2]+" "+newpts[k-1]);
      }
      z = z + dz;
    } // end i

    System.out.println("edge3 ");
    z = 0.0;
    for(int i=0; i<nz; i++) // edge3  39..57  2*nz+1 .. 3*nz
    {
      x = sxy[i];
      y = sxy[i];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("edge3 xyz="+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
      if(i<nz-1)
      {
        newpts[k++] = 3;
        newpts[k++] = i+1+2*nz;
        newpts[k++] = i+2+2*nz;
        newpts[k++] = i+1+3*nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
                           " "+newpts[k-2]+" "+newpts[k-1]);
        newpts[k++] = 3;
        newpts[k++] = i+1+3*nz;
        newpts[k++] = i+2+2*nz;
        newpts[k++] = i+2+3*nz;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
       	                     " "+newpts[k-2]+" "+newpts[k-1]);
      }
      z = z + dz;
    } // end i

    System.out.println("edge4 "); 
    z = 0.0;
    for(int i=0; i<nz; i++) // edge4   58..76  3*nz+1 .. 4*nz
    {
      x = sxy[i];
      y = -sxy[i];
      xyz[new_points][0] = x;
      xyz[new_points][1] = y;
      xyz[new_points][2] = z;
      new_points = new_points +1;
      System.out.println("edge3 xyz="+xyz[new_points-1][0]+" "+
               	         xyz[new_points-1][1]+" "+xyz[new_points-1][2]);
      if(i<nz-1)
      {
        newpts[k++] = 3;
        newpts[k++] = i+1+3*nz;
        newpts[k++] = i+2;
        newpts[k++] = i+1;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
                           " "+newpts[k-2]+" "+newpts[k-1]);
        newpts[k++] = 3;
        newpts[k++] = i+1;
        newpts[k++] = i+2+3*nz;
        newpts[k++] = i+2;
        System.out.println("i="+i+" newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
       	                     " "+newpts[k-2]+" "+newpts[k-1]);
      }
      z = z + dz;
    } // end i

    // bottom and top square
    System.out.println("top and bottom ");
    newpts[k++] = 3;
    newpts[k++] = 1;
    newpts[k++] = 20;
    newpts[k++] = 39;
    System.out.println("bot1 newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
		       " "+newpts[k-2]+" "+newpts[k-1]); 
    newpts[k++] = 3;
    newpts[k++] = 39;
    newpts[k++] = 58;
    newpts[k++] = 1;
    System.out.println("bot2 newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
		       " "+newpts[k-2]+" "+newpts[k-1]); 

    newpts[k++] = 3;
    newpts[k++] = 19;
    newpts[k++] = 38;
    newpts[k++] = 57;
    System.out.println("top1 newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
		       " "+newpts[k-2]+" "+newpts[k-1]); 
    newpts[k++] = 3;
    newpts[k++] = 57;
    newpts[k++] = 76;
    newpts[k++] = 19;
    System.out.println("top2 newpts[k-4]="+newpts[k-4]+" "+newpts[k-3]+
		       " "+newpts[k-2]+" "+newpts[k-1]); 

    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("shape3.java finished");

  } // shape3 constructor

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