ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.14
Committed: Tue Feb 25 02:47:21 2003 UTC (21 years, 2 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.13: +1 -56 lines
Log Message:
Replaced inline copyright notice with #include "copyright.h"

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2 greg 2.13 static const char RCSid[] = "$Id$";
3 greg 1.1 #endif
4     /*
5     * readoct.c - routines to read octree information.
6 greg 2.13 */
7    
8 greg 2.14 #include "copyright.h"
9 greg 1.1
10     #include "standard.h"
11    
12     #include "octree.h"
13    
14     #include "object.h"
15    
16     #include "otypes.h"
17    
18 greg 2.4 static double ogetflt();
19     static long ogetint();
20     static char *ogetstr();
21 greg 2.12 static int nonsurfinset();
22 greg 2.10 static int getobj(), octerror(), skiptree();
23 greg 1.7 static OCTREE getfullnode(), gettree();
24 greg 1.1
25 greg 2.13 static char *infn; /* input file specification */
26 greg 1.1 static FILE *infp; /* input file stream */
27 greg 1.14 static int objsize; /* size of stored OBJECT's */
28 greg 1.1 static OBJECT objorig; /* zeroeth object */
29 greg 2.12 static OBJECT fnobjects; /* number of objects in this file */
30 greg 1.1 static short otypmap[NUMOTYPE+8]; /* object type map */
31    
32    
33     int
34 greg 2.13 readoct(inpspec, load, scene, ofn) /* read in octree file or stream */
35     char *inpspec;
36 greg 1.1 int load;
37     CUBE *scene;
38     char *ofn[];
39     {
40 greg 1.11 char sbuf[512];
41 greg 1.1 int nf;
42     register int i;
43 greg 1.15 long m;
44 greg 1.1
45 greg 2.13 if (inpspec == NULL) {
46 greg 1.1 infn = "standard input";
47     infp = stdin;
48 greg 2.13 } else if (inpspec[0] == '!') {
49     infn = inpspec;
50     if ((infp = popen(inpspec+1, "r")) == NULL) {
51     sprintf(errmsg, "cannot execute \"%s\"", inpspec);
52     error(SYSTEM, errmsg);
53     }
54 greg 1.1 } else {
55 greg 2.13 infn = inpspec;
56     if ((infp = fopen(inpspec, "r")) == NULL) {
57 greg 1.1 sprintf(errmsg, "cannot open octree file \"%s\"",
58 greg 2.13 inpspec);
59 greg 1.1 error(SYSTEM, errmsg);
60     }
61     }
62 greg 2.7 #ifdef MSDOS
63     setmode(fileno(infp), O_BINARY);
64     #endif
65 greg 1.1 /* get header */
66 greg 2.11 if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
67 greg 1.12 octerror(USER, "not an octree");
68 greg 1.1 /* check format */
69 greg 2.4 if ((objsize = ogetint(2)-OCTMAGIC) <= 0 ||
70 greg 1.14 objsize > MAXOBJSIZ || objsize > sizeof(long))
71 greg 1.12 octerror(USER, "incompatible octree format");
72 greg 1.1 /* get boundaries */
73     if (load & IO_BOUNDS) {
74     for (i = 0; i < 3; i++)
75 greg 2.4 scene->cuorg[i] = atof(ogetstr(sbuf));
76     scene->cusize = atof(ogetstr(sbuf));
77 greg 1.1 } else {
78     for (i = 0; i < 4; i++)
79 greg 2.4 ogetstr(sbuf);
80 greg 1.1 }
81     objorig = nobjects; /* set object offset */
82     nf = 0; /* get object files */
83 greg 2.4 while (*ogetstr(sbuf)) {
84 greg 1.1 if (load & IO_SCENE)
85     readobj(sbuf);
86     if (load & IO_FILES)
87     ofn[nf] = savqstr(sbuf);
88     nf++;
89     }
90     if (load & IO_FILES)
91     ofn[nf] = NULL;
92 greg 1.8 /* get number of objects */
93 greg 2.4 fnobjects = m = ogetint(objsize);
94 greg 1.15 if (fnobjects != m)
95     octerror(USER, "too many objects");
96 greg 1.1
97 greg 2.13 if (load & IO_TREE) /* get the octree */
98 greg 1.1 scene->cutree = gettree();
99 greg 2.5 else if (load & IO_SCENE && nf == 0)
100     skiptree();
101    
102     if (load & IO_SCENE) /* get the scene */
103     if (nf == 0) {
104     for (i = 0; *ogetstr(sbuf); i++)
105     if ((otypmap[i] = otype(sbuf)) < 0) {
106     sprintf(errmsg, "unknown type \"%s\"", sbuf);
107     octerror(WARNING, errmsg);
108     }
109     while (getobj() != OVOID)
110     ;
111     } else { /* consistency checks */
112     /* check object count */
113     if (nobjects != objorig+fnobjects)
114     octerror(USER, "bad object count; octree stale?");
115     /* check for non-surfaces */
116 greg 2.12 if (dosets(nonsurfinset))
117 greg 2.5 octerror(USER, "modifier in tree; octree stale?");
118     }
119 greg 2.13 /* close the input */
120     if (infn[0] == '!')
121     pclose(infp);
122     else
123     fclose(infp);
124 greg 1.1 return(nf);
125     }
126    
127    
128     static char *
129 greg 2.4 ogetstr(s) /* get null-terminated string */
130 greg 1.1 char *s;
131     {
132 greg 2.4 extern char *getstr();
133 greg 1.1
134 greg 2.4 if (getstr(s, infp) == NULL)
135     octerror(USER, "truncated octree");
136     return(s);
137 greg 1.1 }
138    
139    
140     static OCTREE
141     getfullnode() /* get a set, return fullnode */
142     {
143 greg 2.7 OBJECT set[MAXSET+1];
144 greg 1.1 register int i;
145 greg 1.14 register long m;
146 greg 1.1
147 greg 2.4 if ((set[0] = ogetint(objsize)) > MAXSET)
148 greg 1.1 octerror(USER, "bad set in getfullnode");
149 greg 1.14 for (i = 1; i <= set[0]; i++) {
150 greg 2.4 m = ogetint(objsize) + objorig;
151 greg 1.15 if ((set[i] = m) != m)
152 greg 1.14 octerror(USER, "too many objects");
153     }
154 greg 1.1 return(fullnode(set));
155     }
156    
157    
158     static long
159 greg 2.4 ogetint(siz) /* get a siz-byte integer */
160     int siz;
161 greg 1.1 {
162 greg 2.4 extern long getint();
163 greg 1.3 register long r;
164 greg 1.1
165 greg 2.4 r = getint(siz, infp);
166     if (feof(infp))
167     octerror(USER, "truncated octree");
168 greg 1.1 return(r);
169     }
170    
171    
172     static double
173 greg 2.4 ogetflt() /* get a floating point number */
174 greg 1.1 {
175 greg 2.4 extern double getflt();
176 greg 2.7 double r;
177 greg 1.1
178 greg 2.4 r = getflt(infp);
179     if (feof(infp))
180     octerror(USER, "truncated octree");
181     return(r);
182 greg 1.1 }
183    
184    
185     static OCTREE
186     gettree() /* get a pre-ordered octree */
187     {
188 greg 2.7 register OCTREE ot;
189 greg 1.1 register int i;
190    
191     switch (getc(infp)) {
192     case OT_EMPTY:
193     return(EMPTY);
194     case OT_FULL:
195     return(getfullnode());
196     case OT_TREE:
197     if ((ot = octalloc()) == EMPTY)
198     octerror(SYSTEM, "out of tree space in gettree");
199     for (i = 0; i < 8; i++)
200     octkid(ot, i) = gettree();
201     return(ot);
202 greg 2.5 case EOF:
203     octerror(USER, "truncated octree");
204     default:
205     octerror(USER, "damaged octree");
206     }
207 greg 2.12 }
208    
209    
210     static int
211     nonsurfinset(os) /* check set for modifier */
212     register OBJECT *os;
213     {
214     register OBJECT s;
215     register int i;
216    
217     for (i = *os; i-- > 0; )
218     if ((s = *++os) >= objorig && s < objorig+fnobjects &&
219     ismodifier(objptr(s)->otype))
220     return(1);
221     return(0);
222 greg 2.5 }
223    
224    
225     static
226     skiptree() /* skip octree on input */
227     {
228     register int i;
229    
230     switch (getc(infp)) {
231     case OT_EMPTY:
232     return;
233     case OT_FULL:
234 greg 2.6 for (i = ogetint(objsize)*objsize; i-- > 0; )
235     if (getc(infp) == EOF)
236     octerror(USER, "truncated octree");
237 greg 2.5 return;
238     case OT_TREE:
239     for (i = 0; i < 8; i++)
240     skiptree();
241     return;
242 greg 1.1 case EOF:
243     octerror(USER, "truncated octree");
244     default:
245     octerror(USER, "damaged octree");
246     }
247     }
248    
249    
250     static
251     getobj() /* get next object */
252     {
253     char sbuf[MAXSTR];
254     int obj;
255     register int i;
256 greg 1.14 register long m;
257 greg 2.7 register OBJREC *objp;
258 greg 1.1
259 greg 2.4 i = ogetint(1);
260 greg 1.3 if (i == -1)
261 greg 1.1 return(OVOID); /* terminator */
262     if ((obj = newobject()) == OVOID)
263     error(SYSTEM, "out of object space");
264     objp = objptr(obj);
265     if ((objp->otype = otypmap[i]) < 0)
266     octerror(USER, "reference to unknown type");
267 greg 2.4 if ((m = ogetint(objsize)) != OVOID) {
268 greg 1.14 m += objorig;
269     if ((OBJECT)m != m)
270     octerror(USER, "too many objects");
271     }
272     objp->omod = m;
273 greg 2.4 objp->oname = savqstr(ogetstr(sbuf));
274     if (objp->oargs.nsargs = ogetint(2)) {
275 greg 2.13 objp->oargs.sarg = (char **)malloc
276 greg 1.1 (objp->oargs.nsargs*sizeof(char *));
277     if (objp->oargs.sarg == NULL)
278     goto memerr;
279     for (i = 0; i < objp->oargs.nsargs; i++)
280 greg 2.4 objp->oargs.sarg[i] = savestr(ogetstr(sbuf));
281 greg 1.1 } else
282     objp->oargs.sarg = NULL;
283 greg 2.7 #ifdef IARGS
284 greg 2.4 if (objp->oargs.niargs = ogetint(2)) {
285 greg 2.13 objp->oargs.iarg = (long *)malloc
286 greg 1.1 (objp->oargs.niargs*sizeof(long));
287     if (objp->oargs.iarg == NULL)
288     goto memerr;
289     for (i = 0; i < objp->oargs.niargs; i++)
290 greg 2.4 objp->oargs.iarg[i] = ogetint(4);
291 greg 1.1 } else
292     objp->oargs.iarg = NULL;
293     #endif
294 greg 2.4 if (objp->oargs.nfargs = ogetint(2)) {
295 greg 2.13 objp->oargs.farg = (FLOAT *)malloc
296 greg 1.15 (objp->oargs.nfargs*sizeof(FLOAT));
297 greg 1.1 if (objp->oargs.farg == NULL)
298     goto memerr;
299     for (i = 0; i < objp->oargs.nfargs; i++)
300 greg 2.4 objp->oargs.farg[i] = ogetflt();
301 greg 1.1 } else
302     objp->oargs.farg = NULL;
303     /* initialize */
304     objp->os = NULL;
305     /* insert */
306     insertobject(obj);
307     return(obj);
308     memerr:
309     error(SYSTEM, "out of memory in getobj");
310     }
311    
312    
313     static
314     octerror(etyp, msg) /* octree error */
315     int etyp;
316     char *msg;
317     {
318     char msgbuf[128];
319    
320     sprintf(msgbuf, "(%s): %s", infn, msg);
321     error(etyp, msgbuf);
322     }