ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.37
Committed: Sat Jun 7 05:09:45 2025 UTC (33 hours, 37 minutes ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: HEAD
Changes since 2.36: +1 -2 lines
Log Message:
refactor: Put some declarations into "paths.h" and included in "platform.h"

File Contents

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