// stl_to_dat.java   read  .stl file  3D into  .dat file

import java.io.*;
import java.util.*;

// needs  datwrite.java

// solid rect  sample  .stl  format
//   facet normal 0.0 0.0 -1.0
//     outer loop
//       vertex 1.0 1.0 1.0
//       vertex 5.0 4.0 1.0
//       vertex 1.0 4.0 1.0
//     endloop
//   endfacet
//   facet normal 0.0 0.0 -1.0
//     outer loop
//       vertex 1.0 1.0 1.0
//       vertex 5.0 1.0 1.0
//       vertex 5.0 4.0 1.0
//    endloop
//   endfacet
//...
// endsolid

public class stl_to_dat
{
  boolean debug = false;
  int num_points = 0;
  int num_polys = 0; // all will be triangles  3,pt1,pt2,pt3
  double x[] = new double[10000];
  double y[] = new double[10000];
  double z[] = new double[10000];
  double xyz[][] = new double[10000][3];
  double xmin, xmax, ymin, ymax, zmin, zmax;
  int data_polys[] = new int[10000];
  int npt = 0;
  int pts[] = new int[3];
  int ipt; // unique     
  StringTokenizer T;
  String token;
  String solid = "solid";
  String facet = "facet";
  String  normal = "normal";
  String  outer = "outer";
  String  loop = "loop";
  String  vertex = "vertex";
  String endloop = "endloop";
  String endfacet = "endfacet";
  String endsolid = "endsolid";
      
    
  stl_to_dat(String filein, String fileout) // callable
  {
    System.out.println("stl_to_dat.java reading "+filein);
    try
    {
      FileReader fin = new FileReader(filein);
      BufferedReader in = new BufferedReader(fin);
      String input_line;
      int len;

      while(true) // to  endsolid
      {
        input_line = in.readLine(); // 
        len = input_line.length();
        if(debug) System.out.println("input:"+input_line+"  len="+len);
        T = new StringTokenizer(input_line," ");
        token = T.nextToken(); // ?
        if(debug) System.out.println("token:"+token);
	if(token.equals(endsolid)) break;
	if(token.equals(endloop))
	{
	  data_polys[num_polys] = 3;
	  data_polys[num_polys+1] = pts[0];
	  data_polys[num_polys+2] = pts[1];
	  data_polys[num_polys+3] = pts[2];
	  num_polys += 4;
	  npt = 0;
	  continue;
	}
	if(token.equals(endfacet)) continue;
	if(token.equals(vertex))
	{
          token = T.nextToken(); // x
          if(debug) System.out.println("token:"+token);
          x[num_points] = Double.parseDouble(token);
          token = T.nextToken(); // y
          if(debug) System.out.println("token:"+token);
          y[num_points] = Double.parseDouble(token);
          token = T.nextToken(); // z
          if(debug) System.out.println("token:"+token);
          z[num_points] = Double.parseDouble(token);
	  num_points++;
	  ipt = unique();
	  pts[npt] = ipt; // one big per  .dat
	  npt++;
	}
	
      } // end while      
      in.close();
    } // end try
    catch(IOException exception)
    {
      System.out.println("datread.java "+exception);
    }

    System.out.println("unique data points "+num_points);
    xmin = x[0];
    xmax = x[0];
    ymin = y[0];
    ymax = y[0];
    zmin = z[0];
    zmax = z[0];
    
    for(int i=0; i<num_points; i++)
    {
      if(debug) System.out.println("x="+x[i]+", y="+y[i]+", z="+z[i]);
      xmin = Math.min(xmin, x[i]);	
      xmax = Math.max(xmax, x[i]);	
      ymin = Math.min(ymin, y[i]);	
      ymax = Math.max(ymax, y[i]);	
      zmin = Math.min(zmin, z[i]);	
      zmax = Math.max(zmax, z[i]);	
    }
    System.out.println("xmin= "+xmin+" xmax= "+xmax);
    System.out.println("ymin= "+ymin+" ymax= "+ymax);
    System.out.println("zmin= "+zmin+" zmax= "+zmax);
    System.out.println("polys "+(num_polys/4));
    for(int i=0; i<num_polys; i=i+4)
    {
	if(debug) System.out.println("pt "+(i/4)+"  "+data_polys[i]+"  "+
				     data_polys[i+1]+"  "+data_polys[i+2]+
				     "  "+data_polys[i+3]);
    } // end i
    for(int i=0; i<num_points; i++)
    {
      xyz[i][0] = x[i];
      xyz[i][1] = y[i];
      xyz[i][2] = z[i];
    }
    
    new datwrite(fileout, xyz, num_points, data_polys, num_polys/4);
  } // end stl_to_dat

  int unique() // keep unique x,y,z
  {
    if(num_points<2) return 1;
    for(int i=0; i<num_points-1; i++)
    {
      if(Math.abs(x[num_points-1]-x[i])<0.001 &&
	 Math.abs(y[num_points-1]-y[i])<0.001 &&
	 Math.abs(z[num_points-1]-z[i])<0.001)
      {
	num_points--;
	return i+1; // one too big
      }
    }
    return num_points;
  } // end unique
    
  public static void main (String[] args)
  {
    new stl_to_dat(args[0],args[1]); // .stl  .dat
  }
} // end class stl_to_dat