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.16 by greg, Mon Oct 28 09:47:22 1991 UTC vs.
Revision 2.27 by greg, Wed Dec 15 22:02:04 2004 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1991 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.
9 *
10 *     7/30/85
6   */
7  
8 < #include  "standard.h"
8 > #include "copyright.h"
9  
10 < #include  "octree.h"
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
18   #include  "otypes.h"
19 + #include  "resolu.h"
20  
21 < extern double  atof();
22 < static double  getflt();
23 < static long  getint();
24 < static char  *getstr();
25 < static int  getobj(), octerror();
26 < static OCTREE  getfullnode(), gettree();
21 > #ifdef getc_unlocked            /* avoid horrendous overhead of flockfile */
22 > #undef getc
23 > #define getc    getc_unlocked
24 > #endif
25  
26 < static char  *infn;                     /* input file name */
26 > static double  ogetflt(void);
27 > static long  ogetint(int);
28 > static char  *ogetstr(char *);
29 > static int  nonsurfintree(OCTREE ot);
30 > static void  octerror(int  etyp, char  *msg);
31 > static void  skiptree(void);
32 > static OCTREE  getfullnode(void), gettree(void);
33 >
34 > static char  *infn;                     /* input file specification */
35   static FILE  *infp;                     /* input file stream */
36   static int  objsize;                    /* size of stored OBJECT's */
37   static OBJECT  objorig;                 /* zeroeth object */
38 < static short  otypmap[NUMOTYPE+8];      /* object type map */
38 > static OBJECT  fnobjects;               /* number of objects in this file */
39  
40  
41   int
42 < readoct(fname, load, scene, ofn)        /* read in octree from file */
43 < char  *fname;
44 < int  load;
45 < CUBE  *scene;
46 < char  *ofn[];
42 > readoct(                                /* read in octree file or stream */
43 >        char  *inpspec,
44 >        int  load,
45 >        CUBE  *scene,
46 >        char  *ofn[]
47 > )
48   {
42        extern int  fputs();
49          char  sbuf[512];
50          int  nf;
51 <        OBJECT  fnobjects;
46 <        register int  i;
51 >        int  i;
52          long  m;
53          
54 <        if (fname == NULL) {
54 >        if (inpspec == NULL) {
55                  infn = "standard input";
56                  infp = stdin;
57 +        } else if (inpspec[0] == '!') {
58 +                infn = inpspec;
59 +                if ((infp = popen(inpspec+1, "r")) == NULL) {
60 +                        sprintf(errmsg, "cannot execute \"%s\"", inpspec);
61 +                        error(SYSTEM, errmsg);
62 +                }
63          } else {
64 <                infn = fname;
65 <                if ((infp = fopen(fname, "r")) == NULL) {
64 >                infn = inpspec;
65 >                if ((infp = fopen(inpspec, "r")) == NULL) {
66                          sprintf(errmsg, "cannot open octree file \"%s\"",
67 <                                        fname);
67 >                                        inpspec);
68                          error(SYSTEM, errmsg);
69                  }
70          }
71 +        SET_FILE_BINARY(infp);
72                                          /* get header */
73 <        if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : NULL) < 0)
73 >        if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
74                  octerror(USER, "not an octree");
75                                          /* check format */
76 <        if ((objsize = getint(2)-OCTMAGIC) <= 0 ||
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(getstr(sbuf));
83 <                scene->cusize = atof(getstr(sbuf));
82 >                        scene->cuorg[i] = atof(ogetstr(sbuf));
83 >                scene->cusize = atof(ogetstr(sbuf));
84          } else {
85                  for (i = 0; i < 4; i++)
86 <                        getstr(sbuf);
86 >                        ogetstr(sbuf);
87          }
88          objorig = nobjects;             /* set object offset */
89          nf = 0;                         /* get object files */
90 <        while (*getstr(sbuf)) {
90 >        while (*ogetstr(sbuf)) {
91                  if (load & IO_SCENE)
92                          readobj(sbuf);
93                  if (load & IO_FILES)
# Line 85 | Line 97 | char  *ofn[];
97          if (load & IO_FILES)
98                  ofn[nf] = NULL;
99                                          /* get number of objects */
100 <        fnobjects = m = getint(objsize);
100 >        fnobjects = m = ogetint(objsize);
101          if (fnobjects != m)
102                  octerror(USER, "too many objects");
103  
104 <        if (load & IO_TREE) {
93 <                                                /* get the octree */
104 >        if (load & IO_TREE)             /* get the octree */
105                  scene->cutree = gettree();
106 <                                                /* get the scene */
107 <                if (nf == 0 && load & IO_SCENE) {
108 <                        for (i = 0; *getstr(sbuf); i++)
109 <                                if ((otypmap[i] = otype(sbuf)) < 0) {
110 <                                        sprintf(errmsg, "unknown type \"%s\"",
111 <                                                        sbuf);
112 <                                        octerror(WARNING, errmsg);
113 <                                }
114 <                        while (getobj() != OVOID)
115 <                                ;
116 <                } else if (load & IO_SCENE) {           /* consistency checks */
117 <                                                /* check object count */
118 <                        if (nobjects != objorig+fnobjects)
119 <                                octerror(USER, "bad object count; octree stale?");
120 <                                                /* check for non-surfaces */
121 <                        if (nonsurfinset(objorig, fnobjects))
111 <                                octerror(USER, "non-surface in set; octree stale?");
112 <                }
106 >        else if (load & IO_SCENE && nf == 0)
107 >                skiptree();
108 >                
109 >        if (load & IO_SCENE) {          /* get the scene */
110 >            if (nf == 0) {
111 >                                        /* load binary scene data */
112 >                readscene(infp, objsize);
113 >
114 >            } else {                    /* consistency checks */
115 >                                /* check object count */
116 >                if (nobjects != objorig+fnobjects)
117 >                        octerror(USER, "bad object count; octree stale?");
118 >                                /* check for non-surfaces */
119 >                if (nonsurfintree(scene->cutree))
120 >                        octerror(USER, "modifier in tree; octree stale?");
121 >            }
122          }
123 <        fclose(infp);
123 >                                /* close the input */
124 >        if (infn[0] == '!')
125 >                pclose(infp);
126 >        else
127 >                fclose(infp);
128          return(nf);
129   }
130  
131  
132   static char *
133 < getstr(s)                       /* get null-terminated string */
121 < char  *s;
133 > ogetstr(char *s)                /* get null-terminated string */
134   {
135 <        register char  *cp;
124 <        register int  c;
135 >        extern char  *getstr();
136  
137 <        cp = s;
138 <        while ((c = getc(infp)) != EOF)
139 <                if ((*cp++ = c) == '\0')
129 <                        return(s);
130 <                
131 <        octerror(USER, "truncated octree");
137 >        if (getstr(s, infp) == NULL)
138 >                octerror(USER, "truncated octree");
139 >        return(s);
140   }
141  
142  
143   static OCTREE
144   getfullnode()                   /* get a set, return fullnode */
145   {
146 <        OBJECT  set[MAXSET+1];
146 >        OBJECT  set[MAXSET+1];
147          register int  i;
148          register long  m;
149  
150 <        if ((set[0] = getint(objsize)) > MAXSET)
150 >        if ((set[0] = ogetint(objsize)) > MAXSET)
151                  octerror(USER, "bad set in getfullnode");
152          for (i = 1; i <= set[0]; i++) {
153 <                m = getint(objsize) + objorig;
153 >                m = ogetint(objsize) + objorig;
154                  if ((set[i] = m) != m)
155                          octerror(USER, "too many objects");
156          }
# Line 151 | Line 159 | getfullnode()                  /* get a set, return fullnode */
159  
160  
161   static long
162 < getint(siz)                     /* get a siz-byte integer */
155 < register int  siz;
162 > ogetint(int siz)                /* get a siz-byte integer */
163   {
164 <        register int  c;
164 >        extern long  getint();
165          register long  r;
166  
167 <        if ((c = getc(infp)) == EOF)
168 <                goto end_file;
169 <        r = 0x80&c ? -1<<8|c : c;               /* sign extend */
163 <        while (--siz > 0) {
164 <                if ((c = getc(infp)) == EOF)
165 <                        goto end_file;
166 <                r <<= 8;
167 <                r |= c;
168 <        }
167 >        r = getint(siz, infp);
168 >        if (feof(infp))
169 >                octerror(USER, "truncated octree");
170          return(r);
170 end_file:
171        octerror(USER, "truncated octree");
171   }
172  
173  
174   static double
175 < getflt()                        /* get a floating point number */
175 > ogetflt()                       /* get a floating point number */
176   {
177 <        extern double  ldexp();
178 <        double  d;
177 >        extern double  getflt();
178 >        double  r;
179  
180 <        d = (double)getint(4)/0x7fffffff;
181 <        return(ldexp(d, (int)getint(1)));
180 >        r = getflt(infp);
181 >        if (feof(infp))
182 >                octerror(USER, "truncated octree");
183 >        return(r);
184   }
185          
186  
187   static OCTREE
188   gettree()                       /* get a pre-ordered octree */
189   {
190 <        register OCTREE  ot;
190 >        register OCTREE  ot;
191          register int  i;
192          
193          switch (getc(infp)) {
# Line 205 | Line 206 | gettree()                      /* get a pre-ordered octree */
206          default:
207                  octerror(USER, "damaged octree");
208          }
209 +        return EMPTY; /* pro forma return */
210   }
211  
212  
213 < static
214 < getobj()                                /* get next object */
213 > static int
214 > nonsurfintree(OCTREE ot)                /* check tree for modifiers */
215   {
216 <        char  sbuf[MAXSTR];
215 <        int  obj;
216 >        OBJECT  set[MAXSET+1];
217          register int  i;
217        register long  m;
218        register OBJREC  *objp;
218  
219 <        i = getint(1);
220 <        if (i == -1)
221 <                return(OVOID);          /* terminator */
222 <        if ((obj = newobject()) == OVOID)
223 <                error(SYSTEM, "out of object space");
224 <        objp = objptr(obj);
225 <        if ((objp->otype = otypmap[i]) < 0)
227 <                octerror(USER, "reference to unknown type");
228 <        if ((m = getint(objsize)) != OVOID) {
229 <                m += objorig;
230 <                if ((OBJECT)m != m)
231 <                        octerror(USER, "too many objects");
219 >        if (isempty(ot))
220 >                return(0);
221 >        if (istree(ot)) {
222 >                for (i = 0; i < 8; i++)
223 >                        if (nonsurfintree(octkid(ot, i)))
224 >                                return(1);
225 >                return(0);
226          }
227 <        objp->omod = m;
228 <        objp->oname = savqstr(getstr(sbuf));
229 <        if (objp->oargs.nsargs = getint(2)) {
230 <                objp->oargs.sarg = (char **)bmalloc
231 <                                (objp->oargs.nsargs*sizeof(char *));
238 <                if (objp->oargs.sarg == NULL)
239 <                        goto memerr;
240 <                for (i = 0; i < objp->oargs.nsargs; i++)
241 <                        objp->oargs.sarg[i] = savestr(getstr(sbuf));
242 <        } else
243 <                objp->oargs.sarg = NULL;
244 < #ifdef  IARGS
245 <        if (objp->oargs.niargs = getint(2)) {
246 <                objp->oargs.iarg = (long *)bmalloc
247 <                                (objp->oargs.niargs*sizeof(long));
248 <                if (objp->oargs.iarg == NULL)
249 <                        goto memerr;
250 <                for (i = 0; i < objp->oargs.niargs; i++)
251 <                        objp->oargs.iarg[i] = getint(4);
252 <        } else
253 <                objp->oargs.iarg = NULL;
254 < #endif
255 <        if (objp->oargs.nfargs = getint(2)) {
256 <                objp->oargs.farg = (FLOAT *)bmalloc
257 <                                (objp->oargs.nfargs*sizeof(FLOAT));
258 <                if (objp->oargs.farg == NULL)
259 <                        goto memerr;
260 <                for (i = 0; i < objp->oargs.nfargs; i++)
261 <                        objp->oargs.farg[i] = getflt();
262 <        } else
263 <                objp->oargs.farg = NULL;
264 <                                                /* initialize */
265 <        objp->os = NULL;
266 <        objp->lastrno = -1;
267 <                                                /* insert */
268 <        insertobject(obj);
269 <        return(obj);
270 < memerr:
271 <        error(SYSTEM, "out of memory in getobj");
227 >        objset(set, ot);
228 >        for (i = set[0]; i-- > 0; )
229 >                if (ismodifier(objptr(set[i])->otype))
230 >                        return(1);
231 >        return(0);
232   }
233  
234  
235 < static
236 < octerror(etyp, msg)                     /* octree error */
237 < int  etyp;
238 < char  *msg;
235 > static void
236 > skiptree(void)                          /* skip octree on input */
237 > {
238 >        register int  i;
239 >        
240 >        switch (getc(infp)) {
241 >        case OT_EMPTY:
242 >                return;
243 >        case OT_FULL:
244 >                for (i = ogetint(objsize)*objsize; i-- > 0; )
245 >                        if (getc(infp) == EOF)
246 >                                octerror(USER, "truncated octree");
247 >                return;
248 >        case OT_TREE:
249 >                for (i = 0; i < 8; i++)
250 >                        skiptree();
251 >                return;
252 >        case EOF:
253 >                octerror(USER, "truncated octree");
254 >        default:
255 >                octerror(USER, "damaged octree");
256 >        }
257 > }
258 >
259 >
260 > static void
261 > octerror(int etyp, char *msg)           /* octree error */
262   {
263          char  msgbuf[128];
264  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines