ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
Revision: 2.15
Committed: Fri Mar 14 21:27:46 2003 UTC (21 years, 1 month ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad3R5
Changes since 2.14: +5 -73 lines
Log Message:
Added -a option to obj2mesh to incorporate materials in mesh

File Contents

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