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

Comparing ray/src/ot/oconv.c (file contents):
Revision 2.2 by greg, Thu Dec 19 15:06:31 1991 UTC vs.
Revision 2.13 by greg, Wed Mar 12 04:59:04 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   *  oconv.c - main program for object to octree conversion.
6   *
# Line 18 | Line 15 | static char SCCSid[] = "$SunId$ LBL";
15  
16   #include  "otypes.h"
17  
18 < #ifndef  DEFPATH
22 < #define  DEFPATH        ":/usr/local/lib/ray"
23 < #endif
18 > #include  "paths.h"
19  
20 < #define  OMARGIN        (10*FTINY)      /* margin around global cube */
20 > #define  OMARGIN        (10*FTINY)      /* margin around global cube */
21  
22 < #define  MAXOBJFIL      63              /* maximum number of scene files */
22 > #define  MAXOBJFIL      127             /* maximum number of scene files */
23  
24   char  *progname;                        /* argv[0] */
25  
31 – char  *libpath;                         /* library search path */
32 –
26   int  nowarn = 0;                        /* supress warnings? */
27  
28   int  objlim = 5;                        /* # of objects before split */
29  
30 < int  resolu = 1024;                     /* octree resolution limit */
30 > int  resolu = 8192;                     /* octree resolution limit */
31  
32   CUBE  thescene = {EMPTY, {0.0, 0.0, 0.0}, 0.0};         /* our scene */
33  
34   char  *ofname[MAXOBJFIL+1];             /* object file names */
35   int  nfiles = 0;                        /* number of object files */
36  
37 < double  mincusize;                      /* minimum cube size from resolu */
37 > double  mincusize;                      /* minimum cube size from resolu */
38  
39 < int  (*addobjnotify[])() = {NULL};      /* new object notifier functions */
39 > void  (*addobjnotify[])() = {NULL};     /* new object notifier functions */
40  
41  
42   main(argc, argv)                /* convert object files to an octree */
43   int  argc;
44 < char  **argv;
44 > char  *argv[];
45   {
53 –        extern char  *getenv();
46          FVECT  bbmin, bbmax;
47          char  *infile = NULL;
48 +        int  inpfrozen = 0;
49          int  outflags = IO_ALL;
50 <        OBJECT  startobj;
50 >        OBJECT  startobj;
51          int  i;
52  
53 <        progname = argv[0];
53 >        progname = argv[0] = fixargv0(argv[0]);
54  
62 –        if ((libpath = getenv("RAYPATH")) == NULL)
63 –                libpath = DEFPATH;
64 –
55          initotypes();
56  
57          for (i = 1; i < argc && argv[i][0] == '-'; i++)
# Line 95 | Line 85 | char  **argv;
85                          break;
86                  }
87   breakopt:
88 + #ifdef MSDOS
89 +        setmode(fileno(stdout), O_BINARY);
90 + #endif
91          if (infile != NULL) {           /* get old octree & objects */
92                  if (thescene.cusize > FTINY)
93                          error(USER, "only one of '-b' or '-i'");
94                  nfiles = readoct(infile, IO_ALL, &thescene, ofname);
95 <                if (nfiles == 0 && outflags & IO_FILES) {
96 <                        error(WARNING, "frozen octree");
97 <                        outflags &= ~IO_FILES;
98 <                }
99 <        }
107 <
108 <        printargs(argc, argv, stdout);  /* info. header */
95 >                if (nfiles == 0)
96 >                        inpfrozen++;
97 >        } else
98 >                newheader("RADIANCE", stdout);  /* new binary file header */
99 >        printargs(argc, argv, stdout);
100          fputformat(OCTFMT, stdout);
101          printf("\n");
102  
103          startobj = nobjects;            /* previous objects already converted */
104 <        
104 >
105          for ( ; i < argc; i++)          /* read new scene descriptions */
106                  if (!strcmp(argv[i], "-")) {    /* from stdin */
107                          readobj(NULL);
# Line 122 | Line 113 | breakopt:
113                  }
114  
115          ofname[nfiles] = NULL;
116 +
117 +        if (inpfrozen && outflags & IO_FILES) {
118 +                error(WARNING, "frozen octree");
119 +                outflags &= ~IO_FILES;
120 +        }
121                                                  /* find bounding box */
122          bbmin[0] = bbmin[1] = bbmin[2] = FHUGE;
123          bbmax[0] = bbmax[1] = bbmax[2] = -FHUGE;
# Line 149 | Line 145 | breakopt:
145          }
146  
147          mincusize = thescene.cusize / resolu - FTINY;
148 <                
148 >
149          for (i = startobj; i < nobjects; i++)           /* add new objects */
150                  addobject(&thescene, i);
151 <        
151 >
152          thescene.cutree = combine(thescene.cutree);     /* optimize */
153  
154          writeoct(outflags, &thescene, ofname);  /* write structures to stdout */
# Line 161 | Line 157 | breakopt:
157   }
158  
159  
160 + void
161   quit(code)                              /* exit program */
162   int  code;
163   {
# Line 168 | Line 165 | int  code;
165   }
166  
167  
168 + void
169   cputs()                                 /* interactive error */
170   {
171          /* referenced, but not used */
172   }
173  
174  
175 + void
176   wputs(s)                                /* warning message */
177   char  *s;
178   {
# Line 182 | Line 181 | char  *s;
181   }
182  
183  
184 + void
185   eputs(s)                                /* put string to stderr */
186   register char  *s;
187   {
# Line 197 | Line 197 | register char  *s;
197   }
198  
199  
200 < #define  bitop(f,i,op)          (f[((i)>>3)] op (1<<((i)&7)))
201 < #define  tstbit(f,i)            bitop(f,i,&)
202 < #define  setbit(f,i)            bitop(f,i,|=)
203 < #define  clrbit(f,i)            bitop(f,i,&=~)
204 < #define  tglbit(f,i)            bitop(f,i,^=)
200 > #define  bitop(f,i,op)          (f[((i)>>3)] op (1<<((i)&7)))
201 > #define  tstbit(f,i)            bitop(f,i,&)
202 > #define  setbit(f,i)            bitop(f,i,|=)
203 > #define  clrbit(f,i)            bitop(f,i,&=~)
204 > #define  tglbit(f,i)            bitop(f,i,^=)
205  
206  
207   addobject(cu, obj)                      /* add an object to a cube */
208   register CUBE  *cu;
209 < OBJECT  obj;
209 > OBJECT  obj;
210   {
211 <        CUBE  cukid;
212 <        OCTREE  ot;
213 <        OBJECT  oset[MAXSET+1];
214 <        unsigned char  inflg[MAXSET/8], volflg[MAXSET/8];
215 <        int  in;
216 <        register int  i, j;
211 >        int  inc;
212  
213 <        in = (*ofun[objptr(obj)->otype].funp)(objptr(obj), cu);
213 >        inc = (*ofun[objptr(obj)->otype].funp)(objptr(obj), cu);
214  
215 <        if (in == O_MISS)
215 >        if (inc == O_MISS)
216                  return;                         /* no intersection */
217 <        
217 >
218          if (istree(cu->cutree)) {
219 <                                                /* do children */
219 >                CUBE  cukid;                    /* do children */
220 >                int  i, j;
221                  cukid.cusize = cu->cusize * 0.5;
222                  for (i = 0; i < 8; i++) {
223                          cukid.cutree = octkid(cu->cutree, i);
# Line 236 | Line 232 | OBJECT  obj;
232                  return;
233          }
234          if (isempty(cu->cutree)) {
235 <                                                /* singular set */
235 >                OBJECT  oset[2];                /* singular set */
236                  oset[0] = 1; oset[1] = obj;
237                  cu->cutree = fullnode(oset);
238                  return;
239          }
240                                          /* add to full node */
241 +        add2full(cu, obj, inc);
242 + }
243 +
244 +
245 + add2full(cu, obj, inc)                  /* add object to full node */
246 + register CUBE  *cu;
247 + OBJECT  obj;
248 + int  inc;
249 + {
250 +        OCTREE  ot;
251 +        OBJECT  oset[MAXSET+1];
252 +        CUBE  cukid;
253 +        unsigned char  inflg[(MAXSET+7)/8], volflg[(MAXSET+7)/8];
254 +        register int  i, j;
255 +
256          objset(oset, cu->cutree);
257          cukid.cusize = cu->cusize * 0.5;
258 <        
259 <        if (in==O_IN || oset[0] < objlim || cukid.cusize < mincusize) {
258 >
259 >        if (inc==O_IN || oset[0] < objlim || cukid.cusize < mincusize) {
260                                                  /* add to set */
261                  if (oset[0] >= MAXSET) {
262 <                        sprintf(errmsg,
252 <                                "set overflow in addobject (%s)",
262 >                        sprintf(errmsg, "set overflow in addobject (%s)",
263                                          objptr(obj)->oname);
264                          error(INTERNAL, errmsg);
265                  }
# Line 268 | Line 278 | OBJECT  obj;
278                  if (isvolume(objptr(oset[j])->otype)) {
279                          setbit(volflg,j-1);
280                          if ((*ofun[objptr(oset[j])->otype].funp)
281 <                                        (objptr(oset[j]),cu) == O_IN)
281 >                                        (objptr(oset[j]), cu) == O_IN)
282                                  setbit(inflg,j-1);
283                  }
284                                          /* assign subcubes */
# Line 285 | Line 295 | OBJECT  obj;
295                                  addobject(&cukid, oset[j]);
296                                          /* then this object */
297                  addobject(&cukid, obj);
298 <                                        /* partial volumes */
298 >                                        /* then partial volumes */
299                  for (j = 1; j <= oset[0]; j++)
300                          if (tstbit(volflg,j-1) &&
301                                          !tstbit(inflg,j-1))
302                                  addobject(&cukid, oset[j]);
303 <                                        /* full volumes */
303 >                                        /* full volumes last */
304                  for (j = 1; j <= oset[0]; j++)
305                          if (tstbit(inflg,j-1))
306                                  addobject(&cukid, oset[j]);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines