--- ray/src/cv/mgflib/parser.h 1994/06/22 15:33:42 1.2 +++ ray/src/cv/mgflib/parser.h 1995/05/11 21:48:42 1.25 @@ -1,4 +1,4 @@ -/* Copyright (c) 1994 Regents of the University of California */ +/* Copyright (c) 1995 Regents of the University of California */ /* SCCSid "$SunId$ LBL" */ @@ -8,47 +8,65 @@ /* must include stdio.h before us */ - /* Entities (order doesn't really matter) */ -#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_VMAJOR 1 /* major version number */ +#define MG_VMINOR 0 /* minor version number */ -#define MG_NENTITIES 22 + /* Entities (list is only appended, never modified) */ +#define MG_E_COMMENT 0 /* # */ +#define MG_E_COLOR 1 /* c */ +#define MG_E_CCT 2 /* cct */ +#define MG_E_CONE 3 /* cone */ +#define MG_E_CMIX 4 /* cmix */ +#define MG_E_CSPEC 5 /* cspec */ +#define MG_E_CXY 6 /* cxy */ +#define MG_E_CYL 7 /* cyl */ +#define MG_E_ED 8 /* ed */ +#define MG_E_FACE 9 /* f */ +#define MG_E_INCLUDE 10 /* i */ +#define MG_E_IES 11 /* ies */ +#define MG_E_IR 12 /* ir */ +#define MG_E_MATERIAL 13 /* m */ +#define MG_E_NORMAL 14 /* n */ +#define MG_E_OBJECT 15 /* o */ +#define MG_E_POINT 16 /* p */ +#define MG_E_PRISM 17 /* prism */ +#define MG_E_RD 18 /* rd */ +#define MG_E_RING 19 /* ring */ +#define MG_E_RS 20 /* rs */ +#define MG_E_SIDES 21 /* sides */ +#define MG_E_SPH 22 /* sph */ +#define MG_E_TD 23 /* td */ +#define MG_E_TORUS 24 /* torus */ +#define MG_E_TS 25 /* ts */ +#define MG_E_VERTEX 26 /* v */ +#define MG_E_XF 27 /* xf */ + /* end of Version 1 entities */ -#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_NENTITIES 28 /* total # entities */ +#define MG_NELIST {28} /* entity count for version 1 and up */ + +#define MG_NAMELIST {"#","c","cct","cone","cmix","cspec","cxy","cyl","ed",\ + "f","i","ies","ir","m","n","o","p","prism","rd",\ + "ring","rs","sides","sph","td","torus","ts","v","xf"} + #define MG_MAXELEN 6 extern char mg_ename[MG_NENTITIES][MG_MAXELEN]; - /* Handler routines for each entity */ - + /* Handler routines for each entity and unknown ones */ #ifdef NOPROTO extern int (*mg_ehand[MG_NENTITIES])(); +extern int (*mg_uhand)(); +extern int mg_defuhand(); #else extern int (*mg_ehand[MG_NENTITIES])(int argc, char **argv); +extern int (*mg_uhand)(int argc, char **argv); +extern int mg_defuhand(int, char **); #endif +extern unsigned mg_nunknown; /* count of unknown entities */ + /* Error codes */ #define MG_OK 0 /* normal return value */ #define MG_EUNK 1 /* unknown entity */ @@ -64,7 +82,7 @@ extern int (*mg_ehand[MG_NENTITIES])(int argc, char ** #define MG_NERRS 11 -extern char *mg_err[MG_NERRS]; +extern char *mg_err[MG_NERRS]; /* list of error messages */ /* * The general process for running the parser is to fill in the mg_ehand @@ -75,24 +93,34 @@ extern char *mg_err[MG_NERRS]; * To read from standard input, use NULL as the file name. * For additional control over error reporting and file management, * use mg_open, mg_read, mg_parse and mg_close instead of mg_load. + * To pass an entity of your own construction to the parser, use + * the mg_handle function rather than the mg_ehand routines directly. + * (The first argument to mg_handle is the entity #, or -1.) * To free any data structures and clear the parser, use mg_clear. - * If there is an error, mg_load, mg_open, mg_parse, and mg_rewind - * will return an error from the list above. In addition, mg_load - * will report the error to stderr. The mg_read routine returns 0 - * when the end of file has been reached. + * If there is an error, mg_load, mg_open, mg_parse, mg_handle and + * mg_fgoto will return an error from the list above. In addition, + * mg_load will report the error to stderr. The mg_read routine + * returns 0 when the end of file has been reached. */ -#define MG_MAXLINE 512 /* maximum input line length */ +#define MG_MAXLINE 4096 /* maximum input line length */ #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 */ + int fid; /* unique file context id */ char inpline[MG_MAXLINE]; /* input line */ int lineno; /* line number */ struct mg_fctxt *prev; /* previous context */ } MG_FCTXT; +typedef struct { + int fid; /* file this position is for */ + int lineno; /* line number in file */ + long offset; /* offset from beginning */ +} MG_FPOS; + extern MG_FCTXT *mg_file; /* current file context */ #ifdef NOPROTO @@ -101,20 +129,22 @@ extern int mg_load(); /* parse a file */ extern int mg_open(); /* open new input file */ extern int mg_read(); /* read next line */ extern int mg_parse(); /* parse current line */ -extern int mg_rewind(); /* rewind input file */ +extern void mg_fgetpos(); /* get position on input file */ +extern int mg_fgoto(); /* go to position on input file */ extern void mg_close(); /* close input file */ extern void mg_clear(); /* clear parser */ -extern int mg_iterate(); +extern int mg_handle(); /* handle an entity */ #else extern void mg_init(void); /* fill in mg_ehand array */ extern int mg_load(char *); /* parse a file */ extern int mg_open(MG_FCTXT *, char *); /* open new input file */ extern int mg_read(void); /* read next line */ extern int mg_parse(void); /* parse current line */ -extern int mg_rewind(void); /* rewind input file */ +extern void mg_fgetpos(MG_FPOS *); /* get position on input file */ +extern int mg_fgoto(MG_FPOS *); /* go to position on input file */ extern void mg_close(void); /* close input file */ extern void mg_clear(void); /* clear parser */ -extern int mg_iterate(int, char **, int (*)(void)); +extern int mg_handle(int, int, char **); /* handle an entity */ #endif #ifndef MG_NQCD @@ -131,24 +161,39 @@ extern int mg_nqcdivs; /* divisions per quarter circl extern int mg_entity(); /* get entity number from its name */ extern int isint(); /* non-zero if integer format */ extern int isflt(); /* non-zero if floating point format */ +extern int isname(); /* non-zero if legal identifier name */ #else extern int mg_entity(char *); /* get entity number from its name */ extern int isint(char *); /* non-zero if integer format */ extern int isflt(char *); /* non-zero if floating point format */ +extern int isname(char *); /* non-zero if legal identifier name */ #endif /************************************************************************ * Definitions for 3-d vector manipulation functions */ -#ifdef SMLFLT -#define FLOAT float +#ifdef SMLFLT +#define FLOAT float +#define FTINY (1e-3) #else -#define FLOAT double +#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 @@ -160,13 +205,46 @@ 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_CLPWM (683./C_CMAXV) /* peak lumens/watt multiplier */ + +#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 */ +#define C_CSEFF 020 /* flag if efficacy set */ + typedef struct { - float cx, cy; /* XY chromaticity coordinates */ -} C_COLOR; /* color context */ + int clock; /* incremented each change */ + 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 */ + float eff; /* efficacy (lumens/watt) */ +} C_COLOR; +#define C_DEFCOLOR { 1, C_CDXY|C_CSXY|C_CSSPEC|C_CSEFF,\ + {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., 178.006 } + +#define c_cval(c,l) ((double)(c)->ssamp[((l)-C_MINWL)/C_CWLI] / (c)->ssum) + +#define C_1SIDEDTHICK 0.005 /* assumed thickness of 1-sided mat. */ + typedef struct { - char *name; /* material name */ int clock; /* incremented each change -- resettable */ + int sided; /* 1 if surface is 1-sided, 0 for 2-sided */ + float nr, ni; /* index of refraction, real and imaginary */ float rd; /* diffuse reflectance */ C_COLOR rd_c; /* diffuse reflectance color */ float td; /* diffuse transmittance */ @@ -182,30 +260,41 @@ typedef struct { } C_MATERIAL; /* material context */ typedef struct { + int clock; /* incremented each change -- resettable */ FVECT p, n; /* point and normal */ } C_VERTEX; /* vertex context */ -#define C_DEFCOLOR {.333,.333} -#define C_DEFMATERIAL {NULL,1,0.,C_DEFCOLOR,0.,C_DEFCOLOR,0.,C_DEFCOLOR,\ +#define C_DEFMATERIAL {1,0,1.,0.,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.}} +#define C_DEFVERTEX {1,{0.,0.,0.},{0.,0.,0.}} extern C_COLOR *c_ccolor; /* the current color */ +extern char *c_ccname; /* current color name */ extern C_MATERIAL *c_cmaterial; /* the current material */ +extern char *c_cmname; /* current material name */ extern C_VERTEX *c_cvertex; /* the current vertex */ +extern char *c_cvname; /* current vertex name */ #ifdef NOPROTO extern int c_hcolor(); /* handle color entity */ extern int c_hmaterial(); /* handle material entity */ extern int c_hvertex(); /* handle vertex entity */ extern void c_clearall(); /* clear context tables */ +extern C_MATERIAL *c_getmaterial(); /* get a named material */ 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_MATERIAL *c_getmaterial(char *); /* get a named material */ 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 /************************************************************************* @@ -233,7 +322,6 @@ typedef FLOAT MAT4[4][4]; #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.}, \ @@ -245,23 +333,41 @@ extern MAT4 m4ident; /* regular transformation */ typedef struct { - MAT4 xfm; /* transform matrix */ - FLOAT sca; /* scalefactor */ + MAT4 xfm; /* transform matrix */ + FLOAT sca; /* scalefactor */ } XF; #define identxf(xp) (void)(setident4((xp)->xfm),(xp)->sca=1.0) +#define XF_MAXDIM 8 /* maximum array dimensions */ + +struct xf_array { + MG_FPOS spos; /* starting position on input */ + int ndim; /* number of array dimensions */ + struct { + short i, n; /* current count and maximum */ + char arg[8]; /* string argument value */ + } aarg[XF_MAXDIM]; +}; + typedef struct xf_spec { - short xac; /* transform argument count */ - short xav0; /* zeroeth argument in xf_argv array */ - XF xf; /* cumulative transformation */ - struct xf_spec *prev; /* previous transformation context */ -} XF_SPEC; + long xid; /* unique transform id */ + short xac; /* context argument count */ + short rev; /* boolean true if vertices reversed */ + XF xf; /* cumulative transformation */ + struct xf_array *xarr; /* transformation array pointer */ + struct xf_spec *prev; /* previous transformation context */ +} XF_SPEC; /* followed by argument buffer */ -extern int xf_argc; /* total # transform args. */ -extern char **xf_argv; /* transform arguments */ -extern XF_SPEC *xf_context; /* current context */ +extern XF_SPEC *xf_context; /* current transform context */ +extern char **xf_argend; /* last transform argument */ +#define xf_ac(xf) ((xf)->xac) +#define xf_av(xf) (xf_argend - (xf)->xac) + +#define xf_argc (xf_context==NULL ? 0 : xf_ac(xf_context)) +#define xf_argv xf_av(xf_context) + /* * The transformation handler should do most of the work that needs * doing. Just pass it any xf entities, then use the associated @@ -304,7 +410,7 @@ extern void xf_clear(void); /* clear xf stack */ 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 @@ -329,3 +435,4 @@ extern double atof(); extern MEM_PTR malloc(); extern MEM_PTR calloc(); extern MEM_PTR realloc(); +extern void free();