ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.28
Committed: Thu Dec 16 19:31:36 2004 UTC (19 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad4R2P2, rad5R0, rad3R7P2, rad3R7P1, rad4R2, rad4R1, rad4R0, rad3R8, rad3R9, rad4R2P1
Changes since 2.27: +2 -2 lines
Log Message:
Fixed problem with last bug-fix (serious)

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: readoct.c,v 2.27 2004/12/15 22:02:04 greg Exp $";
3 #endif
4 /*
5 * readoct.c - routines to read octree information.
6 */
7
8 #include "copyright.h"
9
10 #include <stdio.h>
11 #include <time.h>
12
13 #include "platform.h"
14 #include "rtprocess.h"
15 #include "standard.h"
16 #include "octree.h"
17 #include "object.h"
18 #include "otypes.h"
19 #include "resolu.h"
20
21 #ifdef getc_unlocked /* avoid horrendous overhead of flockfile */
22 #undef getc
23 #define getc getc_unlocked
24 #endif
25
26 static double ogetflt(void);
27 static long ogetint(int);
28 static char *ogetstr(char *);
29 static int nonsurfintree(OCTREE ot);
30 static void octerror(int etyp, char *msg);
31 static void skiptree(void);
32 static OCTREE getfullnode(void), gettree(void);
33
34 static char *infn; /* input file specification */
35 static FILE *infp; /* input file stream */
36 static int objsize; /* size of stored OBJECT's */
37 static OBJECT objorig; /* zeroeth object */
38 static OBJECT fnobjects; /* number of objects in this file */
39
40
41 int
42 readoct( /* read in octree file or stream */
43 char *inpspec,
44 int load,
45 CUBE *scene,
46 char *ofn[]
47 )
48 {
49 char sbuf[512];
50 int nf;
51 int i;
52 long m;
53
54 if (inpspec == NULL) {
55 infn = "standard input";
56 infp = stdin;
57 } else if (inpspec[0] == '!') {
58 infn = inpspec;
59 if ((infp = popen(inpspec+1, "r")) == NULL) {
60 sprintf(errmsg, "cannot execute \"%s\"", inpspec);
61 error(SYSTEM, errmsg);
62 }
63 } else {
64 infn = inpspec;
65 if ((infp = fopen(inpspec, "r")) == NULL) {
66 sprintf(errmsg, "cannot open octree file \"%s\"",
67 inpspec);
68 error(SYSTEM, errmsg);
69 }
70 }
71 SET_FILE_BINARY(infp);
72 /* get header */
73 if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
74 octerror(USER, "not an octree");
75 /* check format */
76 if ((objsize = ogetint(2)-OCTMAGIC) <= 0 ||
77 objsize > MAXOBJSIZ || objsize > sizeof(long))
78 octerror(USER, "incompatible octree format");
79 /* get boundaries */
80 if (load & IO_BOUNDS) {
81 for (i = 0; i < 3; i++)
82 scene->cuorg[i] = atof(ogetstr(sbuf));
83 scene->cusize = atof(ogetstr(sbuf));
84 } else {
85 for (i = 0; i < 4; i++)
86 ogetstr(sbuf);
87 }
88 objorig = nobjects; /* set object offset */
89 nf = 0; /* get object files */
90 while (*ogetstr(sbuf)) {
91 if (load & IO_SCENE)
92 readobj(sbuf);
93 if (load & IO_FILES)
94 ofn[nf] = savqstr(sbuf);
95 nf++;
96 }
97 if (load & IO_FILES)
98 ofn[nf] = NULL;
99 /* get number of objects */
100 fnobjects = m = ogetint(objsize);
101 if (fnobjects != m)
102 octerror(USER, "too many objects");
103
104 if (load & IO_TREE) /* get the octree */
105 scene->cutree = gettree();
106 else if (load & IO_SCENE && nf == 0)
107 skiptree();
108
109 if (load & IO_SCENE) { /* get the scene */
110 if (nf == 0) {
111 /* load binary scene data */
112 readscene(infp, objsize);
113
114 } else { /* consistency checks */
115 /* check object count */
116 if (nobjects != objorig+fnobjects)
117 octerror(USER, "bad object count; octree stale?");
118 /* check for non-surfaces */
119 if (nonsurfintree(scene->cutree))
120 octerror(USER, "modifier in tree; octree stale?");
121 }
122 }
123 /* close the input */
124 if (infn[0] == '!')
125 pclose(infp);
126 else
127 fclose(infp);
128 return(nf);
129 }
130
131
132 static char *
133 ogetstr(char *s) /* get null-terminated string */
134 {
135 extern char *getstr();
136
137 if (getstr(s, infp) == NULL)
138 octerror(USER, "truncated octree");
139 return(s);
140 }
141
142
143 static OCTREE
144 getfullnode() /* get a set, return fullnode */
145 {
146 OBJECT set[MAXSET+1];
147 register int i;
148 register long m;
149
150 if ((set[0] = ogetint(objsize)) > MAXSET)
151 octerror(USER, "bad set in getfullnode");
152 for (i = 1; i <= set[0]; i++) {
153 m = ogetint(objsize) + objorig;
154 if ((set[i] = m) != m)
155 octerror(USER, "too many objects");
156 }
157 return(fullnode(set));
158 }
159
160
161 static long
162 ogetint(int siz) /* get a siz-byte integer */
163 {
164 extern long getint();
165 register long r;
166
167 r = getint(siz, infp);
168 if (feof(infp))
169 octerror(USER, "truncated octree");
170 return(r);
171 }
172
173
174 static double
175 ogetflt() /* get a floating point number */
176 {
177 extern double getflt();
178 double r;
179
180 r = getflt(infp);
181 if (feof(infp))
182 octerror(USER, "truncated octree");
183 return(r);
184 }
185
186
187 static OCTREE
188 gettree() /* get a pre-ordered octree */
189 {
190 register OCTREE ot;
191 register int i;
192
193 switch (getc(infp)) {
194 case OT_EMPTY:
195 return(EMPTY);
196 case OT_FULL:
197 return(getfullnode());
198 case OT_TREE:
199 if ((ot = octalloc()) == EMPTY)
200 octerror(SYSTEM, "out of tree space in gettree");
201 for (i = 0; i < 8; i++)
202 octkid(ot, i) = gettree();
203 return(ot);
204 case EOF:
205 octerror(USER, "truncated octree");
206 default:
207 octerror(USER, "damaged octree");
208 }
209 return EMPTY; /* pro forma return */
210 }
211
212
213 static int
214 nonsurfintree(OCTREE ot) /* check tree for modifiers */
215 {
216 OBJECT set[MAXSET+1];
217 register int i;
218
219 if (isempty(ot))
220 return(0);
221 if (istree(ot)) {
222 for (i = 0; i < 8; i++)
223 if (nonsurfintree(octkid(ot, i)))
224 return(1);
225 return(0);
226 }
227 objset(set, ot);
228 for (i = set[0]; i > 0; i-- )
229 if (ismodifier(objptr(set[i])->otype))
230 return(1);
231 return(0);
232 }
233
234
235 static void
236 skiptree(void) /* skip octree on input */
237 {
238 register int i;
239
240 switch (getc(infp)) {
241 case OT_EMPTY:
242 return;
243 case OT_FULL:
244 for (i = ogetint(objsize)*objsize; i-- > 0; )
245 if (getc(infp) == EOF)
246 octerror(USER, "truncated octree");
247 return;
248 case OT_TREE:
249 for (i = 0; i < 8; i++)
250 skiptree();
251 return;
252 case EOF:
253 octerror(USER, "truncated octree");
254 default:
255 octerror(USER, "damaged octree");
256 }
257 }
258
259
260 static void
261 octerror(int etyp, char *msg) /* octree error */
262 {
263 char msgbuf[128];
264
265 sprintf(msgbuf, "(%s): %s", infn, msg);
266 error(etyp, msgbuf);
267 }