/*

  ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
  º                                                                        º
  º                        Recursive Tetrahedrons                          º
  º                                                                        º
  º                       By Christopher D. Watkins                        º
  º                                                                        º
  ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ

*/

#include "stdio.h"
#include "stdlib.h"
#include "dos.h"
#include "math.h"
#include "defs.h"
#include "globals.h"
#include "mathb.h"
#include "graphb.h"

typedef char name[33];

FILE *textdiskfile;
name filename;

Word xresolution  = 1024;
Word yresolution  = 768;
Byte numberofframes  = 1;

void writeheader()
{
  fprintf(textdiskfile,"ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»\n");
  fprintf(textdiskfile,"º                                                                        º\n");
  fprintf(textdiskfile,"º                         Recursive Tetrahedrons                         º\n");
  fprintf(textdiskfile,"º                                                                        º\n");
  fprintf(textdiskfile,"º                       by Christopher D. Watkins                        º\n");
  fprintf(textdiskfile,"º                                                                        º\n");
  fprintf(textdiskfile,"ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"STATS\n");
  fprintf(textdiskfile,"   XRES    = %d\n",xresolution);
  fprintf(textdiskfile,"   YRES    = %d\n",yresolution);
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"FRAMES\n");
  fprintf(textdiskfile,"   NUMBER  = %d\n",numberofframes);
}

