ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.30
Committed: Fri Mar 4 00:21:21 2016 UTC (8 years, 1 month ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad5R2, rad5R1
Changes since 2.29: +1 -6 lines
Log Message:
Eliminated redundant #define's for getc and putc

File Contents

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