ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/readoct.c
(Generate patch)

Comparing ray/src/common/readoct.c (file contents):
Revision 1.11 by greg, Thu Feb 7 11:00:42 1991 UTC vs.
Revision 2.13 by greg, Sat Feb 22 02:07:22 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1986 Regents of the University of California */
2
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6
4   /*
5   *  readoct.c - routines to read octree information.
6 + */
7 +
8 + /* ====================================================================
9 + * The Radiance Software License, Version 1.0
10   *
11 < *     7/30/85
11 > * Copyright (c) 1990 - 2002 The Regents of the University of California,
12 > * through Lawrence Berkeley National Laboratory.   All rights reserved.
13 > *
14 > * Redistribution and use in source and binary forms, with or without
15 > * modification, are permitted provided that the following conditions
16 > * are met:
17 > *
18 > * 1. Redistributions of source code must retain the above copyright
19 > *         notice, this list of conditions and the following disclaimer.
20 > *
21 > * 2. Redistributions in binary form must reproduce the above copyright
22 > *       notice, this list of conditions and the following disclaimer in
23 > *       the documentation and/or other materials provided with the
24 > *       distribution.
25 > *
26 > * 3. The end-user documentation included with the redistribution,
27 > *           if any, must include the following acknowledgment:
28 > *             "This product includes Radiance software
29 > *                 (http://radsite.lbl.gov/)
30 > *                 developed by the Lawrence Berkeley National Laboratory
31 > *               (http://www.lbl.gov/)."
32 > *       Alternately, this acknowledgment may appear in the software itself,
33 > *       if and wherever such third-party acknowledgments normally appear.
34 > *
35 > * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
36 > *       and "The Regents of the University of California" must
37 > *       not be used to endorse or promote products derived from this
38 > *       software without prior written permission. For written
39 > *       permission, please contact [email protected].
40 > *
41 > * 5. Products derived from this software may not be called "Radiance",
42 > *       nor may "Radiance" appear in their name, without prior written
43 > *       permission of Lawrence Berkeley National Laboratory.
44 > *
45 > * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
46 > * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
47 > * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
48 > * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
49 > * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
50 > * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
51 > * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
52 > * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
53 > * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
54 > * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
55 > * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
56 > * SUCH DAMAGE.
57 > * ====================================================================
58 > *
59 > * This software consists of voluntary contributions made by many
60 > * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
61 > * information on Lawrence Berkeley National Laboratory, please see
62 > * <http://www.lbl.gov/>.
63   */
64  
65   #include  "standard.h"
# Line 18 | Line 70 | static char SCCSid[] = "$SunId$ LBL";
70  
71   #include  "otypes.h"
72  
73 < extern double  atof();
74 < static double  getflt();
75 < static long  getint();
76 < static char  *getstr();
73 > static double  ogetflt();
74 > static long  ogetint();
75 > static char  *ogetstr();
76 > static int  nonsurfinset();
77 > static int  getobj(), octerror(), skiptree();
78   static OCTREE  getfullnode(), gettree();
79  
80 < static char  *infn;                     /* input file name */
80 > static char  *infn;                     /* input file specification */
81   static FILE  *infp;                     /* input file stream */
82 + static int  objsize;                    /* size of stored OBJECT's */
83   static OBJECT  objorig;                 /* zeroeth object */
84 + static OBJECT  fnobjects;               /* number of objects in this file */
85   static short  otypmap[NUMOTYPE+8];      /* object type map */
86  
87  
88   int
89 < readoct(fname, load, scene, ofn)        /* read in octree from file */
90 < char  *fname;
89 > readoct(inpspec, load, scene, ofn)      /* read in octree file or stream */
90 > char  *inpspec;
91   int  load;
92   CUBE  *scene;
93   char  *ofn[];
94   {
95          char  sbuf[512];
96          int  nf;
42        OBJECT  fnobjects;
97          register int  i;
98 +        long  m;
99          
100 <        if (fname == NULL) {
100 >        if (inpspec == NULL) {
101                  infn = "standard input";
102                  infp = stdin;
103 +        } else if (inpspec[0] == '!') {
104 +                infn = inpspec;
105 +                if ((infp = popen(inpspec+1, "r")) == NULL) {
106 +                        sprintf(errmsg, "cannot execute \"%s\"", inpspec);
107 +                        error(SYSTEM, errmsg);
108 +                }
109          } else {
110 <                infn = fname;
111 <                if ((infp = fopen(fname, "r")) == NULL) {
110 >                infn = inpspec;
111 >                if ((infp = fopen(inpspec, "r")) == NULL) {
112                          sprintf(errmsg, "cannot open octree file \"%s\"",
113 <                                        fname);
113 >                                        inpspec);
114                          error(SYSTEM, errmsg);
115                  }
116          }
117 + #ifdef MSDOS
118 +        setmode(fileno(infp), O_BINARY);
119 + #endif
120                                          /* get header */
121 <        if (load & IO_INFO)
122 <                copyheader(infp, stdout);
59 <        else
60 <                getheader(infp, NULL);
121 >        if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
122 >                octerror(USER, "not an octree");
123                                          /* check format */
124 <        if (getint(2) != OCTMAGIC)
125 <                octerror(USER, "invalid octree format");
124 >        if ((objsize = ogetint(2)-OCTMAGIC) <= 0 ||
125 >                        objsize > MAXOBJSIZ || objsize > sizeof(long))
126 >                octerror(USER, "incompatible octree format");
127                                          /* get boundaries */
128          if (load & IO_BOUNDS) {
129                  for (i = 0; i < 3; i++)
130 <                        scene->cuorg[i] = atof(getstr(sbuf));
131 <                scene->cusize = atof(getstr(sbuf));
130 >                        scene->cuorg[i] = atof(ogetstr(sbuf));
131 >                scene->cusize = atof(ogetstr(sbuf));
132          } else {
133                  for (i = 0; i < 4; i++)
134 <                        getstr(sbuf);
134 >                        ogetstr(sbuf);
135          }
136          objorig = nobjects;             /* set object offset */
137          nf = 0;                         /* get object files */
138 <        while (*getstr(sbuf)) {
138 >        while (*ogetstr(sbuf)) {
139                  if (load & IO_SCENE)
140                          readobj(sbuf);
141                  if (load & IO_FILES)
# Line 82 | Line 145 | char  *ofn[];
145          if (load & IO_FILES)
146                  ofn[nf] = NULL;
147                                          /* get number of objects */
148 <        fnobjects = getint(sizeof(OBJECT));
148 >        fnobjects = m = ogetint(objsize);
149 >        if (fnobjects != m)
150 >                octerror(USER, "too many objects");
151  
152 <        if (load & IO_TREE) {
88 <                                                /* get the octree */
152 >        if (load & IO_TREE)             /* get the octree */
153                  scene->cutree = gettree();
154 <                                                /* get the scene */
155 <                if (nf == 0 && load & IO_SCENE) {
156 <                        for (i = 0; *getstr(sbuf); i++)
157 <                                if ((otypmap[i] = otype(sbuf)) < 0) {
158 <                                        sprintf(errmsg, "unknown type \"%s\"",
159 <                                                        sbuf);
160 <                                        octerror(WARNING, errmsg);
161 <                                }
162 <                        while (getobj() != OVOID)
163 <                                ;
164 <                }
165 <        }
166 <        fclose(infp);
167 <                                        /* consistency checks */
104 <        if (load & IO_SCENE) {
105 <                                        /* check object count */
154 >        else if (load & IO_SCENE && nf == 0)
155 >                skiptree();
156 >                
157 >        if (load & IO_SCENE)            /* get the scene */
158 >            if (nf == 0) {
159 >                for (i = 0; *ogetstr(sbuf); i++)
160 >                        if ((otypmap[i] = otype(sbuf)) < 0) {
161 >                                sprintf(errmsg, "unknown type \"%s\"", sbuf);
162 >                                octerror(WARNING, errmsg);
163 >                        }
164 >                while (getobj() != OVOID)
165 >                        ;
166 >            } else {                    /* consistency checks */
167 >                                /* check object count */
168                  if (nobjects != objorig+fnobjects)
169                          octerror(USER, "bad object count; octree stale?");
170 <                                        /* check for non-surfaces */
171 <                if (nonsurfinset(objorig, fnobjects))
172 <                        octerror(USER, "non-surface in set; octree stale?");
173 <        }
170 >                                /* check for non-surfaces */
171 >                if (dosets(nonsurfinset))
172 >                        octerror(USER, "modifier in tree; octree stale?");
173 >            }
174 >                                /* close the input */
175 >        if (infn[0] == '!')
176 >                pclose(infp);
177 >        else
178 >                fclose(infp);
179          return(nf);
180   }
181  
182  
183   static char *
184 < getstr(s)                       /* get null-terminated string */
184 > ogetstr(s)                      /* get null-terminated string */
185   char  *s;
186   {
187 <        register char  *cp;
121 <        register int  c;
187 >        extern char  *getstr();
188  
189 <        cp = s;
190 <        while ((c = getc(infp)) != EOF)
191 <                if ((*cp++ = c) == '\0')
126 <                        return(s);
127 <                
128 <        octerror(USER, "truncated octree");
189 >        if (getstr(s, infp) == NULL)
190 >                octerror(USER, "truncated octree");
191 >        return(s);
192   }
193  
194  
195   static OCTREE
196   getfullnode()                   /* get a set, return fullnode */
197   {
198 <        OBJECT  set[MAXSET+1];
198 >        OBJECT  set[MAXSET+1];
199          register int  i;
200 +        register long  m;
201  
202 <        set[0] = getint(sizeof(OBJECT));
139 <        if (set[0] > MAXSET)
202 >        if ((set[0] = ogetint(objsize)) > MAXSET)
203                  octerror(USER, "bad set in getfullnode");
204 <        for (i = 1; i <= set[0]; i++)
205 <                set[i] = getint(sizeof(OBJECT)) + objorig;
204 >        for (i = 1; i <= set[0]; i++) {
205 >                m = ogetint(objsize) + objorig;
206 >                if ((set[i] = m) != m)
207 >                        octerror(USER, "too many objects");
208 >        }
209          return(fullnode(set));
210   }      
211  
212  
213   static long
214 < getint(siz)                     /* get a siz-byte integer */
215 < register int  siz;
214 > ogetint(siz)                    /* get a siz-byte integer */
215 > int  siz;
216   {
217 <        register int  c;
217 >        extern long  getint();
218          register long  r;
219  
220 <        if ((c = getc(infp)) == EOF)
221 <                goto end_file;
222 <        r = 0x80&c ? -1<<8|c : c;               /* sign extend */
157 <        while (--siz > 0) {
158 <                if ((c = getc(infp)) == EOF)
159 <                        goto end_file;
160 <                r <<= 8;
161 <                r |= c;
162 <        }
220 >        r = getint(siz, infp);
221 >        if (feof(infp))
222 >                octerror(USER, "truncated octree");
223          return(r);
164 end_file:
165        octerror(USER, "truncated octree");
224   }
225  
226  
227   static double
228 < getflt()                        /* get a floating point number */
228 > ogetflt()                       /* get a floating point number */
229   {
230 <        extern double  ldexp();
231 <        double  d;
230 >        extern double  getflt();
231 >        double  r;
232  
233 <        d = (double)getint(4)/0x7fffffff;
234 <        return(ldexp(d, (int)getint(1)));
233 >        r = getflt(infp);
234 >        if (feof(infp))
235 >                octerror(USER, "truncated octree");
236 >        return(r);
237   }
238          
239  
240   static OCTREE
241   gettree()                       /* get a pre-ordered octree */
242   {
243 <        register OCTREE  ot;
243 >        register OCTREE  ot;
244          register int  i;
245          
246          switch (getc(infp)) {
# Line 202 | Line 262 | gettree()                      /* get a pre-ordered octree */
262   }
263  
264  
265 + static int
266 + nonsurfinset(os)                        /* check set for modifier */
267 + register OBJECT  *os;
268 + {
269 +        register OBJECT  s;
270 +        register int  i;
271 +
272 +        for (i = *os; i-- > 0; )
273 +                if ((s = *++os) >= objorig && s < objorig+fnobjects &&
274 +                                ismodifier(objptr(s)->otype))
275 +                        return(1);
276 +        return(0);
277 + }
278 +
279 +
280   static
281 + skiptree()                              /* skip octree on input */
282 + {
283 +        register int  i;
284 +        
285 +        switch (getc(infp)) {
286 +        case OT_EMPTY:
287 +                return;
288 +        case OT_FULL:
289 +                for (i = ogetint(objsize)*objsize; i-- > 0; )
290 +                        if (getc(infp) == EOF)
291 +                                octerror(USER, "truncated octree");
292 +                return;
293 +        case OT_TREE:
294 +                for (i = 0; i < 8; i++)
295 +                        skiptree();
296 +                return;
297 +        case EOF:
298 +                octerror(USER, "truncated octree");
299 +        default:
300 +                octerror(USER, "damaged octree");
301 +        }
302 + }
303 +
304 +
305 + static
306   getobj()                                /* get next object */
307   {
308          char  sbuf[MAXSTR];
309          int  obj;
310          register int  i;
311 <        register OBJREC  *objp;
311 >        register long  m;
312 >        register OBJREC  *objp;
313  
314 <        i = getint(1);
314 >        i = ogetint(1);
315          if (i == -1)
316                  return(OVOID);          /* terminator */
317          if ((obj = newobject()) == OVOID)
# Line 218 | Line 319 | getobj()                               /* get next object */
319          objp = objptr(obj);
320          if ((objp->otype = otypmap[i]) < 0)
321                  octerror(USER, "reference to unknown type");
322 <        if ((objp->omod = getint(sizeof(OBJECT))) != OVOID)
323 <                objp->omod += objorig;
324 <        objp->oname = savqstr(getstr(sbuf));
325 <        if (objp->oargs.nsargs = getint(2)) {
326 <                objp->oargs.sarg = (char **)bmalloc
322 >        if ((m = ogetint(objsize)) != OVOID) {
323 >                m += objorig;
324 >                if ((OBJECT)m != m)
325 >                        octerror(USER, "too many objects");
326 >        }
327 >        objp->omod = m;
328 >        objp->oname = savqstr(ogetstr(sbuf));
329 >        if (objp->oargs.nsargs = ogetint(2)) {
330 >                objp->oargs.sarg = (char **)malloc
331                                  (objp->oargs.nsargs*sizeof(char *));
332                  if (objp->oargs.sarg == NULL)
333                          goto memerr;
334                  for (i = 0; i < objp->oargs.nsargs; i++)
335 <                        objp->oargs.sarg[i] = savestr(getstr(sbuf));
335 >                        objp->oargs.sarg[i] = savestr(ogetstr(sbuf));
336          } else
337                  objp->oargs.sarg = NULL;
338 < #ifdef  IARGS
339 <        if (objp->oargs.niargs = getint(2)) {
340 <                objp->oargs.iarg = (long *)bmalloc
338 > #ifdef  IARGS
339 >        if (objp->oargs.niargs = ogetint(2)) {
340 >                objp->oargs.iarg = (long *)malloc
341                                  (objp->oargs.niargs*sizeof(long));
342                  if (objp->oargs.iarg == NULL)
343                          goto memerr;
344                  for (i = 0; i < objp->oargs.niargs; i++)
345 <                        objp->oargs.iarg[i] = getint(4);
345 >                        objp->oargs.iarg[i] = ogetint(4);
346          } else
347                  objp->oargs.iarg = NULL;
348   #endif
349 <        if (objp->oargs.nfargs = getint(2)) {
350 <                objp->oargs.farg = (double *)bmalloc
351 <                                (objp->oargs.nfargs*sizeof(double));
349 >        if (objp->oargs.nfargs = ogetint(2)) {
350 >                objp->oargs.farg = (FLOAT *)malloc
351 >                                (objp->oargs.nfargs*sizeof(FLOAT));
352                  if (objp->oargs.farg == NULL)
353                          goto memerr;
354                  for (i = 0; i < objp->oargs.nfargs; i++)
355 <                        objp->oargs.farg[i] = getflt();
355 >                        objp->oargs.farg[i] = ogetflt();
356          } else
357                  objp->oargs.farg = NULL;
358                                                  /* initialize */
359          objp->os = NULL;
255        objp->lastrno = -1;
360                                                  /* insert */
361          insertobject(obj);
362          return(obj);

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines