ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.33
Committed: Thu Jun 8 17:39:13 2023 UTC (10 months, 3 weeks ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad5R4, HEAD
Changes since 2.32: +10 -2 lines
Log Message:
perf: Avoid calling nonsurfintree() when octree newer than input files

File Contents

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