void writestudio()
{
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"ENVIRONMENT\n");
  fprintf(textdiskfile,"   LOCLWGT =    0.750    0.750    0.750\n");
  fprintf(textdiskfile,"   REFLWGT =    0.250    0.250    0.250\n");
  fprintf(textdiskfile,"   TRANWGT =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   MINWGT  =    0.030    0.030    0.030\n");
  fprintf(textdiskfile,"   MAXWGT  =    1.000    1.000    1.000\n");
  fprintf(textdiskfile,"   RDEPTH  =    3\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"LAMPS\n");
  fprintf(textdiskfile,"   REFLAMP = FALSE\n");
  fprintf(textdiskfile,"   LAMPREF =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   DISTEFF =    0.050\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"OBSERVER\n");
  fprintf(textdiskfile,"   FLENGTH =    3.500\n");
  fprintf(textdiskfile,"   OBSPOS  = -120.000  660.000  220.000\n");
  fprintf(textdiskfile,"   ROTATE  =  168.000\n");
  fprintf(textdiskfile,"   TILT    =   15.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"SKY\n");
  fprintf(textdiskfile,"   HORCOL  =    0.300    0.300    0.400\n");
  fprintf(textdiskfile,"   ZENCOL  =    0.700    0.800    1.000\n");
  fprintf(textdiskfile,"   CLOUDS  = FALSE\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = PLASTICTILE\n");
  fprintf(textdiskfile,"   TEXTURE = CHECKER\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.700    0.700    0.700\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.200    0.200    0.200\n");
  fprintf(textdiskfile,"   GLOSS   =    8.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"   TILE1   =    1.000    0.200    0.200\n");
  fprintf(textdiskfile,"   TILE2   =    0.250    0.250    0.250\n");
  fprintf(textdiskfile,"   TILE    =    0.012\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = CHROME\n");
  fprintf(textdiskfile,"   TEXTURE = SMOOTH\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.600    0.600    0.600\n");
  fprintf(textdiskfile,"   GLOSS   =   30.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = ITALIANMARBLE\n");
  fprintf(textdiskfile,"   TEXTURE = MARBLE\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.700    0.700    0.700\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.200    0.200    0.200\n");
  fprintf(textdiskfile,"   GLOSS   =   10.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = OAKWOOD\n");
  fprintf(textdiskfile,"   TEXTURE = WOOD\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.200    0.200    0.200\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.700    0.400    0.300\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   GLOSS   =    5.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = MIRROR\n");
  fprintf(textdiskfile,"   TEXTURE = SMOOTH\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.700    0.700    0.700\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.300    0.300    0.300\n");
  fprintf(textdiskfile,"   GLOSS   =   50.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"MATERIAL\n");
  fprintf(textdiskfile,"   TYPE    = BRASS\n");
  fprintf(textdiskfile,"   TEXTURE = SMOOTH\n");
  fprintf(textdiskfile,"   AMBRFL  =    0.100    0.100    0.100\n");
  fprintf(textdiskfile,"   DIFRFL  =    0.300    0.300    0.300\n");
  fprintf(textdiskfile,"   SPCRFL  =    0.600    0.500    0.250\n");
  fprintf(textdiskfile,"   GLOSS   =   20.000\n");
  fprintf(textdiskfile,"   TRANS   =    0.000    0.000    0.000\n");
  fprintf(textdiskfile,"   INDEX   =    1.000\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"GROUND\n");
  fprintf(textdiskfile,"   MATL    = PLASTICTILE\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"LAMP\n");
  fprintf(textdiskfile,"   LOC     =  350.000  150.000  300.000\n");
  fprintf(textdiskfile,"   RADIUS  =  100.000\n");
  fprintf(textdiskfile,"   INTENS  =    0.900    0.900    0.750\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"LAMP\n");
  fprintf(textdiskfile,"   LOC     = -350.000  250.000  300.000\n");
  fprintf(textdiskfile,"   RADIUS  =  100.000\n");
  fprintf(textdiskfile,"   INTENS  =    0.910    0.700    0.910\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"LAMP\n");
  fprintf(textdiskfile,"   LOC     =    0.000  -50.000  500.000\n");
  fprintf(textdiskfile,"   RADIUS  =  100.000\n");
  fprintf(textdiskfile,"   INTENS  =    0.810    0.800    0.910\n");
}

TDA v1, v2, avr, temp;

void guaranteeoutwardsurfacenormal(TDA pnt)
{
  TDA outray, nrm;

  VecSub(pnt,avr,outray);
  VecCross(v1,v2,nrm);
  if (!(VecDot(outray,nrm)>0.0))
  {
    VecCopy(v1,temp);
    VecCopy(v2,v1);
    VecCopy(temp,v2);
  }
}

Byte mtlnum;

void getmaterialnum()
{
  switch (mtlnum)
  {
    case 0: fprintf(textdiskfile,"   MATL    = CHROME\n");
	    break;
    case 1: fprintf(textdiskfile,"   MATL    = ITALIANMARBLE\n");
	    break;
    case 2: fprintf(textdiskfile,"   MATL    = OAKWOOD\n");
	    break;
    case 3: fprintf(textdiskfile,"   MATL    = MIRROR\n");
	    break;
    case 4: fprintf(textdiskfile,"   MATL    = BRASS\n");
	    break;
  }
}

Byte col;

void tetra(int x1, int y1, int z1,
	   int x2, int y2, int z2,
	   int x3, int y3, int z3,
	   int x4, int y4, int z4,
	   Byte level)
{
  TDA pnt1, pnt2, pnt3, pnt4;

  if ((level == 0))
  {
    mtlnum = random(5);
    Vec(x1,y1,z1,pnt1);
    Vec(x2,y2,z2,pnt2);
    Vec(x3,y3,z3,pnt3);
    Vec(x4,y4,z4,pnt4);
    VecAdd3(pnt1,pnt2,pnt3,temp);
    VecAdd(pnt4,temp,temp);
    VecScalMult(0.25,temp,avr);   /*  center of tetrahedron  */
    VecSub(pnt1,pnt3,v1);
    VecSub(pnt2,pnt3,v2);
    guaranteeoutwardsurfacenormal(pnt3);
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"TRIANGLE\n");
    fprintf(textdiskfile,"   LOC     = %3.3f %3.3f %3.3f\n",pnt3[0],pnt3[1],pnt3[2]);
    fprintf(textdiskfile,"   V1      = %3.3f %3.3f %3.3f\n",v1[0],v1[1],v1[2]);
    fprintf(textdiskfile,"   V2      = %3.3f %3.3f %3.3f\n",v2[0],v2[1],v2[2]);
    getmaterialnum();
    VecSub(pnt1,pnt2,v1);
    VecSub(pnt4,pnt2,v2);
    guaranteeoutwardsurfacenormal(pnt2);
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"TRIANGLE\n");
    fprintf(textdiskfile,"   LOC     = %3.3f %3.3f %3.3f\n",pnt2[0],pnt2[1],pnt2[2]);
    fprintf(textdiskfile,"   V1      = %3.3f %3.3f %3.3f\n",v1[0],v1[1],v1[2]);
    fprintf(textdiskfile,"   V2      = %3.3f %3.3f %3.3f\n",v2[0],v2[1],v2[2]);
    getmaterialnum();
    VecSub(pnt1,pnt4,v1);
    VecSub(pnt3,pnt4,v2);
    guaranteeoutwardsurfacenormal(pnt4);
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"TRIANGLE\n");
    fprintf(textdiskfile,"   LOC     = %3.3f %3.3f %3.3f\n",pnt4[0],pnt4[1],pnt4[2]);
    fprintf(textdiskfile,"   V1      = %3.3f %3.3f %3.3f\n",v1[0],v1[1],v1[2]);
    fprintf(textdiskfile,"   V2      = %3.3f %3.3f %3.3f\n",v2[0],v2[1],v2[2]);
    getmaterialnum();
    VecSub(pnt2,pnt3,v1);
    VecSub(pnt4,pnt3,v2);
    guaranteeoutwardsurfacenormal(pnt3);
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"\n");
    fprintf(textdiskfile,"TRIANGLE\n");
    fprintf(textdiskfile,"   LOC     = %3.3f %3.3f %3.3f\n",pnt3[0],pnt3[1],pnt3[2]);
    fprintf(textdiskfile,"   V1      = %3.3f %3.3f %3.3f\n",v1[0],v1[1],v1[2]);
    fprintf(textdiskfile,"   V2      = %3.3f %3.3f %3.3f\n",v2[0],v2[1],v2[2]);
    getmaterialnum();
    Draw_Line_3D(pnt1,pnt2,col);
    Draw_Line_3D(pnt2,pnt3,col);
    Draw_Line_3D(pnt3,pnt1,col);
    Draw_Line_3D(pnt1,pnt4,col);
    Draw_Line_3D(pnt4,pnt2,col);
    Draw_Line_3D(pnt3,pnt4,col);
  }
  else
  {
    tetra(x1,y1,z1,(x1 + x2) / 2,(y1 + y2) / 2,(z1 + z2) / 2,(x1 + x3) / 2,(y1 + y3) / 2,(z1 + z3) / 2,(x1 + x4) / 2,(y1 + y4) / 2,(z1 + z4) / 2,level - 1);
    tetra(x2,y2,z2,(x2 + x3) / 2,(y2 + y3) / 2,(z2 + z3) / 2,(x2 + x4) / 2,(y2 + y4) / 2,(z2 + z4) / 2,(x2 + x1) / 2,(y2 + y1) / 2,(z2 + z1) / 2,level - 1);
    tetra(x3,y3,z3,(x3 + x4) / 2,(y3 + y4) / 2,(z3 + z4) / 2,(x3 + x1) / 2,(y3 + y1) / 2,(z3 + z1) / 2,(x3 + x2) / 2,(y3 + y2) / 2,(z3 + z2) / 2,level - 1);
    tetra(x4,y4,z4,(x4 + x1) / 2,(y4 + y1) / 2,(z4 + z1) / 2,(x4 + x2) / 2,(y4 + y2) / 2,(z4 + z2) / 2,(x4 + x3) / 2,(y4 + y3) / 2,(z4 + z3) / 2,level - 1);
    col = (level % 6) * 36 + 35;
  }
}

/*
  ÉÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍ»
  º                        M a i n  P r o g r a m                          º
  ÈÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍÍͼ
*/

int x1, y1, z1;
int x2, y2, z2;
int x3, y3, z3;
int x4, y4, z4;

#define rot 0    /*  rotation around z-axis    0 = x-direction  */
#define ang1 0 + rot
#define ang2 120 + rot
#define ang3 240 + rot
#define scale 150    /*  scaling  */
#define xs scale
#define ys scale
#define zs scale
#define level 3    /*  level of recursion  */   /*  4^(Level+1) triangles  */
		       /*  0 ->   4 triangles  */
		       /*  1 ->  16 triangles  */
		       /*  2 ->  64 triangles  */
		       /*  3 -> 256 triangles  */


void main()
{
  randomize();
  strcpy(filename,"TETRA.RT");
  Title();
  clrscr();
  printf("Generating %s File\n",filename);
  printf("\n");
  printf("By Christopher D. Watkins\n");
  printf("\n");
  printf("\n");
  printf("Writing %d frames\n",numberofframes);
  printf("\n");
  printf("\n");
  printf("Depth of Recursion = %d\n",level);
  printf("\n");
  printf("This Depth yields %d triangles and %d tetrahedrons\n",IntPower(4,level + 1),IntPower(4,level));
  printf("\n");
  printf("\n");
  printf("< Hit any Key >\n");
  Wait_For_Key();
  textdiskfile=fopen(filename, "wt");
  Init_Graphics(45);
  Init_Perspective(false,0,0,500,500);
  Init_Plotting(20,10);
  Put_Axis_And_Palette(true);
  Axis_And_Palette();
  col = 35;
  x1 = Round(xs * CosD(ang1));
  y1 = Round(ys * SinD(ang1));
  z1 = 0;
  x2 = Round(xs * CosD(ang2));
  y2 = Round(ys * SinD(ang2));
  z2 = 0;
  x3 = Round(xs * CosD(ang3));
  y3 = Round(ys * SinD(ang3));
  z3 = 0;
  x4 = 0;
  y4 = 0;
  z4 = zs;
  writeheader();
  writestudio();
  tetra(x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4,level);
  fprintf(textdiskfile,"\n");   /*  extents for bounding tetrahedron  */
  fprintf(textdiskfile,"\n");   /*  BoundingTetraTest automatically turned on  */
  fprintf(textdiskfile,"TETRAEXTENTS\n");
  fprintf(textdiskfile,"   EXTENT1 = %d %d %d\n",x1,y1,z1);
  fprintf(textdiskfile,"   EXTENT2 = %d %d %d\n",x2,y2,z2);
  fprintf(textdiskfile,"   EXTENT3 = %d %d %d\n",x3,y3,z3);
  fprintf(textdiskfile,"   EXTENT4 = %d %d %d\n",x4,y4,z4);
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"BOUNDINGSPHERETEST\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"\n");
  fprintf(textdiskfile,"ENDFRAME\n");
  fclose(textdiskfile);
  Exit_Graphics();
}