// cyl1_slop.java   compile    javac -cp . cyl1_slop.java
//                  execute    java  -cp . cyl1_slop
//                             writes      cyl1_slop.dat
//                  execute    dat_to_stl cyl1_slop.dat cyl1_slop.stl
//
// curved version of cyl1_pipe

public class cyl1_slop  // make cyl1_slop.dat then .stl for 3D printer
{
  public cyl1_slop()
  {
    int num_points = 0;
    double data_points[][] = new double[5000][3];
    int num_polys = 0;
    int k = 0; // k index+1/4 = num_polys
    int data_polys[] = new int[5000]; // 3 pt1 pt2 pt3  +1 index
    double pi = 3.14159265358979323846;
    double theta, x, y, z;
    double phi, xp, yp, zp, del; // rotated points
    double r = 10.0;
    double rc = 30.0; // curv cyl
    int pts = 16;
    double dfct = (double)pts/2.0; // 0.0 to 2 pi
    double dtheta = pi/dfct;
    double dphi   = pi/16.0;
    int pti = 0;
    int pt0=0;      // outside 0 bottom
    int pt1=1*pts;  // outside 1
    int pt2=2*pts;  // outside 2
    int pt3=3*pts;  // outside 3
    int pt4=4*pts;  // outside 4 top
    int pt5=5*pts;  // inside  5 bottom
    int pt6=6*pts;  // inside  6
    int pt7=7*pts;  // inside  7
    int pt8=8*pts;  // inside  8
    int pt9=9*pts;  // inside  9 top
    
    
    System.out.println("cyl1_slop.java writing cyl1_slop.stl");
    System.out.println("r="+r+", pts="+pts+", dtheta="+dtheta);
    System.out.println("rc="+rc+", curv 5, dphi="+dphi);
    
    r = 10.0;
    for(int j=0; j<5; j++) // curv
    {
      phi = (double)j*dphi;
      for(int i=0; i<pts; i++) // generate vertices outside 0-4
      {
        theta = (double)i * dtheta;
        x = r * Math.sin(theta);
        y = r * Math.cos(theta);
        zp = rc*(1.0-Math.cos(phi));
	// zp = zp +0.5*y*Math.sin(phi);
        yp = y + rc*(1.0-Math.cos(phi));
        data_points[i+pti][0] = x;
        data_points[i+pti][1] = yp;
        data_points[i+pti][2] = zp;
        num_points++;
        // debug
	System.out.println("outside phi="+phi+", y="+y+", yp="+yp+",zp="+zp);
	System.out.println("pti="+pti+", x="+x);
	System.out.println("sin(phi)="+Math.sin(phi)+", cos(phi)="+Math.cos(phi));
        System.out.println(" ");
      } // end outside 0-4
      pti = pti + pts;
    } // end generate outside points

    System.out.println("added points, num_points="+num_points);

    System.out.println(" ");
    System.out.println("generate triangles outside");
    for(int i=0; i<pts; i++) // generate triangles outside 0-1
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt0+1;
      data_polys[k++] = i+pt1+2;
      data_polys[k++] = i+pt1+1;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt0+1;
      data_polys[k++] = i+pt0+2;
      data_polys[k++] = i+pt1+2;
    }
    // fix last two triangles
    data_polys[k-2] = pt0+1;
    data_polys[k-1] = pt1+1;
    data_polys[k-6] = pt1+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles outside 1-2
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt1+1;
      data_polys[k++] = i+pt2+2;
      data_polys[k++] = i+pt2+1;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt1+1;
      data_polys[k++] = i+pt1+2;
      data_polys[k++] = i+pt2+2;
    }
    // fix last two triangles
    data_polys[k-2] = pt1+1;
    data_polys[k-1] = pt2+1;
    data_polys[k-6] = pt2+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles outside 2-3
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt2+1;
      data_polys[k++] = i+pt3+2;
      data_polys[k++] = i+pt3+1;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt2+1;
      data_polys[k++] = i+pt2+2;
      data_polys[k++] = i+pt3+2;
    }
    // fix last two triangles
    data_polys[k-2] = pt2+1;
    data_polys[k-1] = pt3+1;
    data_polys[k-6] = pt3+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles outside 3-4
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt3+1;
      data_polys[k++] = i+pt4+2;
      data_polys[k++] = i+pt4+1;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt3+1;
      data_polys[k++] = i+pt3+2;
      data_polys[k++] = i+pt4+2;
    }
    // fix last two triangles
    data_polys[k-2] = pt3+1;
    data_polys[k-1] = pt4+1;
    data_polys[k-6] = pt4+1;
    num_polys = k/4;

    System.out.println("added triangles, k="+k+", num_polys="+num_polys);

    r = 8.0;
    for(int j=0; j<5; j++) // curv
    {
      phi = (double)j*dphi;
      for(int i=0; i<pts; i++) // generate vertices inside 5-9 
      {
        theta = (double)i * dtheta;
        x = r * Math.sin(theta);
        y = r * Math.cos(theta);
        zp = rc*(1.0-Math.cos(phi));
	// zp = zp + 0.5*y*Math.sin(phi);
        yp = y + rc*(1.0-Math.cos(phi));
        data_points[i+pti][0] = x;
        data_points[i+pti][1] = yp;
        data_points[i+pti][2] = zp;
        num_points++;
        // debug
	System.out.println("inside phi="+phi+", y="+y+", yp="+yp+", zp="+zp);
	System.out.println("pti="+pti+", x="+x);
	System.out.println("sin(phi)="+Math.sin(phi)+", cos(phi)="+Math.cos(phi));
        System.out.println(" ");
      }
      pti = pti + pts;
    } // end generate inside points

    System.out.println("added points, num_points="+num_points);

    System.out.println(" ");
    System.out.println("generate triangles inside");
    
    for(int i=0; i<pts; i++) // generate triangles inside 5-6
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt5+1;
      data_polys[k++] = i+pt6+1;
      data_polys[k++] = i+pt6+2;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt5+1;
      data_polys[k++] = i+pt6+2;
      data_polys[k++] = i+pt5+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt5+1;
    data_polys[k-2] = pt6+1;
    data_polys[k-5] = pt6+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles inside 6-7
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt6+1;
      data_polys[k++] = i+pt7+1;
      data_polys[k++] = i+pt7+2;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt6+1;
      data_polys[k++] = i+pt7+2;
      data_polys[k++] = i+pt6+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt6+1;
    data_polys[k-2] = pt7+1;
    data_polys[k-5] = pt7+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles inside 7-8
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt7+1;
      data_polys[k++] = i+pt8+1;
      data_polys[k++] = i+pt8+2;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt7+1;
      data_polys[k++] = i+pt8+2;
      data_polys[k++] = i+pt7+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt7+1;
    data_polys[k-2] = pt8+1;
    data_polys[k-5] = pt8+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles inside 8-9
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt8+1;
      data_polys[k++] = i+pt9+1;
      data_polys[k++] = i+pt9+2;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt8+1;
      data_polys[k++] = i+pt9+2;
      data_polys[k++] = i+pt8+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt8+1;
    data_polys[k-2] = pt9+1;
    data_polys[k-5] = pt9+1;
    num_polys = k/4;

    System.out.println("added triangles, k="+k+", num_polys="+num_polys);

    for(int i=0; i<pts; i++) // generate triangles bottom 0-5
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt0+1;
      data_polys[k++] = i+pt5+1;
      data_polys[k++] = i+pt5+2;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt0+1;
      data_polys[k++] = i+pt5+2;
      data_polys[k++] = i+pt0+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt0+1;
    data_polys[k-2] = pt5+1;
    data_polys[k-5] = pt5+1;
    num_polys = k/4;

    for(int i=0; i<pts; i++) // generate triangles top 4-9
    {
      data_polys[k++] = 3;
      data_polys[k++] = i+pt4+1;
      data_polys[k++] = i+pt9+2;
      data_polys[k++] = i+pt9+1;
      data_polys[k++] = 3;
      data_polys[k++] = i+pt4+1;
      data_polys[k++] = i+pt4+2;
      data_polys[k++] = i+pt9+2;
    }
    // fix last two triangles
    data_polys[k-1] = pt9+1;
    data_polys[k-2] = pt4+1;
    data_polys[k-6] = pt9+1;
    num_polys = k/4;

    System.out.println("added triangles, k="+k+", num_polys="+num_polys);
    System.out.println(" ");
    
    new datwrite("cyl1_slop.dat", data_points, num_points, data_polys, num_polys);
    
  } // end class cyl1_slop

  public static void main (String[] args)
  {
    new cyl1_slop();
  }
} // end  cyl1_slop.java