ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.24
Committed: Mon Oct 27 10:19:31 2003 UTC (20 years, 6 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.23: +2 -1 lines
Log Message:
Added gethomedir.c and various compatibility fixes.

File Contents

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