cv/mgflib/parser.h

/* Copyright (c) 1995 Regents of the University of California */

/* SCCSid "$SunId$ LBL" */

/*
 * Header file for MGF interpreter
 */

/* must include stdio.h before us */

                        /* Entities (order doesn't really matter) */
#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           */

#define MG_NENTITIES    28              /* total # entities */

#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 */

#ifdef NOPROTO
extern int      (*mg_ehand[MG_NENTITIES])();
#else
extern int      (*mg_ehand[MG_NENTITIES])(int argc, char **argv);
#endif

                        /* Error codes */
#define MG_OK           0               /* normal return value */
#define MG_EUNK         1               /* unknown entity */
#define MG_EARGC        2               /* wrong number of arguments */
#define MG_ETYPE        3               /* argument type error */
#define MG_EILL         4               /* illegal argument value */
#define MG_EUNDEF       5               /* undefined reference */
#define MG_ENOFILE      6               /* cannot open input file */
#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        11

extern char     *mg_err[MG_NERRS];

/*
 * The general process for running the parser is to fill in the mg_ehand
 * array with handlers for each entity you know how to handle.
 * Then, call mg_init to fill in the rest.  This function will report
 * an error and quit if you try to support an inconsistent set of entities.
 * For each file you want to parse, call mg_load with the file name.
 * 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, 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      4096            /* maximum input line length */
#define MG_MAXARGC      (MG_MAXLINE/4)  /* maximum argument count */

typedef struct mg_fctxt {
        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
extern void     mg_init();              /* fill in mg_ehand array */
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 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_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 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_handle(int, int, char **);   /* handle an entity */
#endif

#ifndef MG_NQCD
#define MG_NQCD         5               /* default number of divisions */
#endif

extern int      mg_nqcdivs;             /* divisions per quarter circle */

/*
 * The following library routines are included for your convenience:
 */

#ifdef NOPROTO
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 */
#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 */
#endif

/************************************************************************
 *      Definitions for 3-d vector manipulation functions
 */

#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
extern double   normalize(FVECT);       /* normalize a vector */
#endif

/************************************************************************
 *      Definitions for context handling routines
 *      (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 {
        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 {
        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 */
        C_COLOR td_c;           /* diffuse transmittance color */
        float   ed;             /* diffuse emittance */
        C_COLOR ed_c;           /* diffuse emittance color */
        float   rs;             /* specular reflectance */
        C_COLOR rs_c;           /* specular reflectance color */
        float   rs_a;           /* specular reflectance roughness */
        float   ts;             /* specular transmittance */
        C_COLOR ts_c;           /* specular transmittance color */
        float   ts_a;           /* specular transmittance roughness */
} C_MATERIAL;           /* material context */

typedef struct {
        int     clock;          /* incremented each change -- resettable */
        FVECT   p, n;           /* point and normal */
} C_VERTEX;             /* vertex context */

#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     {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

/*************************************************************************
 *      Definitions for hierarchical object name handler
 */

extern int      obj_nnames;             /* depth of name hierarchy */
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 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))
#endif

#define  MAT4IDENT              { {1.,0.,0.,0.}, {0.,1.,0.,0.}, \
                                {0.,0.,1.,0.}, {0.,0.,0.,1.} }

extern MAT4  m4ident;

#define  setident4(m4)          copymat4(m4, m4ident)

                                /* regular transformation */
typedef struct {
        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 {
        long    xid;                    /* unique transform id */
        short   xav0;                   /* zeroeth argument in xf_argv array */
        short   xac;                    /* transform 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 */

/*
 * The transformation handler should do most of the work that needs
 * doing.  Just pass it any xf entities, then use the associated
 * functions to transform and translate points, transform vectors
 * (without translation), rotate vectors (without scaling) and scale
 * values appropriately.
 *
 * The routines xf_xfmpoint, xf_xfmvect and xf_rotvect take two
 * 3-D vectors (which may be identical), transforms the second and
 * puts the result into the first.
 */

#ifdef NOPROTO

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 void     xf_clear();             /* clear xf stack */

/* The following are support routines you probably won't call directly */

extern void     multmat4();             /* m4a = m4b X m4c */
extern void     multv3();               /* v3a = v3b X m4 (vectors) */
extern void     multp3();               /* p3a = p3b X m4 (points) */
extern int      xf();                   /* interpret transform spec. */

#else

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. */

#endif

/************************************************************************
 *      Miscellaneous definitions
 */

#ifdef  M_PI
#define  PI             M_PI
#else
#define  PI             3.14159265358979323846
#endif

#ifdef DCL_ATOF
extern double   atof();
#endif

#ifndef MEM_PTR
#define MEM_PTR         void *
#endif

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