Raytrace_dump.c

Go to the documentation of this file.

#include <setjmp.h>
#include "render.h"

#define UC unsigned short
#define US unsigned long

typedef long ivector2[2];

typedef struct traceNODE {
  long n,            
       mi[3],ma[3];  
  UC *vxi;           
  UC *oid;           
  US *fid;           
  struct traceNODE *octree[8]; 
  long id; struct traceNODE *last; long level; 
} node;

static void WriteVertices(void);
static void WriteEdges(void);
static void WriteFaces(void);
static void createVoxelData(node *n);

static HANDLE fq;
static FILE *fo;
static long newface,newedge,newvert,newglue,minX,minY,minZ,maxX,maxY,maxZ,voxelvid,voxelc,voxeleid,vertbias;
static char str[10];
static long FaceBase,EdgeBase,VertexBase;

static int Nvert=0,Nedge=0,Nface=0,Oid=0;
static ivector  *voxelv=NULL;
static ivector2 *voxele=NULL;
static ivector  *voxelf=NULL;


static void outsht(short sht){
 fputc((char)(sht >>  8),fo);
 fputc((char)(sht      ),fo);
}

static void outusht(unsigned short sht){
 fputc((char)(sht >>  8),fo);
 fputc((char)(sht      ),fo);
}

static void outlng(long lng){
 fputc((char)(lng >> 24),fo);
 fputc((char)(lng >> 16),fo);
 fputc((char)(lng >>  8),fo);
 fputc((char)(lng      ),fo);
}

void SaveDivided(long nvoxels, node *np, int o, char *Filename){  // for object o
  long pos1,pos2,pos3,pos4,pose;
  long lrulerx=1024, lrulery=1024, lrulerz=1024,
   rulerx=1024, rulery=1024, rulerz=1024;
  long i,FORMsize;
  HCURSOR hSave;
  char ruler_name[8]="CMS";
  short id;
  FORMsize=0; newvert=0; newface=0; newedge=0; newglue=0;
  if(nvoxels > 0){
    voxelvid=0; voxelc=0; voxeleid=0;
    voxelv=(ivector *)malloc(nvoxels*8*sizeof(ivector));
    if(voxelv == NULL)return;
    voxele=(ivector2 *)malloc(nvoxels*12*sizeof(ivector2));
    if(voxele == NULL)return;
    createVoxelData(np);
  }
  if(debug != NULL)fprintf(debug,"Saving [%s], expect %d voxels got %d voxels with %d extra vertices\n",
                           Filename, nvoxels,voxelc,voxelvid);
  if(debug != NULL)fprintf(debug,"got %d voxel edges\n",voxeleid);
  Oid=o;
  Nvert=Object[o].NoTraceVertices;
  Nface=Object[o].NoTraceFaces;
  Nedge=Nface*3;
  if( (fo=fopen(Filename,"wb")) == NULL){   /* now open for reading */
    return;
  }
  fprintf(fo,"FORM");
  pos1=ftell(fo); outlng(FORMsize);
  fprintf(fo,"OFXM");
  fprintf(fo,"AXIS"); outlng(48L); /* size of data chunk */
  outlng(0); outlng(0); outlng(0);
  for(i=0;i<9;i++)outlng(0);
  fprintf(fo,"OFFS"); outlng(16L);      /* follow point */
  outlng(3L);
  outlng(0);
  outlng(0);
  outlng(0);
  pos2 = -1; if(Nvert > 0){
    fprintf(fo,"VERT"); pos2=ftell(fo); outlng(0L); /* filled in later */
    WriteVertices();
    if(newvert == 0){ // no vertices - so reset the file pointer
      fseek(fo,pos2-4,0);
    }
  }
  pos3 = -1; if(Nedge > 0 && newvert > 1){
      fprintf(fo,"EDGE"); pos3=ftell(fo); outlng(0L); WriteEdges();
  }
  pos4 = -1; if(Nface > 0 && newvert > 2){
    fprintf(fo,"NFCE"); pos4=ftell(fo); outlng(0L); WriteFaces();
    if(newface == 0)fseek(fo,pos4-4,0);
  }
  fprintf(fo,"UNIT"); outlng((long)(7*sizeof(long)+8));
  outlng((long)1024);
  fwrite(ruler_name,1,8,fo);
  outlng(lrulerx); outlng(lrulery); outlng(lrulerz);
  outlng(rulerx);  outlng(rulery);  outlng(rulerz);
  fprintf(fo,"GRID"); outlng((long)(2*sizeof(long)));
  outlng(4096); outlng(4096); /* second is unused */
  pose=ftell(fo);    /* the end of the object */
  // FIX UP the FILE lengths
  if(pos2 > 0){fseek(fo,pos2,0); outlng(12*newvert); }  // vertex chunk
  if(pos3 > 0){fseek(fo,pos3,0); // size of edge chunk 
    outlng(8*newedge);
  }
  if(pos4   > 0){fseek(fo,pos4,0);   outlng(20*newface);}     // size of face chunk
  fseek(fo,pos1,0); outlng(pose-8);
  fclose(fo);
  if(voxelv != NULL)free(voxelv);
  if(voxele != NULL)free(voxele);
  if(voxelf != NULL)free(voxelf);
  if(debug != NULL)fprintf(debug,"Model with %d verts %d edges %d faces\n",newvert,newedge,newface);
  return;
}

