ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/mgflib/parser.h
Revision: 1.17
Committed: Fri Feb 24 16:50:05 1995 UTC (29 years, 2 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 1.16: +7 -7 lines
Log Message:
fixed context handlers for spectral power distribution data

File Contents

# Content
1 /* Copyright (c) 1994 Regents of the University of California */
2
3 /* SCCSid "$SunId$ LBL" */
4
5 /*
6 * Header file for MGF interpreter
7 */
8
9 /* must include stdio.h before us */
10
11 /* Entities (order doesn't really matter) */
12 #define MG_E_COMMENT 0
13 #define MG_E_COLOR 1
14 #define MG_E_CONE 2
15 #define MG_E_CMIX 3
16 #define MG_E_CSPEC 4
17 #define MG_E_CXY 5
18 #define MG_E_CYL 6
19 #define MG_E_ED 7
20 #define MG_E_FACE 8
21 #define MG_E_INCLUDE 9
22 #define MG_E_IES 10
23 #define MG_E_MATERIAL 11
24 #define MG_E_NORMAL 12
25 #define MG_E_OBJECT 13
26 #define MG_E_POINT 14
27 #define MG_E_PRISM 15
28 #define MG_E_RD 16
29 #define MG_E_RING 17
30 #define MG_E_RS 18
31 #define MG_E_SIDES 19
32 #define MG_E_SPH 20
33 #define MG_E_TD 21
34 #define MG_E_TORUS 22
35 #define MG_E_TS 23
36 #define MG_E_VERTEX 24
37 #define MG_E_XF 25
38
39 #define MG_NENTITIES 26
40
41 #define MG_NAMELIST {"#","c","cone","cmix","cspec","cxy","cyl","ed","f",\
42 "i","ies","m","n","o","p","prism","rd","ring","rs",\
43 "sides","sph","td","torus","ts","v","xf"}
44
45 #define MG_MAXELEN 6
46
47 extern char mg_ename[MG_NENTITIES][MG_MAXELEN];
48
49 /* Handler routines for each entity */
50
51 #ifdef NOPROTO
52 extern int (*mg_ehand[MG_NENTITIES])();
53 #else
54 extern int (*mg_ehand[MG_NENTITIES])(int argc, char **argv);
55 #endif
56
57 /* Error codes */
58 #define MG_OK 0 /* normal return value */
59 #define MG_EUNK 1 /* unknown entity */
60 #define MG_EARGC 2 /* wrong number of arguments */
61 #define MG_ETYPE 3 /* argument type error */
62 #define MG_EILL 4 /* illegal argument value */
63 #define MG_EUNDEF 5 /* undefined reference */
64 #define MG_ENOFILE 6 /* cannot open input file */
65 #define MG_EINCL 7 /* error in included file */
66 #define MG_EMEM 8 /* out of memory */
67 #define MG_ESEEK 9 /* file seek error */
68 #define MG_EBADMAT 10 /* bad material specification */
69
70 #define MG_NERRS 11
71
72 extern char *mg_err[MG_NERRS];
73
74 /*
75 * The general process for running the parser is to fill in the mg_ehand
76 * array with handlers for each entity you know how to handle.
77 * Then, call mg_init to fill in the rest. This function will report
78 * an error and quit if you try to support an inconsistent set of entities.
79 * For each file you want to parse, call mg_load with the file name.
80 * To read from standard input, use NULL as the file name.
81 * For additional control over error reporting and file management,
82 * use mg_open, mg_read, mg_parse and mg_close instead of mg_load.
83 * To pass an entity of your own construction to the parser, use
84 * the mg_handle function rather than the mg_ehand routines directly.
85 * (The first argument to mg_handle is the entity #, or -1.)
86 * To free any data structures and clear the parser, use mg_clear.
87 * If there is an error, mg_load, mg_open, mg_parse, mg_handle and
88 * mg_rewind will return an error from the list above. In addition,
89 * mg_load will report the error to stderr. The mg_read routine
90 * returns 0 when the end of file has been reached.
91 */
92
93 #define MG_MAXLINE 4096 /* maximum input line length */
94 #define MG_MAXARGC (MG_MAXLINE/4) /* maximum argument count */
95
96 typedef struct mg_fctxt {
97 char fname[96]; /* file name */
98 FILE *fp; /* stream pointer */
99 int fid; /* unique file context id */
100 char inpline[MG_MAXLINE]; /* input line */
101 int lineno; /* line number */
102 struct mg_fctxt *prev; /* previous context */
103 } MG_FCTXT;
104
105 typedef struct {
106 int fid; /* file this position is for */
107 int lineno; /* line number in file */
108 long offset; /* offset from beginning */
109 } MG_FPOS;
110
111 extern MG_FCTXT *mg_file; /* current file context */
112
113 #ifdef NOPROTO
114 extern void mg_init(); /* fill in mg_ehand array */
115 extern int mg_load(); /* parse a file */
116 extern int mg_open(); /* open new input file */
117 extern int mg_read(); /* read next line */
118 extern int mg_parse(); /* parse current line */
119 extern void mg_fgetpos(); /* get position on input file */
120 extern int mg_fgoto(); /* go to position on input file */
121 extern void mg_close(); /* close input file */
122 extern void mg_clear(); /* clear parser */
123 extern int mg_handle(); /* handle an entity */
124 #else
125 extern void mg_init(void); /* fill in mg_ehand array */
126 extern int mg_load(char *); /* parse a file */
127 extern int mg_open(MG_FCTXT *, char *); /* open new input file */
128 extern int mg_read(void); /* read next line */
129 extern int mg_parse(void); /* parse current line */
130 extern void mg_fgetpos(MG_FPOS *); /* get position on input file */
131 extern int mg_fgoto(MG_FPOS *); /* go to position on input file */
132 extern void mg_close(void); /* close input file */
133 extern void mg_clear(void); /* clear parser */
134 extern int mg_handle(int, int, char **); /* handle an entity */
135 #endif
136
137 #ifndef MG_NQCD
138 #define MG_NQCD 5 /* default number of divisions */
139 #endif
140
141 extern int mg_nqcdivs; /* divisions per quarter circle */
142
143 /*
144 * The following library routines are included for your convenience:
145 */
146
147 #ifdef NOPROTO
148 extern int mg_entity(); /* get entity number from its name */
149 extern int isint(); /* non-zero if integer format */
150 extern int isflt(); /* non-zero if floating point format */
151 #else
152 extern int mg_entity(char *); /* get entity number from its name */
153 extern int isint(char *); /* non-zero if integer format */
154 extern int isflt(char *); /* non-zero if floating point format */
155 #endif
156
157 /************************************************************************
158 * Definitions for 3-d vector manipulation functions
159 */
160
161 #ifdef SMLFLT
162 #define FLOAT float
163 #define FTINY (1e-3)
164 #else
165 #define FLOAT double
166 #define FTINY (1e-6)
167 #endif
168 #define FHUGE (1e10)
169
170 typedef FLOAT FVECT[3];
171
172 #define VCOPY(v1,v2) ((v1)[0]=(v2)[0],(v1)[1]=(v2)[1],(v1)[2]=(v2)[2])
173 #define DOT(v1,v2) ((v1)[0]*(v2)[0]+(v1)[1]*(v2)[1]+(v1)[2]*(v2)[2])
174 #define VSUM(vr,v1,v2,f) ((vr)[0]=(v1)[0]+(f)*(v2)[0], \
175 (vr)[1]=(v1)[1]+(f)*(v2)[1], \
176 (vr)[2]=(v1)[2]+(f)*(v2)[2])
177
178 #define is0vect(v) (DOT(v,v) <= FTINY*FTINY)
179
180 #define round0(x) if (x <= FTINY && x >= -FTINY) x = 0
181
182 #ifdef NOPROTO
183 extern double normalize(); /* normalize a vector */
184 #else
185 extern double normalize(FVECT); /* normalize a vector */
186 #endif
187
188 /************************************************************************
189 * Definitions for context handling routines
190 * (materials, colors, vectors)
191 */
192
193 #define C_CMINWL 380 /* minimum wavelength */
194 #define C_CMAXWL 780 /* maximum wavelength */
195 #define C_CNSS 41 /* number of spectral samples */
196 #define C_CWLI ((C_CMAXWL-C_CMINWL)/(C_CNSS-1))
197 #define C_CMAXV 10000 /* nominal maximum sample value */
198 #define C_CLPWM (683./C_CMAXV) /* peak lumens/watt multiplier */
199
200 #define C_CSSPEC 01 /* flag if spectrum is set */
201 #define C_CDSPEC 02 /* flag if defined w/ spectrum */
202 #define C_CSXY 04 /* flag if xy is set */
203 #define C_CDXY 010 /* flag if defined w/ xy */
204 #define C_CSEFF 020 /* flag if efficacy set */
205
206 typedef struct {
207 int clock; /* incremented each change */
208 unsigned short flags; /* what's been set */
209 unsigned short ssamp[C_CNSS]; /* spectral samples, min wl to max */
210 unsigned long ssum; /* straight sum of spectral values */
211 float cx, cy; /* xy chromaticity value */
212 float eff; /* efficacy (lumens/watt) */
213 } C_COLOR;
214
215 #define C_DEFCOLOR { 1, C_CDXY|C_CSXY|C_CSSPEC|C_CSEFF,\
216 {C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
217 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
218 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
219 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
220 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
221 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,\
222 C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV,C_CMAXV},\
223 (long)C_CNSS*C_CMAXV, 1./3., 1./3., 178.006 }
224
225 #define C_CIEX { 1, C_CDSPEC|C_CSSPEC|C_CSXY|C_CSEFF,\
226 {14,42,143,435,1344,2839,3483,3362,2908,1954,956,\
227 320,49,93,633,1655,2904,4334,5945,7621,9163,10263,\
228 10622,10026,8544,6424,4479,2835,1649,874,468,227,\
229 114,58,29,14,7,3,2,1,0}, 106836L, .467, .368, 362.230 }
230
231 #define C_CIEY { 1, C_CDSPEC|C_CSSPEC|C_CSXY|C_CSEFF,\
232 {0,1,4,12,40,116,230,380,600,910,1390,2080,3230,\
233 5030,7100,8620,9540,9950,9950,9520,8700,7570,6310,\
234 5030,3810,2650,1750,1070,610,320,170,82,41,21,10,\
235 5,2,1,1,0,0}, 106856L, .398, .542, 493.525 }
236
237 #define C_CIEZ { 1, C_CDSPEC|C_CSSPEC|C_CSXY|C_CSEFF,\
238 {65,201,679,2074,6456,13856,17471,17721,16692,\
239 12876,8130,4652,2720,1582,782,422,203,87,39,21,17,\
240 11,8,3,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},\
241 106770L, .147, .077, 54.363 }
242
243 #define c_cval(c,l) ((double)(c)->ssamp[((l)-C_MINWL)/C_CWLI] / (c)->ssum)
244
245 typedef struct {
246 int clock; /* incremented each change -- resettable */
247 int sided; /* 1 if surface is 1-sided, 0 for 2-sided */
248 float rd; /* diffuse reflectance */
249 C_COLOR rd_c; /* diffuse reflectance color */
250 float td; /* diffuse transmittance */
251 C_COLOR td_c; /* diffuse transmittance color */
252 float ed; /* diffuse emittance */
253 C_COLOR ed_c; /* diffuse emittance color */
254 float rs; /* specular reflectance */
255 C_COLOR rs_c; /* specular reflectance color */
256 float rs_a; /* specular reflectance roughness */
257 float ts; /* specular transmittance */
258 C_COLOR ts_c; /* specular transmittance color */
259 float ts_a; /* specular transmittance roughness */
260 } C_MATERIAL; /* material context */
261
262 typedef struct {
263 int clock; /* incremented each change -- resettable */
264 FVECT p, n; /* point and normal */
265 } C_VERTEX; /* vertex context */
266
267 #define C_DEFMATERIAL {1,0,0.,C_DEFCOLOR,0.,C_DEFCOLOR,0.,C_DEFCOLOR,\
268 0.,C_DEFCOLOR,0.,0.,C_DEFCOLOR,0.}
269 #define C_DEFVERTEX {1,{0.,0.,0.},{0.,0.,0.}}
270
271 extern C_COLOR *c_ccolor; /* the current color */
272 extern char *c_ccname; /* current color name */
273 extern C_MATERIAL *c_cmaterial; /* the current material */
274 extern char *c_cmname; /* current material name */
275 extern C_VERTEX *c_cvertex; /* the current vertex */
276 extern char *c_cvname; /* current vertex name */
277
278 #ifdef NOPROTO
279 extern int c_hcolor(); /* handle color entity */
280 extern int c_hmaterial(); /* handle material entity */
281 extern int c_hvertex(); /* handle vertex entity */
282 extern void c_clearall(); /* clear context tables */
283 extern C_MATERIAL *c_getmaterial(); /* get a named material */
284 extern C_VERTEX *c_getvert(); /* get a named vertex */
285 extern C_COLOR *c_getcolor(); /* get a named color */
286 extern void c_ccvt(); /* fix color representation */
287 extern int c_isgrey(); /* check if color is grey */
288 #else
289 extern int c_hcolor(int, char **); /* handle color entity */
290 extern int c_hmaterial(int, char **); /* handle material entity */
291 extern int c_hvertex(int, char **); /* handle vertex entity */
292 extern void c_clearall(void); /* clear context tables */
293 extern C_MATERIAL *c_getmaterial(char *); /* get a named material */
294 extern C_VERTEX *c_getvert(char *); /* get a named vertex */
295 extern C_COLOR *c_getcolor(char *); /* get a named color */
296 extern void c_ccvt(C_COLOR *, int); /* fix color representation */
297 extern int c_isgrey(C_COLOR *); /* check if color is grey */
298 #endif
299
300 /*************************************************************************
301 * Definitions for hierarchical object name handler
302 */
303
304 extern int obj_nnames; /* depth of name hierarchy */
305 extern char **obj_name; /* names in hierarchy */
306
307 #ifdef NOPROTO
308 extern int obj_handler(); /* handle an object entity */
309 extern void obj_clear(); /* clear object stack */
310 #else
311 extern int obj_handler(int, char **); /* handle an object entity */
312 extern void obj_clear(void); /* clear object stack */
313 #endif
314
315 /**************************************************************************
316 * Definitions for hierarchical transformation handler
317 */
318
319 typedef FLOAT MAT4[4][4];
320
321 #ifdef BSD
322 #define copymat4(m4a,m4b) bcopy((char *)m4b,(char *)m4a,sizeof(MAT4))
323 #else
324 #define copymat4(m4a,m4b) (void)memcpy((char *)m4a,(char *)m4b,sizeof(MAT4))
325 #endif
326
327 #define MAT4IDENT { {1.,0.,0.,0.}, {0.,1.,0.,0.}, \
328 {0.,0.,1.,0.}, {0.,0.,0.,1.} }
329
330 extern MAT4 m4ident;
331
332 #define setident4(m4) copymat4(m4, m4ident)
333
334 /* regular transformation */
335 typedef struct {
336 MAT4 xfm; /* transform matrix */
337 FLOAT sca; /* scalefactor */
338 } XF;
339
340 #define identxf(xp) (void)(setident4((xp)->xfm),(xp)->sca=1.0)
341
342 #define XF_MAXDIM 8 /* maximum array dimensions */
343
344 struct xf_array {
345 MG_FPOS spos; /* starting position on input */
346 int ndim; /* number of array dimensions */
347 struct {
348 short i, n; /* current count and maximum */
349 char arg[8]; /* string argument value */
350 } aarg[XF_MAXDIM];
351 };
352
353 typedef struct xf_spec {
354 long xid; /* unique transform id */
355 short xav0; /* zeroeth argument in xf_argv array */
356 short xac; /* transform argument count */
357 short rev; /* boolean true if vertices reversed */
358 XF xf; /* cumulative transformation */
359 struct xf_array *xarr; /* transformation array pointer */
360 struct xf_spec *prev; /* previous transformation context */
361 } XF_SPEC; /* followed by argument buffer */
362
363 extern int xf_argc; /* total # transform args. */
364 extern char **xf_argv; /* transform arguments */
365 extern XF_SPEC *xf_context; /* current context */
366
367 /*
368 * The transformation handler should do most of the work that needs
369 * doing. Just pass it any xf entities, then use the associated
370 * functions to transform and translate points, transform vectors
371 * (without translation), rotate vectors (without scaling) and scale
372 * values appropriately.
373 *
374 * The routines xf_xfmpoint, xf_xfmvect and xf_rotvect take two
375 * 3-D vectors (which may be identical), transforms the second and
376 * puts the result into the first.
377 */
378
379 #ifdef NOPROTO
380
381 extern int xf_handler(); /* handle xf entity */
382 extern void xf_xfmpoint(); /* transform point */
383 extern void xf_xfmvect(); /* transform vector */
384 extern void xf_rotvect(); /* rotate vector */
385 extern double xf_scale(); /* scale a value */
386 extern void xf_clear(); /* clear xf stack */
387
388 /* The following are support routines you probably won't call directly */
389
390 extern void multmat4(); /* m4a = m4b X m4c */
391 extern void multv3(); /* v3a = v3b X m4 (vectors) */
392 extern void multp3(); /* p3a = p3b X m4 (points) */
393 extern int xf(); /* interpret transform spec. */
394
395 #else
396
397 extern int xf_handler(int, char **); /* handle xf entity */
398 extern void xf_xfmpoint(FVECT, FVECT); /* transform point */
399 extern void xf_xfmvect(FVECT, FVECT); /* transform vector */
400 extern void xf_rotvect(FVECT, FVECT); /* rotate vector */
401 extern double xf_scale(double); /* scale a value */
402 extern void xf_clear(void); /* clear xf stack */
403
404 /* The following are support routines you probably won't call directly */
405
406 extern void multmat4(MAT4, MAT4, MAT4); /* m4a = m4b X m4c */
407 extern void multv3(FVECT, FVECT, MAT4); /* v3a = v3b X m4 (vectors) */
408 extern void multp3(FVECT, FVECT, MAT4); /* p3a = p3b X m4 (points) */
409 extern int xf(XF *, int, char **); /* interpret transform spec. */
410
411 #endif
412
413 /************************************************************************
414 * Miscellaneous definitions
415 */
416
417 #ifdef M_PI
418 #define PI M_PI
419 #else
420 #define PI 3.14159265358979323846
421 #endif
422
423 #ifdef DCL_ATOF
424 extern double atof();
425 #endif
426
427 #ifndef MEM_PTR
428 #define MEM_PTR void *
429 #endif
430
431 extern MEM_PTR malloc();
432 extern MEM_PTR calloc();
433 extern MEM_PTR realloc();
434 extern void free();