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.15 by greg, Tue Jun 21 00:26:44 2005 UTC vs.
Revision 2.29 by greg, Thu Feb 2 19:13:13 2023 UTC

# 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  "data.h"
102   #include  "font.h"
103 + #include  "pmapray.h"
104  
105   char    *progname = "unknown_app";      /* caller sets to argv[0] */
106  
# Line 117 | 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 139 | 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 = 2e-3;               /* 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 154 | Line 158 | char   *amblist[AMBLLEN+1];            /* ambient include/exclude
158   int     ambincl = -1;                   /* include == 1, exclude == 0 */
159  
160  
161 < extern void
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
177   )
# Line 165 | Line 182 | ray_init(                      /* initialize ray-tracing calculation */
182          if (ofun[OBJ_SPHERE].funp == o_default)
183                  initotypes();
184                                          /* initialize urand */
185 <        srandom(rand_samp ? (long)time(0) : 0L);
186 <        initurand(2048);
187 <                                        /* read scene octree */
188 <        readoct(octname = otnm, ~(IO_FILES|IO_INFO), &thescene, NULL);
185 >        reset_random();
186 >
187 >        octname = savqstr(otnm);        /* read scene octree */
188 >        readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
189          nsceneobjs = nobjects;
190 <                                        /* find and mark sources */
191 <        marksources();
192 <                                        /* initialize ambient calculation */
193 <        setambient();
194 <                                        /* ready to go... */
190 >
191 >        if (!castonly) {                /* any actual ray traversal to do? */
192 >
193 >                ray_init_pmap();        /* PMAP: set up & load photon maps */
194 >
195 >                marksources();          /* find and mark sources */
196 >
197 >                setambient();           /* initialize ambient calculation */
198 >        } else
199 >                distantsources();       /* else mark only distant sources */
200   }
201  
202 < extern void
202 >
203 > void
204   ray_trace(                      /* trace a primary ray */
205          RAY     *r
206   )
207   {
208          rayorigin(r, PRIMARY, NULL, NULL);
209 <        samplendx = rand_samp ? random() : samplendx+1;
209 >        samplendx++;
210          rayvalue(r);            /* assumes origin and direction are set */
211   }
212  
213  
214 < extern void
214 > void
215   ray_done(               /* free ray-tracing data */
216          int     freall
217   )
# Line 201 | Line 224 | ray_done(              /* free ray-tracing data */
224          donesets();
225          octdone();
226          thescene.cutree = EMPTY;
227 <        octname = NULL;
227 >        freeqstr(octname); octname = NULL;
228 >        retainfonts = 0;
229          if (freall) {
206 –                retainfonts = 0;
230                  freefont(NULL);
231                  freedata(NULL);
232 +                SDfreeCache(NULL);
233                  initurand(0);
234          }
235          if (nobjects > 0) {
236                  sprintf(errmsg, "%ld objects left after call to ray_done()",
237 <                                nobjects);
237 >                                (long)nobjects);
238                  error(WARNING, errmsg);
239          }
240 +        
241 +        ray_done_pmap();
242   }
243  
244  
245 < extern void
245 > void
246   ray_save(                       /* save current parameter settings */
247          RAYPARAMS       *rp
248   )
# Line 226 | Line 252 | ray_save(                      /* save current parameter settings */
252          if (rp == NULL)
253                  return;
254          rp->do_irrad = do_irrad;
255 +        rp->rand_samp = rand_samp;
256          rp->dstrsrc = dstrsrc;
257          rp->shadthresh = shadthresh;
258          rp->shadcert = shadcert;
# Line 242 | Line 269 | ray_save(                      /* save current parameter settings */
269          rp->backvis = backvis;
270          rp->maxdepth = maxdepth;
271          rp->minweight = minweight;
245 –        copycolor(rp->ambval, ambval);
246 –        memset(rp->ambfile, '\0', sizeof(rp->ambfile));
272          if (ambfile != NULL)
273                  strncpy(rp->ambfile, ambfile, sizeof(rp->ambfile)-1);
274 +        else
275 +                memset(rp->ambfile, '\0', sizeof(rp->ambfile));
276 +        copycolor(rp->ambval, ambval);
277          rp->ambvwt = ambvwt;
278          rp->ambacc = ambacc;
279          rp->ambres = ambres;
# Line 260 | Line 288 | ray_save(                      /* save current parameter settings */
288                  if (ndx+len >= sizeof(rp->amblval))
289                          break;
290                  strcpy(rp->amblval+ndx, amblist[i]);
291 +                rp->amblndx[i] = ndx;
292                  ndx += len+1;
293          }
294          while (i <= AMBLLEN)
295                  rp->amblndx[i++] = -1;
296 +                
297 +        /* PMAP: save photon mapping params */
298 +        ray_save_pmap(rp);
299   }
300  
301  
302 < extern void
302 > void
303   ray_restore(                    /* restore parameter settings */
304          RAYPARAMS       *rp
305   )
306   {
307 <        register int    i;
307 >        int     i;
308  
309          if (rp == NULL) {               /* restore defaults */
310                  RAYPARAMS       dflt;
# Line 282 | Line 314 | ray_restore(                   /* restore parameter settings */
314          }
315                                          /* restore saved settings */
316          do_irrad = rp->do_irrad;
317 +        if (!rand_samp != !rp->rand_samp) {
318 +                rand_samp = rp->rand_samp;
319 +                reset_random();
320 +        }
321          dstrsrc = rp->dstrsrc;
322          shadthresh = rp->shadthresh;
323          shadcert = rp->shadcert;
# Line 303 | Line 339 | ray_restore(                   /* restore parameter settings */
339          ambdiv = rp->ambdiv;
340          ambssamp = rp->ambssamp;
341          ambounce = rp->ambounce;
342 +                                        /* a bit dangerous if not static */
343          for (i = 0; rp->amblndx[i] >= 0; i++)
344                  amblist[i] = rp->amblval + rp->amblndx[i];
345          while (i <= AMBLLEN)
# Line 310 | Line 347 | ray_restore(                   /* restore parameter settings */
347          ambincl = rp->ambincl;
348                                          /* update ambient calculation */
349          ambnotify(OVOID);
350 <        if (thescene.cutree != EMPTY) {
350 >        if ((thescene.cutree != EMPTY) & !castonly) {
351                  int     newamb = (ambfile == NULL) ?  rp->ambfile[0] :
352                                          strcmp(ambfile, rp->ambfile) ;
353  
# Line 332 | Line 369 | ray_restore(                   /* restore parameter settings */
369                  ambres = rp->ambres;
370                  ambacc = rp->ambacc;
371          }
372 +        
373 +        /* PMAP: restore photon mapping params */
374 +        ray_restore_pmap(rp);
375   }
376  
377  
378 < extern void
378 > void
379   ray_defaults(           /* get default parameter values */
380          RAYPARAMS       *rp
381   )
# Line 346 | Line 386 | ray_defaults(          /* get default parameter values */
386                  return;
387  
388          rp->do_irrad = 0;
389 +        rp->rand_samp = 1;
390          rp->dstrsrc = 0.0;
391 <        rp->shadthresh = .03;
392 <        rp->shadcert = .75;
391 >        rp->shadthresh = 0.03;
392 >        rp->shadcert = 0.75;
393          rp->directrelay = 2;
394          rp->vspretest = 512;
395          rp->directvis = 1;
# Line 357 | Line 398 | ray_defaults(          /* get default parameter values */
398          setcolor(rp->salbedo, 0., 0., 0.);
399          rp->seccg = 0.;
400          rp->ssampdist = 0.;
401 <        rp->specthresh = .15;
401 >        rp->specthresh = 0.15;
402          rp->specjitter = 1.;
403          rp->backvis = 1;
404 <        rp->maxdepth = 8;
405 <        rp->minweight = 2e-3;
365 <        setcolor(rp->ambval, 0., 0., 0.);
404 >        rp->maxdepth = -10;
405 >        rp->minweight = 1e-4;
406          memset(rp->ambfile, '\0', sizeof(rp->ambfile));
407 +        setcolor(rp->ambval, 0., 0., 0.);
408          rp->ambvwt = 0;
409          rp->ambres = 256;
410 <        rp->ambacc = 0.15;
410 >        rp->ambacc = 0.1;
411          rp->ambdiv = 1024;
412          rp->ambssamp = 512;
413          rp->ambounce = 0;
# Line 374 | Line 415 | ray_defaults(          /* get default parameter values */
415          memset(rp->amblval, '\0', sizeof(rp->amblval));
416          for (i = AMBLLEN+1; i--; )
417                  rp->amblndx[i] = -1;
418 +        
419 +        /* PMAP: restore photon mapping defaults */
420 +        ray_defaults_pmap(rp);
421   }

Diff Legend

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