static void WriteVertices(void){
 long i,j;
 long x,y,z;
 for(j=0;j<Nvert;j++){
  x=Object[Oid].Vtracebase[j].ip[0];
  y=Object[Oid].Vtracebase[j].ip[1];
  z=Object[Oid].Vtracebase[j].ip[2];
  outlng(x); outlng(y); outlng(z);
  newvert++;
 }
 if(voxelvid > 0){
   vertbias=newvert;
   for(j=0;j<voxelvid;j++){
     outlng(voxelv[j][0]); outlng(voxelv[j][1]); outlng(voxelv[j][2]);
     newvert++;
   }
 }
}

static void WriteEdges(){
 long i,id1,id2,id3;
 char color[3]={255,255,255};
 char bSmooth=0; 
 for(i=0;i<Nface;i++){
   id1=Object[Oid].Ftracebase[i].V[0];
   id2=Object[Oid].Ftracebase[i].V[1];
   id3=Object[Oid].Ftracebase[i].V[2];
   outlng(id1); outlng(id2);
   newedge++;
   outlng(id2); outlng(id3);
   newedge++;
   outlng(id3); outlng(id1);
   newedge++;
 }
 if(voxeleid > 0){
   for(i=0;i<voxeleid;i++){
     outlng(voxele[i][0]+vertbias); outlng(voxele[i][1]+vertbias);
     newedge++;
   }
 }
 
}

static void WriteFaces(void){
 long i,id1,id2,id3,vx,vx1;
 unsigned char color[3]={255,255,255};
 char bSmooth=0; 
 for(i=0;i<Nface;i++){
   id1=Object[Oid].Ftracebase[i].V[0];
   id2=Object[Oid].Ftracebase[i].V[1];
   id3=Object[Oid].Ftracebase[i].V[2];
   vx=Object[Oid].Ftracebase[i].x;
   vx1=vx%512;
   color[0] = vx1/2;
   vx1=vx%64;
   color[1] = vx1*4;
   vx=vx*8; 
   color[2] = vx1*32;
   outlng(id1);
   outlng(id2);
   outlng(id3);
   fwrite(color,1,3,fo);
   fwrite(&bSmooth,1,1,fo);
   outsht(-1);
   outsht(-1);
   newface++;
 }
}

static void createVoxelData(node *n){
 voxele[voxeleid][0]=voxelvid; voxele[voxeleid][1]=voxelvid+1; voxeleid++;
 voxele[voxeleid][0]=voxelvid+1; voxele[voxeleid][1]=voxelvid+2; voxeleid++;
 voxele[voxeleid][0]=voxelvid+2; voxele[voxeleid][1]=voxelvid+3; voxeleid++;
 voxele[voxeleid][0]=voxelvid+3; voxele[voxeleid][1]=voxelvid; voxeleid++;
 voxele[voxeleid][0]=voxelvid+4; voxele[voxeleid][1]=voxelvid+5; voxeleid++;
 voxele[voxeleid][0]=voxelvid+5; voxele[voxeleid][1]=voxelvid+6; voxeleid++;
 voxele[voxeleid][0]=voxelvid+6; voxele[voxeleid][1]=voxelvid+7; voxeleid++;
 voxele[voxeleid][0]=voxelvid+7; voxele[voxeleid][1]=voxelvid+4; voxeleid++;
 voxele[voxeleid][0]=voxelvid; voxele[voxeleid][1]=voxelvid+4; voxeleid++;
 voxele[voxeleid][0]=voxelvid+1; voxele[voxeleid][1]=voxelvid+5; voxeleid++;
 voxele[voxeleid][0]=voxelvid+2; voxele[voxeleid][1]=voxelvid+6; voxeleid++;
 voxele[voxeleid][0]=voxelvid+3; voxele[voxeleid][1]=voxelvid+7; voxeleid++;

 voxelv[voxelvid][0]=n->mi[0];  voxelv[voxelvid][1]=n->mi[1];  voxelv[voxelvid][2]=n->mi[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->ma[0];  voxelv[voxelvid][1]=n->mi[1];  voxelv[voxelvid][2]=n->mi[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->ma[0];  voxelv[voxelvid][1]=n->ma[1];  voxelv[voxelvid][2]=n->mi[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->mi[0];  voxelv[voxelvid][1]=n->ma[1];  voxelv[voxelvid][2]=n->mi[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->mi[0];  voxelv[voxelvid][1]=n->mi[1];  voxelv[voxelvid][2]=n->ma[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->ma[0];  voxelv[voxelvid][1]=n->mi[1];  voxelv[voxelvid][2]=n->ma[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->ma[0];  voxelv[voxelvid][1]=n->ma[1];  voxelv[voxelvid][2]=n->ma[2];  voxelvid++; 
 voxelv[voxelvid][0]=n->mi[0];  voxelv[voxelvid][1]=n->ma[1];  voxelv[voxelvid][2]=n->ma[2];  voxelvid++; 
 voxelc++;
 if(n->octree[0] != NULL)createVoxelData(n->octree[0]);
 if(n->octree[1] != NULL)createVoxelData(n->octree[1]);
 if(n->octree[2] != NULL)createVoxelData(n->octree[2]);
 if(n->octree[3] != NULL)createVoxelData(n->octree[3]);
 if(n->octree[4] != NULL)createVoxelData(n->octree[4]);
 if(n->octree[5] != NULL)createVoxelData(n->octree[5]);
 if(n->octree[6] != NULL)createVoxelData(n->octree[6]);
 if(n->octree[7] != NULL)createVoxelData(n->octree[7]);
 return;
}