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.16 by greg, Wed Apr 5 06:22:56 2006 UTC vs.
Revision 2.31 by greg, Wed Aug 21 20:42:20 2024 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  
107   char    *octname;                       /* octree name we are given */
108  
107 – char    *shm_boundary = NULL;           /* boundary of shared memory */
108 –
109   CUBE    thescene;                       /* our scene */
110   OBJECT  nsceneobjs;                     /* number of objects in our scene */
111  
# Line 117 | Line 117 | void   (*trace)() = NULL;              /* trace call */
117  
118   void    (*addobjnotify[8])() = {ambnotify, NULL};
119  
120 + int     castonly = 0;                   /* only doing ray-casting? */
121 +
122   int     do_irrad = 0;                   /* compute irradiance? */
123  
124 < int     rand_samp = 0;                  /* pure Monte Carlo sampling? */
124 > int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
125  
126   double  dstrsrc = 0.0;                  /* square source distribution */
127   double  shadthresh = .03;               /* shadow threshold */
# Line 139 | Line 141 | double specjitter = 1.;                /* specular sampling jitter *
141  
142   int     backvis = 1;                    /* back face visibility */
143  
144 < int     maxdepth = 8;                   /* maximum recursion depth */
145 < double  minweight = 2e-3;               /* minimum ray weight */
144 > int     maxdepth = -10;                 /* maximum recursion depth */
145 > double  minweight = 1e-4;               /* minimum ray weight */
146  
147   char    *ambfile = NULL;                /* ambient file name */
148   COLOR   ambval = BLKCOLOR;              /* ambient value */
# Line 154 | Line 156 | char   *amblist[AMBLLEN+1];            /* ambient include/exclude
156   int     ambincl = -1;                   /* include == 1, exclude == 0 */
157  
158  
159 < extern void
159 > static void
160 > reset_random(void)              /* re-initialize random number generator */
161 > {
162 >        if (rand_samp) {
163 >                srandom((long)time(0));
164 >                initurand(0);
165 >        } else {
166 >                srandom(0L);
167 >                initurand(2048);
168 >        }
169 > }
170 >
171 >
172 > void
173   ray_init(                       /* initialize ray-tracing calculation */
174          char    *otnm
175   )
# Line 165 | Line 180 | ray_init(                      /* initialize ray-tracing calculation */
180          if (ofun[OBJ_SPHERE].funp == o_default)
181                  initotypes();
182                                          /* initialize urand */
183 <        if (rand_samp) {
184 <                srandom((long)time(0));
185 <                initurand(0);
186 <        } else {
187 <                srandom(0L);
188 <                initurand(2048);
189 <        }
175 <                                        /* read scene octree */
176 <        readoct(octname = otnm, ~(IO_FILES|IO_INFO), &thescene, NULL);
183 >        reset_random();
184 >                                        /* initialize spectral sampling */
185 >        if (setspectrsamp(CNDX, WLPART) < 0)
186 >                error(USER, "unsupported spectral sampling");
187 >
188 >        octname = savqstr(otnm);        /* read scene octree */
189 >        readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
190          nsceneobjs = nobjects;
191 <                                        /* find and mark sources */
192 <        marksources();
193 <                                        /* initialize ambient calculation */
194 <        setambient();
195 <                                        /* ready to go... */
191 >
192 >        if (!castonly) {                /* any actual ray traversal to do? */
193 >
194 >                ray_init_pmap();        /* PMAP: set up & load photon maps */
195 >
196 >                marksources();          /* find and mark sources */
197 >
198 >                setambient();           /* initialize ambient calculation */
199 >        } else
200 >                distantsources();       /* else mark only distant sources */
201   }
202  
203 < extern void
203 >
204 > void
205   ray_trace(                      /* trace a primary ray */
206          RAY     *r
207   )
# Line 193 | Line 212 | ray_trace(                     /* trace a primary ray */
212   }
213  
214  
215 < extern void
215 > void
216   ray_done(               /* free ray-tracing data */
217          int     freall
218   )
# Line 206 | Line 225 | ray_done(              /* free ray-tracing data */
225          donesets();
226          octdone();
227          thescene.cutree = EMPTY;
228 <        octname = NULL;
228 >        freeqstr(octname); octname = NULL;
229 >        retainfonts = 0;
230          if (freall) {
211 –                retainfonts = 0;
231                  freefont(NULL);
232                  freedata(NULL);
233 +                SDfreeCache(NULL);
234                  initurand(0);
235          }
236          if (nobjects > 0) {
237                  sprintf(errmsg, "%ld objects left after call to ray_done()",
238 <                                nobjects);
238 >                                (long)nobjects);
239                  error(WARNING, errmsg);
240          }
241 +        
242 +        ray_done_pmap();
243   }
244  
245  
246 < extern void
246 > void
247   ray_save(                       /* save current parameter settings */
248          RAYPARAMS       *rp
249   )
# Line 231 | Line 253 | ray_save(                      /* save current parameter settings */
253          if (rp == NULL)
254                  return;
255          rp->do_irrad = do_irrad;
256 +        rp->rand_samp = rand_samp;
257          rp->dstrsrc = dstrsrc;
258          rp->shadthresh = shadthresh;
259          rp->shadcert = shadcert;
# Line 247 | Line 270 | ray_save(                      /* save current parameter settings */
270          rp->backvis = backvis;
271          rp->maxdepth = maxdepth;
272          rp->minweight = minweight;
250 –        copycolor(rp->ambval, ambval);
251 –        memset(rp->ambfile, '\0', sizeof(rp->ambfile));
273          if (ambfile != NULL)
274                  strncpy(rp->ambfile, ambfile, sizeof(rp->ambfile)-1);
275 +        else
276 +                memset(rp->ambfile, '\0', sizeof(rp->ambfile));
277 +        copycolor(rp->ambval, ambval);
278          rp->ambvwt = ambvwt;
279          rp->ambacc = ambacc;
280          rp->ambres = ambres;
# Line 265 | Line 289 | ray_save(                      /* save current parameter settings */
289                  if (ndx+len >= sizeof(rp->amblval))
290                          break;
291                  strcpy(rp->amblval+ndx, amblist[i]);
292 +                rp->amblndx[i] = ndx;
293                  ndx += len+1;
294          }
295          while (i <= AMBLLEN)
296                  rp->amblndx[i++] = -1;
297 +                
298 +        /* PMAP: save photon mapping params */
299 +        ray_save_pmap(rp);
300   }
301  
302  
303 < extern void
303 > void
304   ray_restore(                    /* restore parameter settings */
305          RAYPARAMS       *rp
306   )
307   {
308 <        register int    i;
308 >        int     i;
309  
310          if (rp == NULL) {               /* restore defaults */
311                  RAYPARAMS       dflt;
# Line 287 | Line 315 | ray_restore(                   /* restore parameter settings */
315          }
316                                          /* restore saved settings */
317          do_irrad = rp->do_irrad;
318 +        if (!rand_samp != !rp->rand_samp) {
319 +                rand_samp = rp->rand_samp;
320 +                reset_random();
321 +        }
322          dstrsrc = rp->dstrsrc;
323          shadthresh = rp->shadthresh;
324          shadcert = rp->shadcert;
# Line 308 | Line 340 | ray_restore(                   /* restore parameter settings */
340          ambdiv = rp->ambdiv;
341          ambssamp = rp->ambssamp;
342          ambounce = rp->ambounce;
343 +                                        /* a bit dangerous if not static */
344          for (i = 0; rp->amblndx[i] >= 0; i++)
345                  amblist[i] = rp->amblval + rp->amblndx[i];
346          while (i <= AMBLLEN)
# Line 315 | Line 348 | ray_restore(                   /* restore parameter settings */
348          ambincl = rp->ambincl;
349                                          /* update ambient calculation */
350          ambnotify(OVOID);
351 <        if (thescene.cutree != EMPTY) {
351 >        if ((thescene.cutree != EMPTY) & !castonly) {
352                  int     newamb = (ambfile == NULL) ?  rp->ambfile[0] :
353                                          strcmp(ambfile, rp->ambfile) ;
354  
# Line 337 | Line 370 | ray_restore(                   /* restore parameter settings */
370                  ambres = rp->ambres;
371                  ambacc = rp->ambacc;
372          }
373 +        
374 +        /* PMAP: restore photon mapping params */
375 +        ray_restore_pmap(rp);
376   }
377  
378  
379 < extern void
379 > void
380   ray_defaults(           /* get default parameter values */
381          RAYPARAMS       *rp
382   )
# Line 351 | Line 387 | ray_defaults(          /* get default parameter values */
387                  return;
388  
389          rp->do_irrad = 0;
390 +        rp->rand_samp = 1;
391          rp->dstrsrc = 0.0;
392 <        rp->shadthresh = .03;
393 <        rp->shadcert = .75;
392 >        rp->shadthresh = 0.03;
393 >        rp->shadcert = 0.75;
394          rp->directrelay = 2;
395          rp->vspretest = 512;
396          rp->directvis = 1;
# Line 362 | Line 399 | ray_defaults(          /* get default parameter values */
399          setcolor(rp->salbedo, 0., 0., 0.);
400          rp->seccg = 0.;
401          rp->ssampdist = 0.;
402 <        rp->specthresh = .15;
402 >        rp->specthresh = 0.15;
403          rp->specjitter = 1.;
404          rp->backvis = 1;
405 <        rp->maxdepth = 8;
406 <        rp->minweight = 2e-3;
370 <        setcolor(rp->ambval, 0., 0., 0.);
405 >        rp->maxdepth = -10;
406 >        rp->minweight = 1e-4;
407          memset(rp->ambfile, '\0', sizeof(rp->ambfile));
408 +        setcolor(rp->ambval, 0., 0., 0.);
409          rp->ambvwt = 0;
410          rp->ambres = 256;
411 <        rp->ambacc = 0.15;
411 >        rp->ambacc = 0.1;
412          rp->ambdiv = 1024;
413          rp->ambssamp = 512;
414          rp->ambounce = 0;
# Line 379 | Line 416 | ray_defaults(          /* get default parameter values */
416          memset(rp->amblval, '\0', sizeof(rp->amblval));
417          for (i = AMBLLEN+1; i--; )
418                  rp->amblndx[i] = -1;
419 +        
420 +        /* PMAP: restore photon mapping defaults */
421 +        ray_defaults_pmap(rp);
422   }

Diff Legend

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