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.20 by greg, Fri Feb 25 06:51:50 2011 UTC vs.
Revision 2.30 by greg, Thu May 2 15:02:08 2024 UTC

# Line 100 | Line 100 | static const char      RCSid[] = "$Id$";
100   #include  "random.h"
101   #include  "data.h"
102   #include  "font.h"
103 + #include  "pmapray.h"
104  
105   char    *progname = "unknown_app";      /* caller sets to argv[0] */
106  
# Line 118 | Line 119 | void   (*trace)() = NULL;              /* trace call */
119  
120   void    (*addobjnotify[8])() = {ambnotify, NULL};
121  
122 + int     castonly = 0;                   /* only doing ray-casting? */
123 +
124   int     do_irrad = 0;                   /* compute irradiance? */
125  
126 < int     rand_samp = 0;                  /* pure Monte Carlo sampling? */
126 > int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
127  
128   double  dstrsrc = 0.0;                  /* square source distribution */
129   double  shadthresh = .03;               /* shadow threshold */
# Line 140 | Line 143 | double specjitter = 1.;                /* specular sampling jitter *
143  
144   int     backvis = 1;                    /* back face visibility */
145  
146 < int     maxdepth = 8;                   /* maximum recursion depth */
147 < double  minweight = 5e-4;               /* minimum ray weight */
146 > int     maxdepth = -10;                 /* maximum recursion depth */
147 > double  minweight = 1e-4;               /* minimum ray weight */
148  
149   char    *ambfile = NULL;                /* ambient file name */
150   COLOR   ambval = BLKCOLOR;              /* ambient value */
# Line 155 | Line 158 | char   *amblist[AMBLLEN+1];            /* ambient include/exclude
158   int     ambincl = -1;                   /* include == 1, exclude == 0 */
159  
160  
161 + static void
162 + reset_random(void)              /* re-initialize random number generator */
163 + {
164 +        if (rand_samp) {
165 +                srandom((long)time(0));
166 +                initurand(0);
167 +        } else {
168 +                srandom(0L);
169 +                initurand(2048);
170 +        }
171 + }
172 +
173 +
174   void
175   ray_init(                       /* initialize ray-tracing calculation */
176          char    *otnm
# Line 166 | Line 182 | ray_init(                      /* initialize ray-tracing calculation */
182          if (ofun[OBJ_SPHERE].funp == o_default)
183                  initotypes();
184                                          /* initialize urand */
185 <        if (rand_samp) {
186 <                srandom((long)time(0));
187 <                initurand(0);
188 <        } else {
189 <                srandom(0L);
190 <                initurand(2048);
191 <        }
176 <                                        /* read scene octree */
177 <        readoct(octname = otnm, ~(IO_FILES|IO_INFO), &thescene, NULL);
185 >        reset_random();
186 >                                        /* initialize spectral sampling */
187 >        if (setspectrsamp(CNDX, WLPART) < 0)
188 >                error(USER, "unsupported spectral sampling");
189 >
190 >        octname = savqstr(otnm);        /* read scene octree */
191 >        readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
192          nsceneobjs = nobjects;
193 <                                        /* find and mark sources */
194 <        marksources();
195 <                                        /* initialize ambient calculation */
196 <        setambient();
197 <                                        /* ready to go... */
193 >
194 >        if (!castonly) {                /* any actual ray traversal to do? */
195 >
196 >                ray_init_pmap();        /* PMAP: set up & load photon maps */
197 >
198 >                marksources();          /* find and mark sources */
199 >
200 >                setambient();           /* initialize ambient calculation */
201 >        } else
202 >                distantsources();       /* else mark only distant sources */
203   }
204  
205 +
206   void
207   ray_trace(                      /* trace a primary ray */
208          RAY     *r
# Line 207 | Line 227 | ray_done(              /* free ray-tracing data */
227          donesets();
228          octdone();
229          thescene.cutree = EMPTY;
230 <        octname = NULL;
230 >        freeqstr(octname); octname = NULL;
231          retainfonts = 0;
232          if (freall) {
233                  freefont(NULL);
# Line 220 | Line 240 | ray_done(              /* free ray-tracing data */
240                                  (long)nobjects);
241                  error(WARNING, errmsg);
242          }
243 +        
244 +        ray_done_pmap();
245   }
246  
247  
# Line 233 | Line 255 | ray_save(                      /* save current parameter settings */
255          if (rp == NULL)
256                  return;
257          rp->do_irrad = do_irrad;
258 +        rp->rand_samp = rand_samp;
259          rp->dstrsrc = dstrsrc;
260          rp->shadthresh = shadthresh;
261          rp->shadcert = shadcert;
# Line 249 | Line 272 | ray_save(                      /* save current parameter settings */
272          rp->backvis = backvis;
273          rp->maxdepth = maxdepth;
274          rp->minweight = minweight;
252 –        copycolor(rp->ambval, ambval);
253 –        memset(rp->ambfile, '\0', sizeof(rp->ambfile));
275          if (ambfile != NULL)
276                  strncpy(rp->ambfile, ambfile, sizeof(rp->ambfile)-1);
277 +        else
278 +                memset(rp->ambfile, '\0', sizeof(rp->ambfile));
279 +        copycolor(rp->ambval, ambval);
280          rp->ambvwt = ambvwt;
281          rp->ambacc = ambacc;
282          rp->ambres = ambres;
# Line 267 | Line 291 | ray_save(                      /* save current parameter settings */
291                  if (ndx+len >= sizeof(rp->amblval))
292                          break;
293                  strcpy(rp->amblval+ndx, amblist[i]);
294 +                rp->amblndx[i] = ndx;
295                  ndx += len+1;
296          }
297          while (i <= AMBLLEN)
298                  rp->amblndx[i++] = -1;
299 +                
300 +        /* PMAP: save photon mapping params */
301 +        ray_save_pmap(rp);
302   }
303  
304  
# Line 279 | Line 307 | ray_restore(                   /* restore parameter settings */
307          RAYPARAMS       *rp
308   )
309   {
310 <        register int    i;
310 >        int     i;
311  
312          if (rp == NULL) {               /* restore defaults */
313                  RAYPARAMS       dflt;
# Line 289 | Line 317 | ray_restore(                   /* restore parameter settings */
317          }
318                                          /* restore saved settings */
319          do_irrad = rp->do_irrad;
320 +        if (!rand_samp != !rp->rand_samp) {
321 +                rand_samp = rp->rand_samp;
322 +                reset_random();
323 +        }
324          dstrsrc = rp->dstrsrc;
325          shadthresh = rp->shadthresh;
326          shadcert = rp->shadcert;
# Line 310 | Line 342 | ray_restore(                   /* restore parameter settings */
342          ambdiv = rp->ambdiv;
343          ambssamp = rp->ambssamp;
344          ambounce = rp->ambounce;
345 +                                        /* a bit dangerous if not static */
346          for (i = 0; rp->amblndx[i] >= 0; i++)
347                  amblist[i] = rp->amblval + rp->amblndx[i];
348          while (i <= AMBLLEN)
# Line 317 | Line 350 | ray_restore(                   /* restore parameter settings */
350          ambincl = rp->ambincl;
351                                          /* update ambient calculation */
352          ambnotify(OVOID);
353 <        if (thescene.cutree != EMPTY) {
353 >        if ((thescene.cutree != EMPTY) & !castonly) {
354                  int     newamb = (ambfile == NULL) ?  rp->ambfile[0] :
355                                          strcmp(ambfile, rp->ambfile) ;
356  
# Line 339 | Line 372 | ray_restore(                   /* restore parameter settings */
372                  ambres = rp->ambres;
373                  ambacc = rp->ambacc;
374          }
375 +        
376 +        /* PMAP: restore photon mapping params */
377 +        ray_restore_pmap(rp);
378   }
379  
380  
# Line 353 | Line 389 | ray_defaults(          /* get default parameter values */
389                  return;
390  
391          rp->do_irrad = 0;
392 +        rp->rand_samp = 1;
393          rp->dstrsrc = 0.0;
394 <        rp->shadthresh = .03;
395 <        rp->shadcert = .75;
394 >        rp->shadthresh = 0.03;
395 >        rp->shadcert = 0.75;
396          rp->directrelay = 2;
397          rp->vspretest = 512;
398          rp->directvis = 1;
# Line 364 | Line 401 | ray_defaults(          /* get default parameter values */
401          setcolor(rp->salbedo, 0., 0., 0.);
402          rp->seccg = 0.;
403          rp->ssampdist = 0.;
404 <        rp->specthresh = .15;
404 >        rp->specthresh = 0.15;
405          rp->specjitter = 1.;
406          rp->backvis = 1;
407 <        rp->maxdepth = 8;
408 <        rp->minweight = 2e-3;
372 <        setcolor(rp->ambval, 0., 0., 0.);
407 >        rp->maxdepth = -10;
408 >        rp->minweight = 1e-4;
409          memset(rp->ambfile, '\0', sizeof(rp->ambfile));
410 +        setcolor(rp->ambval, 0., 0., 0.);
411          rp->ambvwt = 0;
412          rp->ambres = 256;
413 <        rp->ambacc = 0.15;
413 >        rp->ambacc = 0.1;
414          rp->ambdiv = 1024;
415          rp->ambssamp = 512;
416          rp->ambounce = 0;
# Line 381 | Line 418 | ray_defaults(          /* get default parameter values */
418          memset(rp->amblval, '\0', sizeof(rp->amblval));
419          for (i = AMBLLEN+1; i--; )
420                  rp->amblndx[i] = -1;
421 +        
422 +        /* PMAP: restore photon mapping defaults */
423 +        ray_defaults_pmap(rp);
424   }

Diff Legend

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