// test_datread.c
#include <stdio.h>
#include <stdlib.h>
#include "datread.h"

static int num_points;            // for datread and datwrite 
static int num_polys;
static dpts * data_points;        // vertex points
static int *  data_polys;         // polygons
static float size;
static int status = -1;

int main(int argc, char *argv[])
{
  int i, j, k, pts, pt;

  printf("test_datread.c running \n"); 
  // read a .dat file 
  if(argc>1)
  {
    status = datread(argv[1], &data_points, &num_points,
                              &data_polys, &num_polys, &size);
    printf("read %s, %d vert, %d poly, %f size\n", 
	   argv[1], num_points, num_polys, size); 
  }
  else
  {
    printf("a filename of a .dat file is required.\n");
    return 0;
  }
  if(status!=0)
  {
    printf("%s file not read as a .dat file\n", argv[1]);
    return 0;
  }

  printf("\n vertex points \n");
  for(i=0; i<num_points; i++)
  {
    printf("%f %f %f \n",
           data_points[i].x, data_points[i].y, data_points[i].z); 
  }

  printf("\n surface polygons \n");
  k = 0; // polygons
  for(i=0; i<num_polys; i++)
  {
    pts = data_polys[k++];
    printf("%3d ", pts); // no cr
    for(j=0; j<pts; j++)
    {
      pt = data_polys[k++];
      printf("%d ", pt);
    }
    printf(" \n");
  } // end print polys

  status = datclosed(&data_points, &num_points,
                     &data_polys, &num_polys);
  if(status==1)
    printf("datclosed found polyhedra closed \n");
  else
    printf("datclosed found polyhedra not closed \n");

  printf("test_datread.c finished \n"); 
  return 0;
} // end main test_datread.c