--- ray/src/cv/mgflib/parser.h 1994/06/21 14:45:41 1.1 +++ ray/src/cv/mgflib/parser.h 1994/06/25 09:47:58 1.7 @@ -12,30 +12,34 @@ #define MG_E_COMMENT 0 #define MG_E_COLOR 1 #define MG_E_CONE 2 -#define MG_E_CXY 3 -#define MG_E_CYL 4 -#define MG_E_ED 5 -#define MG_E_FACE 6 -#define MG_E_INCLUDE 7 -#define MG_E_IES 8 -#define MG_E_MATERIAL 9 -#define MG_E_NORMAL 10 -#define MG_E_OBJECT 11 -#define MG_E_POINT 12 -#define MG_E_RD 13 -#define MG_E_RING 14 -#define MG_E_RS 15 -#define MG_E_SPH 16 -#define MG_E_TD 17 -#define MG_E_TORUS 18 -#define MG_E_TS 19 -#define MG_E_VERTEX 20 -#define MG_E_XF 21 +#define MG_E_CMIX 3 +#define MG_E_CSPEC 4 +#define MG_E_CXY 5 +#define MG_E_CYL 6 +#define MG_E_ED 7 +#define MG_E_FACE 8 +#define MG_E_INCLUDE 9 +#define MG_E_IES 10 +#define MG_E_MATERIAL 11 +#define MG_E_NORMAL 12 +#define MG_E_OBJECT 13 +#define MG_E_POINT 14 +#define MG_E_PRISM 15 +#define MG_E_RD 16 +#define MG_E_RING 17 +#define MG_E_RS 18 +#define MG_E_SPH 19 +#define MG_E_TD 20 +#define MG_E_TORUS 21 +#define MG_E_TS 22 +#define MG_E_VERTEX 23 +#define MG_E_XF 24 -#define MG_NENTITIES 22 +#define MG_NENTITIES 25 -#define MG_NAMELIST {"#","c","cone","cxy","cyl","ed","f","i","ies",\ - "m","n","o","p","rd","ring","rs","sph","td","torus","ts","v","xf"} +#define MG_NAMELIST {"#","c","cone","cmix","cspec","cxy","cyl","ed","f",\ + "i","ies","m","n","o","p","prism","rd","ring","rs",\ + "sph","td","torus","ts","v","xf"} #define MG_MAXELEN 6 @@ -60,8 +64,9 @@ extern int (*mg_ehand[MG_NENTITIES])(int argc, char ** #define MG_EINCL 7 /* error in included file */ #define MG_EMEM 8 /* out of memory */ #define MG_ESEEK 9 /* file seek error */ +#define MG_EBADMAT 10 /* bad material specification */ -#define MG_NERRS 10 +#define MG_NERRS 11 extern char *mg_err[MG_NERRS]; @@ -85,7 +90,7 @@ extern char *mg_err[MG_NERRS]; #define MG_MAXARGC (MG_MAXLINE/4) /* maximum argument count */ typedef struct mg_fctxt { - char *fname; /* file name */ + char fname[96]; /* file name */ FILE *fp; /* stream pointer */ char inpline[MG_MAXLINE]; /* input line */ int lineno; /* line number */ @@ -140,8 +145,27 @@ extern int isflt(char *); /* non-zero if floating poi * Definitions for 3-d vector manipulation functions */ -typedef double FVECT[3]; +#ifdef SMLFLT +#define FLOAT float +#define FTINY (1e-3) +#else +#define FLOAT double +#define FTINY (1e-6) +#endif +#define FHUGE (1e10) +typedef FLOAT FVECT[3]; + +#define VCOPY(v1,v2) ((v1)[0]=(v2)[0],(v1)[1]=(v2)[1],(v1)[2]=(v2)[2]) +#define DOT(v1,v2) ((v1)[0]*(v2)[0]+(v1)[1]*(v2)[1]+(v1)[2]*(v2)[2]) +#define VSUM(vr,v1,v2,f) ((vr)[0]=(v1)[0]+(f)*(v2)[0], \ + (vr)[1]=(v1)[1]+(f)*(v2)[1], \ + (vr)[2]=(v1)[2]+(f)*(v2)[2]) + +#define is0vect(v) (DOT(v,v) <= FTINY*FTINY) + +#define round0(x) if (x <= FTINY && x >= -FTINY) x = 0 + #ifdef NOPROTO extern double normalize(); /* normalize a vector */ #else @@ -153,32 +177,77 @@ extern double normalize(FVECT); /* normalize a vector * (materials, colors, vectors) */ +#define C_CMINWL 380 /* minimum wavelength */ +#define C_CMAXWL 780 /* maximum wavelength */ +#define C_CNSS 41 /* number of spectral samples */ +#define C_CWLI ((C_CMAXWL-C_CMINWL)/(C_CNSS-1)) +#define C_CMAXV 10000 /* nominal maximum sample value */ + +#define C_CSSPEC 01 /* flag if spectrum is set */ +#define C_CDSPEC 02 /* flag if defined w/ spectrum */ +#define C_CSXY 04 /* flag if xy is set */ +#define C_CDXY 010 /* flag if defined w/ xy */ + typedef struct { - double cx, cy; /* XY chromaticity coordinates */ -} C_COLOR; /* color context */ + short flags; /* what's been set */ + short ssamp[C_CNSS]; /* spectral samples, min wl to max */ + long ssum; /* straight sum of spectral values */ + float cx, cy; /* xy chromaticity value */ +} C_COLOR; +#define C_DEFCOLOR { C_CDXY|C_CSXY|C_CSSPEC,\ + {C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\ + C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV},\ + (long)C_CNSS*C_CMAXV, 1./3., 1./3. } + +#define C_CIEX { C_CDSPEC|C_CSSPEC|C_CSXY,\ + {14,42,143,435,1344,2839,3483,3362,2908,1954,956,\ + 320,49,93,633,1655,2904,4334,5945,7621,9163,10263,\ + 10622,10026,8544,6424,4479,2835,1649,874,468,227,\ + 114,58,29,14,7,3,2,1,0}, 106836L, .735, .265 } + +#define C_CIEY { C_CDSPEC|C_CSSPEC|C_CSXY,\ + {0,1,4,12,40,116,230,380,600,910,1390,2080,3230,\ + 5030,7100,8620,9540,9950,9950,9520,8700,7570,6310,\ + 5030,3810,2650,1750,1070,610,320,170,82,41,21,10,\ + 5,2,1,1,0,0}, 106856L, .274, .717 } + +#define C_CIEZ { C_CDSPEC|C_CSSPEC|C_CSXY,\ + {65,201,679,2074,6456,13856,17471,17721,16692,\ + 12876,8130,4652,2720,1582,782,422,203,87,39,21,17,\ + 11,8,3,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},\ + 106770L, .167, .009 } + +#define c_cval(c,l) ((double)(c)->ssamp[((l)-C_MINWL)/C_CWLI] / (c)->sum) + typedef struct { - double rd; /* diffuse reflectance */ + char *name; /* material name */ + int clock; /* incremented each change -- resettable */ + float rd; /* diffuse reflectance */ C_COLOR rd_c; /* diffuse reflectance color */ - double td; /* diffuse transmittance */ + float td; /* diffuse transmittance */ C_COLOR td_c; /* diffuse transmittance color */ - double ed; /* diffuse emittance */ + float ed; /* diffuse emittance */ C_COLOR ed_c; /* diffuse emittance color */ - double rs; /* specular reflectance */ + float rs; /* specular reflectance */ C_COLOR rs_c; /* specular reflectance color */ - double rs_a; /* specular reflectance roughness */ - double ts; /* specular transmittance */ + float rs_a; /* specular reflectance roughness */ + float ts; /* specular transmittance */ C_COLOR ts_c; /* specular transmittance color */ - double ts_a; /* specular transmittance roughness */ + float ts_a; /* specular transmittance roughness */ } C_MATERIAL; /* material context */ typedef struct { FVECT p, n; /* point and normal */ } C_VERTEX; /* vertex context */ -#define C_DEFCOLOR {.333,.333} -#define C_DEFMATERIAL {0.,C_DEFCOLOR,0.,C_DEFCOLOR,0.,C_DEFCOLOR,\ - 0.,C_DEFCOLOR,0.,0.,C_DEFCOLOR,0.} +#define C_DEFMATERIAL {NULL,1,0.,C_DEFCOLOR,0.,C_DEFCOLOR,0.,C_DEFCOLOR,\ + 0.,C_DEFCOLOR,0.,0.,C_DEFCOLOR,0.} #define C_DEFVERTEX {{0.,0.,0.},{0.,0.,0.}} extern C_COLOR *c_ccolor; /* the current color */ @@ -191,12 +260,18 @@ extern int c_hmaterial(); /* handle material entity extern int c_hvertex(); /* handle vertex entity */ extern void c_clearall(); /* clear context tables */ extern C_VERTEX *c_getvert(); /* get a named vertex */ +extern C_COLOR *c_getcolor(); /* get a named color */ +extern void c_ccvt(); /* fix color representation */ +extern int c_isgrey(); /* check if color is grey */ #else extern int c_hcolor(int, char **); /* handle color entity */ extern int c_hmaterial(int, char **); /* handle material entity */ extern int c_hvertex(int, char **); /* handle vertex entity */ extern void c_clearall(void); /* clear context tables */ extern C_VERTEX *c_getvert(char *); /* get a named vertex */ +extern C_COLOR *c_getcolor(char *); /* get a named color */ +extern void c_ccvt(C_COLOR *, int); /* fix color representation */ +extern int c_isgrey(C_COLOR *); /* check if color is grey */ #endif /************************************************************************* @@ -208,21 +283,22 @@ extern char **obj_name; /* names in hierarchy */ #ifdef NOPROTO extern int obj_handler(); /* handle an object entity */ +extern void obj_clear(); /* clear object stack */ #else extern int obj_handler(int, char **); /* handle an object entity */ +extern void obj_clear(void); /* clear object stack */ #endif /************************************************************************** * Definitions for hierarchical transformation handler */ -typedef double MAT4[4][4]; +typedef FLOAT MAT4[4][4]; #ifdef BSD #define copymat4(m4a,m4b) bcopy((char *)m4b,(char *)m4a,sizeof(MAT4)) #else #define copymat4(m4a,m4b) (void)memcpy((char *)m4a,(char *)m4b,sizeof(MAT4)) -extern char *memcpy(); #endif #define MAT4IDENT { {1.,0.,0.,0.}, {0.,1.,0.,0.}, \ @@ -235,7 +311,7 @@ extern MAT4 m4ident; /* regular transformation */ typedef struct { MAT4 xfm; /* transform matrix */ - double sca; /* scalefactor */ + FLOAT sca; /* scalefactor */ } XF; #define identxf(xp) (void)(setident4((xp)->xfm),(xp)->sca=1.0) @@ -270,6 +346,7 @@ extern void xf_xfmpoint(); /* transform point */ extern void xf_xfmvect(); /* transform vector */ extern void xf_rotvect(); /* rotate vector */ extern double xf_scale(); /* scale a value */ +extern void xf_clear(); /* clear xf stack */ /* The following are support routines you probably won't call directly */ @@ -280,18 +357,19 @@ extern int xf(); /* interpret transform spec. */ #else -extern int xf_handler(); /* handle xf entity */ -extern void xf_xfmpoint(); /* transform point */ -extern void xf_xfmvect(); /* transform vector */ -extern void xf_rotvect(); /* rotate vector */ -extern double xf_scale(); /* scale a value */ +extern int xf_handler(int, char **); /* handle xf entity */ +extern void xf_xfmpoint(FVECT, FVECT); /* transform point */ +extern void xf_xfmvect(FVECT, FVECT); /* transform vector */ +extern void xf_rotvect(FVECT, FVECT); /* rotate vector */ +extern double xf_scale(double); /* scale a value */ +extern void xf_clear(void); /* clear xf stack */ /* The following are support routines you probably won't call directly */ extern void multmat4(MAT4, MAT4, MAT4); /* m4a = m4b X m4c */ extern void multv3(FVECT, FVECT, MAT4); /* v3a = v3b X m4 (vectors) */ extern void multp3(FVECT, FVECT, MAT4); /* p3a = p3b X m4 (points) */ -extern int xf(XF, int, char **); /* interpret transform spec. */ +extern int xf(XF *, int, char **); /* interpret transform spec. */ #endif @@ -316,3 +394,4 @@ extern double atof(); extern MEM_PTR malloc(); extern MEM_PTR calloc(); extern MEM_PTR realloc(); +extern void free();