// make_spike_inp.c    read  spike.dat  write spike.inp UCD pde format

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include "datread.h"
#undef  abs
#define abs(x) ((x>=0.0)?x:(-x))
#undef  max
#define max(x,y) ((x)>(y)?(x):(y))

static int num_points;            // for datread and datwrite 
static int num_polys;
static dpts * data_points;
static int *  data_verts;
static float size;
static int status = -1;
static char* ctyp[7] = {" ", "pt", "line", "tri", "quad", "pyr", "prism"};

int main(int argc, char *argv[])
{
  int i, j, k, pts, pt;
  int matl = 0;
  FILE * outp;

  printf("make_spike_inp.c running \n"); 
  // read a .dat file 
  if(argc>1)
  {
    status = datread(argv[1], &data_points, &num_points,
                              &data_verts, &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;
  }
  outp = fopen("spike.inp", "w");
  fprintf(outp, "%3d %3d %3d %3d %3d \n", num_points, num_polys, 1, 0, 0); 
  for(i=0; i<num_points; i++)
  {
    fprintf(outp, "%3d %15f %15f %15f \n", i+1,
            data_points[i].x, data_points[i].y, data_points[i].z);
  }

  k = 0; // reading data_verts
  for(i=0; i<num_polys; i++)
  {
    pts = data_verts[k++];
    fprintf(outp, "%3d %3d %s ", i+1, matl, ctyp[pts]); // no cr
    for(j=0; j<pts; j++)
    {
      pt = data_verts[k++];
      fprintf(outp, "%d ", pt);
    }
    fprintf(outp, " \n");
  } // end write cells

  fprintf(outp, "%3d %3d \n", 1, 1);
  fprintf(outp, "Dirchilet boundary values \n");

  for(i=0; i<num_points; i++)
  {
    fprintf(outp, "%3d %15f \n", i+1, 
            sin(data_points[i].x+data_points[i].y+data_points[i].z+1.0));
  }
  fflush(outp);
  fclose(outp);
  printf("make_spike_inp.c wrote spike.inp \n"); 
  return 0;
} // end main of light_normal.c