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

Comparing ray/src/rt/raycalls.c (file contents):
Revision 2.17 by greg, Sat Dec 12 05:20:10 2009 UTC vs.
Revision 2.33 by greg, Tue Jun 3 21:31:51 2025 UTC

# Line 36 | Line 36 | static const char      RCSid[] = "$Id$";
36   *  as simple as we can.  Global variables and their defaults
37   *  are defined below, and including "ray.h" declares these
38   *  along with all the routines you are likely to need.  First,
39 < *  assign the global variable progname to your argv[0], then
40 < *  change the rendering parameters as you like.  If you have a set
41 < *  of option arguments you are working from, the getrenderopt(ac,av)
39 > *  assign the global variable progname by calling fixargv0(argv[0]),
40 > *  then change the rendering parameters as you like.  If you have a
41 > *  set of option arguments you are working from, the getrenderopt()
42   *  call should be very useful.  Before tracing any rays, you
43   *  must read in the octree with a call to ray_init(oct).
44   *  Passing NULL for the file name causes ray_init() to read
# Line 94 | Line 94 | static const char      RCSid[] = "$Id$";
94  
95   #include  "ray.h"
96   #include  "source.h"
97 + #include  "bsdf.h"
98   #include  "ambient.h"
99   #include  "otypes.h"
100   #include  "random.h"
101 + #include  "func.h"
102   #include  "data.h"
103   #include  "font.h"
104 + #include  "pmapray.h"
105  
103 – char    *progname = "unknown_app";      /* caller sets to argv[0] */
104 –
106   char    *octname;                       /* octree name we are given */
107  
107 – char    *shm_boundary = NULL;           /* boundary of shared memory */
108 –
108   CUBE    thescene;                       /* our scene */
109   OBJECT  nsceneobjs;                     /* number of objects in our scene */
110  
# Line 117 | Line 116 | void   (*trace)() = NULL;              /* trace call */
116  
117   void    (*addobjnotify[8])() = {ambnotify, NULL};
118  
119 + int     castonly = 0;                   /* only doing ray-casting? */
120 +
121   int     do_irrad = 0;                   /* compute irradiance? */
122  
123 < int     rand_samp = 0;                  /* pure Monte Carlo sampling? */
123 > int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
124  
125   double  dstrsrc = 0.0;                  /* square source distribution */
126   double  shadthresh = .03;               /* shadow threshold */
# Line 139 | Line 140 | double specjitter = 1.;                /* specular sampling jitter *
140  
141   int     backvis = 1;                    /* back face visibility */
142  
143 < int     maxdepth = 8;                   /* maximum recursion depth */
144 < double  minweight = 2e-3;               /* minimum ray weight */
143 > int     maxdepth = -10;                 /* maximum recursion depth */
144 > double  minweight = 1e-4;               /* minimum ray weight */
145  
146   char    *ambfile = NULL;                /* ambient file name */
147   COLOR   ambval = BLKCOLOR;              /* ambient value */
# Line 154 | Line 155 | char   *amblist[AMBLLEN+1];            /* ambient include/exclude
155   int     ambincl = -1;                   /* include == 1, exclude == 0 */
156  
157  
158 + static void
159 + reset_random(void)              /* re-initialize random number generator */
160 + {
161 +        if (rand_samp) {
162 +                srandom((long)time(0));
163 +                initurand(0);
164 +        } else {
165 +                srandom(0L);
166 +                initurand(2048);
167 +        }
168 + }
169 +
170 +
171   void
172   ray_init(                       /* initialize ray-tracing calculation */
173          char    *otnm
# Line 164 | Line 178 | ray_init(                      /* initialize ray-tracing calculation */
178                                          /* initialize object types */
179          if (ofun[OBJ_SPHERE].funp == o_default)
180                  initotypes();
181 +                                        /* initialize calcomp routines */
182 +        initfunc();
183                                          /* initialize urand */
184 <        if (rand_samp) {
185 <                srandom((long)time(0));
186 <                initurand(0);
187 <        } else {
188 <                srandom(0L);
189 <                initurand(2048);
190 <        }
175 <                                        /* read scene octree */
176 <        readoct(octname = otnm, ~(IO_FILES|IO_INFO), &thescene, NULL);
184 >        reset_random();
185 >                                        /* initialize spectral sampling */
186 >        if (setspectrsamp(CNDX, WLPART) < 0)
187 >                error(USER, "unsupported spectral sampling");
188 >
189 >        octname = savqstr(otnm);        /* read scene octree */
190 >        readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
191          nsceneobjs = nobjects;
192 <                                        /* find and mark sources */
193 <        marksources();
194 <                                        /* initialize ambient calculation */
195 <        setambient();
196 <                                        /* ready to go... */
192 >
193 >        if (!castonly) {                /* any actual ray traversal to do? */
194 >
195 >                ray_init_pmap();        /* PMAP: set up & load photon maps */
196 >
197 >                marksources();          /* find and mark sources */
198 >
199 >                setambient();           /* initialize ambient calculation */
200 >        } else
201 >                distantsources();       /* else mark only distant sources */
202   }
203  
204 +
205   void
206   ray_trace(                      /* trace a primary ray */
207          RAY     *r
# Line 206 | Line 226 | ray_done(              /* free ray-tracing data */
226          donesets();
227          octdone();
228          thescene.cutree = EMPTY;
229 <        octname = NULL;
229 >        freeqstr(octname); octname = NULL;
230 >        retainfonts = 0;
231          if (freall) {
211 –                retainfonts = 0;
232                  freefont(NULL);
233                  freedata(NULL);
234 +                SDfreeCache(NULL);
235                  initurand(0);
236          }
237          if (nobjects > 0) {
238                  sprintf(errmsg, "%ld objects left after call to ray_done()",
239 <                                nobjects);
239 >                                (long)nobjects);
240                  error(WARNING, errmsg);
241          }
242 +        
243 +        ray_done_pmap();
244   }
245  
246  
# Line 231 | Line 254 | ray_save(                      /* save current parameter settings */
254          if (rp == NULL)
255                  return;
256          rp->do_irrad = do_irrad;
257 +        rp->rand_samp = rand_samp;
258          rp->dstrsrc = dstrsrc;
259          rp->shadthresh = shadthresh;
260          rp->shadcert = shadcert;
# Line 247 | Line 271 | ray_save(                      /* save current parameter settings */
271          rp->backvis = backvis;
272          rp->maxdepth = maxdepth;
273          rp->minweight = minweight;
250 –        copycolor(rp->ambval, ambval);
251 –        memset(rp->ambfile, '\0', sizeof(rp->ambfile));
274          if (ambfile != NULL)
275                  strncpy(rp->ambfile, ambfile, sizeof(rp->ambfile)-1);
276 +        else
277 +                memset(rp->ambfile, '\0', sizeof(rp->ambfile));
278 +        copycolor(rp->ambval, ambval);
279          rp->ambvwt = ambvwt;
280          rp->ambacc = ambacc;
281          rp->ambres = ambres;
# Line 265 | Line 290 | ray_save(                      /* save current parameter settings */
290                  if (ndx+len >= sizeof(rp->amblval))
291                          break;
292                  strcpy(rp->amblval+ndx, amblist[i]);
293 +                rp->amblndx[i] = ndx;
294                  ndx += len+1;
295          }
296          while (i <= AMBLLEN)
297                  rp->amblndx[i++] = -1;
298 +                
299 +        /* PMAP: save photon mapping params */
300 +        ray_save_pmap(rp);
301   }
302  
303  
# Line 277 | Line 306 | ray_restore(                   /* restore parameter settings */
306          RAYPARAMS       *rp
307   )
308   {
309 <        register int    i;
309 >        int     i;
310  
311          if (rp == NULL) {               /* restore defaults */
312                  RAYPARAMS       dflt;
# Line 287 | Line 316 | ray_restore(                   /* restore parameter settings */
316          }
317                                          /* restore saved settings */
318          do_irrad = rp->do_irrad;
319 +        if (!rand_samp != !rp->rand_samp) {
320 +                rand_samp = rp->rand_samp;
321 +                reset_random();
322 +        }
323          dstrsrc = rp->dstrsrc;
324          shadthresh = rp->shadthresh;
325          shadcert = rp->shadcert;
# Line 308 | Line 341 | ray_restore(                   /* restore parameter settings */
341          ambdiv = rp->ambdiv;
342          ambssamp = rp->ambssamp;
343          ambounce = rp->ambounce;
344 +                                        /* a bit dangerous if not static */
345          for (i = 0; rp->amblndx[i] >= 0; i++)
346                  amblist[i] = rp->amblval + rp->amblndx[i];
347          while (i <= AMBLLEN)
# Line 315 | Line 349 | ray_restore(                   /* restore parameter settings */
349          ambincl = rp->ambincl;
350                                          /* update ambient calculation */
351          ambnotify(OVOID);
352 <        if (thescene.cutree != EMPTY) {
352 >        if ((thescene.cutree != EMPTY) & !castonly) {
353                  int     newamb = (ambfile == NULL) ?  rp->ambfile[0] :
354                                          strcmp(ambfile, rp->ambfile) ;
355  
# Line 337 | Line 371 | ray_restore(                   /* restore parameter settings */
371                  ambres = rp->ambres;
372                  ambacc = rp->ambacc;
373          }
374 +        
375 +        /* PMAP: restore photon mapping params */
376 +        ray_restore_pmap(rp);
377   }
378  
379  
# Line 351 | Line 388 | ray_defaults(          /* get default parameter values */
388                  return;
389  
390          rp->do_irrad = 0;
391 +        rp->rand_samp = 1;
392          rp->dstrsrc = 0.0;
393 <        rp->shadthresh = .03;
394 <        rp->shadcert = .75;
393 >        rp->shadthresh = 0.03;
394 >        rp->shadcert = 0.75;
395          rp->directrelay = 2;
396          rp->vspretest = 512;
397          rp->directvis = 1;
# Line 362 | Line 400 | ray_defaults(          /* get default parameter values */
400          setcolor(rp->salbedo, 0., 0., 0.);
401          rp->seccg = 0.;
402          rp->ssampdist = 0.;
403 <        rp->specthresh = .15;
403 >        rp->specthresh = 0.15;
404          rp->specjitter = 1.;
405          rp->backvis = 1;
406 <        rp->maxdepth = 8;
407 <        rp->minweight = 2e-3;
370 <        setcolor(rp->ambval, 0., 0., 0.);
406 >        rp->maxdepth = -10;
407 >        rp->minweight = 1e-4;
408          memset(rp->ambfile, '\0', sizeof(rp->ambfile));
409 +        setcolor(rp->ambval, 0., 0., 0.);
410          rp->ambvwt = 0;
411          rp->ambres = 256;
412 <        rp->ambacc = 0.15;
412 >        rp->ambacc = 0.1;
413          rp->ambdiv = 1024;
414          rp->ambssamp = 512;
415          rp->ambounce = 0;
# Line 379 | Line 417 | ray_defaults(          /* get default parameter values */
417          memset(rp->amblval, '\0', sizeof(rp->amblval));
418          for (i = AMBLLEN+1; i--; )
419                  rp->amblndx[i] = -1;
420 +        
421 +        /* PMAP: restore photon mapping defaults */
422 +        ray_defaults_pmap(rp);
423   }

Diff Legend

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