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.14 by greg, Wed Jun 15 15:36:52 2005 UTC vs.
Revision 2.32 by greg, Wed Apr 23 02:35:26 2025 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  "func.h"
102   #include  "data.h"
103   #include  "font.h"
104 + #include  "pmapray.h"
105  
106   char    *progname = "unknown_app";      /* caller sets to argv[0] */
107  
108   char    *octname;                       /* octree name we are given */
109  
107 – char    *shm_boundary = NULL;           /* boundary of shared memory */
108 –
110   CUBE    thescene;                       /* our scene */
111   OBJECT  nsceneobjs;                     /* number of objects in our scene */
112  
# Line 117 | Line 118 | void   (*trace)() = NULL;              /* trace call */
118  
119   void    (*addobjnotify[8])() = {ambnotify, NULL};
120  
121 + int     castonly = 0;                   /* only doing ray-casting? */
122 +
123   int     do_irrad = 0;                   /* compute irradiance? */
124  
125 < int     rand_samp = 0;                  /* pure Monte Carlo sampling? */
125 > int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
126  
127   double  dstrsrc = 0.0;                  /* square source distribution */
128   double  shadthresh = .03;               /* shadow threshold */
# Line 139 | Line 142 | double specjitter = 1.;                /* specular sampling jitter *
142  
143   int     backvis = 1;                    /* back face visibility */
144  
145 < int     maxdepth = 8;                   /* maximum recursion depth */
146 < double  minweight = 2e-3;               /* minimum ray weight */
145 > int     maxdepth = -10;                 /* maximum recursion depth */
146 > double  minweight = 1e-4;               /* minimum ray weight */
147  
148   char    *ambfile = NULL;                /* ambient file name */
149   COLOR   ambval = BLKCOLOR;              /* ambient value */
# Line 154 | Line 157 | char   *amblist[AMBLLEN+1];            /* ambient include/exclude
157   int     ambincl = -1;                   /* include == 1, exclude == 0 */
158  
159  
160 < extern void
160 > static void
161 > reset_random(void)              /* re-initialize random number generator */
162 > {
163 >        if (rand_samp) {
164 >                srandom((long)time(0));
165 >                initurand(0);
166 >        } else {
167 >                srandom(0L);
168 >                initurand(2048);
169 >        }
170 > }
171 >
172 >
173 > void
174   ray_init(                       /* initialize ray-tracing calculation */
175          char    *otnm
176   )
# Line 164 | Line 180 | ray_init(                      /* initialize ray-tracing calculation */
180                                          /* initialize object types */
181          if (ofun[OBJ_SPHERE].funp == o_default)
182                  initotypes();
183 +                                        /* initialize calcomp routines */
184 +        initfunc();
185                                          /* initialize urand */
186 <        initurand(2048);
187 <        srandom(rand_samp ? (long)time(0) : 0L);
188 <                                        /* read scene octree */
189 <        readoct(octname = otnm, ~(IO_FILES|IO_INFO), &thescene, NULL);
186 >        reset_random();
187 >                                        /* initialize spectral sampling */
188 >        if (setspectrsamp(CNDX, WLPART) < 0)
189 >                error(USER, "unsupported spectral sampling");
190 >
191 >        octname = savqstr(otnm);        /* read scene octree */
192 >        readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
193          nsceneobjs = nobjects;
194 <                                        /* find and mark sources */
195 <        marksources();
196 <                                        /* initialize ambient calculation */
197 <        setambient();
198 <                                        /* ready to go... */
194 >
195 >        if (!castonly) {                /* any actual ray traversal to do? */
196 >
197 >                ray_init_pmap();        /* PMAP: set up & load photon maps */
198 >
199 >                marksources();          /* find and mark sources */
200 >
201 >                setambient();           /* initialize ambient calculation */
202 >        } else
203 >                distantsources();       /* else mark only distant sources */
204   }
205  
206 < extern void
206 >
207 > void
208   ray_trace(                      /* trace a primary ray */
209          RAY     *r
210   )
211   {
212          rayorigin(r, PRIMARY, NULL, NULL);
213 <        samplendx = rand_samp ? random() : samplendx+1;
213 >        samplendx++;
214          rayvalue(r);            /* assumes origin and direction are set */
215   }
216  
217  
218 < extern void
218 > void
219   ray_done(               /* free ray-tracing data */
220          int     freall
221   )
# Line 201 | Line 228 | ray_done(              /* free ray-tracing data */
228          donesets();
229          octdone();
230          thescene.cutree = EMPTY;
231 <        octname = NULL;
231 >        freeqstr(octname); octname = NULL;
232 >        retainfonts = 0;
233          if (freall) {
206 –                retainfonts = 0;
234                  freefont(NULL);
235                  freedata(NULL);
236 +                SDfreeCache(NULL);
237                  initurand(0);
238          }
239          if (nobjects > 0) {
240                  sprintf(errmsg, "%ld objects left after call to ray_done()",
241 <                                nobjects);
241 >                                (long)nobjects);
242                  error(WARNING, errmsg);
243          }
244 +        
245 +        ray_done_pmap();
246   }
247  
248  
249 < extern void
249 > void
250   ray_save(                       /* save current parameter settings */
251          RAYPARAMS       *rp
252   )
# Line 226 | Line 256 | ray_save(                      /* save current parameter settings */
256          if (rp == NULL)
257                  return;
258          rp->do_irrad = do_irrad;
259 +        rp->rand_samp = rand_samp;
260          rp->dstrsrc = dstrsrc;
261          rp->shadthresh = shadthresh;
262          rp->shadcert = shadcert;
# Line 242 | Line 273 | ray_save(                      /* save current parameter settings */
273          rp->backvis = backvis;
274          rp->maxdepth = maxdepth;
275          rp->minweight = minweight;
245 –        copycolor(rp->ambval, ambval);
246 –        memset(rp->ambfile, '\0', sizeof(rp->ambfile));
276          if (ambfile != NULL)
277                  strncpy(rp->ambfile, ambfile, sizeof(rp->ambfile)-1);
278 +        else
279 +                memset(rp->ambfile, '\0', sizeof(rp->ambfile));
280 +        copycolor(rp->ambval, ambval);
281          rp->ambvwt = ambvwt;
282          rp->ambacc = ambacc;
283          rp->ambres = ambres;
# Line 260 | Line 292 | ray_save(                      /* save current parameter settings */
292                  if (ndx+len >= sizeof(rp->amblval))
293                          break;
294                  strcpy(rp->amblval+ndx, amblist[i]);
295 +                rp->amblndx[i] = ndx;
296                  ndx += len+1;
297          }
298          while (i <= AMBLLEN)
299                  rp->amblndx[i++] = -1;
300 +                
301 +        /* PMAP: save photon mapping params */
302 +        ray_save_pmap(rp);
303   }
304  
305  
306 < extern void
306 > void
307   ray_restore(                    /* restore parameter settings */
308          RAYPARAMS       *rp
309   )
310   {
311 <        register int    i;
311 >        int     i;
312  
313          if (rp == NULL) {               /* restore defaults */
314                  RAYPARAMS       dflt;
# Line 282 | Line 318 | ray_restore(                   /* restore parameter settings */
318          }
319                                          /* restore saved settings */
320          do_irrad = rp->do_irrad;
321 +        if (!rand_samp != !rp->rand_samp) {
322 +                rand_samp = rp->rand_samp;
323 +                reset_random();
324 +        }
325          dstrsrc = rp->dstrsrc;
326          shadthresh = rp->shadthresh;
327          shadcert = rp->shadcert;
# Line 303 | Line 343 | ray_restore(                   /* restore parameter settings */
343          ambdiv = rp->ambdiv;
344          ambssamp = rp->ambssamp;
345          ambounce = rp->ambounce;
346 +                                        /* a bit dangerous if not static */
347          for (i = 0; rp->amblndx[i] >= 0; i++)
348                  amblist[i] = rp->amblval + rp->amblndx[i];
349          while (i <= AMBLLEN)
# Line 310 | Line 351 | ray_restore(                   /* restore parameter settings */
351          ambincl = rp->ambincl;
352                                          /* update ambient calculation */
353          ambnotify(OVOID);
354 <        if (thescene.cutree != EMPTY) {
354 >        if ((thescene.cutree != EMPTY) & !castonly) {
355                  int     newamb = (ambfile == NULL) ?  rp->ambfile[0] :
356                                          strcmp(ambfile, rp->ambfile) ;
357  
# Line 332 | Line 373 | ray_restore(                   /* restore parameter settings */
373                  ambres = rp->ambres;
374                  ambacc = rp->ambacc;
375          }
376 +        
377 +        /* PMAP: restore photon mapping params */
378 +        ray_restore_pmap(rp);
379   }
380  
381  
382 < extern void
382 > void
383   ray_defaults(           /* get default parameter values */
384          RAYPARAMS       *rp
385   )
# Line 346 | Line 390 | ray_defaults(          /* get default parameter values */
390                  return;
391  
392          rp->do_irrad = 0;
393 +        rp->rand_samp = 1;
394          rp->dstrsrc = 0.0;
395 <        rp->shadthresh = .03;
396 <        rp->shadcert = .75;
395 >        rp->shadthresh = 0.03;
396 >        rp->shadcert = 0.75;
397          rp->directrelay = 2;
398          rp->vspretest = 512;
399          rp->directvis = 1;
# Line 357 | Line 402 | ray_defaults(          /* get default parameter values */
402          setcolor(rp->salbedo, 0., 0., 0.);
403          rp->seccg = 0.;
404          rp->ssampdist = 0.;
405 <        rp->specthresh = .15;
405 >        rp->specthresh = 0.15;
406          rp->specjitter = 1.;
407          rp->backvis = 1;
408 <        rp->maxdepth = 8;
409 <        rp->minweight = 2e-3;
365 <        setcolor(rp->ambval, 0., 0., 0.);
408 >        rp->maxdepth = -10;
409 >        rp->minweight = 1e-4;
410          memset(rp->ambfile, '\0', sizeof(rp->ambfile));
411 +        setcolor(rp->ambval, 0., 0., 0.);
412          rp->ambvwt = 0;
413          rp->ambres = 256;
414 <        rp->ambacc = 0.15;
414 >        rp->ambacc = 0.1;
415          rp->ambdiv = 1024;
416          rp->ambssamp = 512;
417          rp->ambounce = 0;
# Line 374 | Line 419 | ray_defaults(          /* get default parameter values */
419          memset(rp->amblval, '\0', sizeof(rp->amblval));
420          for (i = AMBLLEN+1; i--; )
421                  rp->amblndx[i] = -1;
422 +        
423 +        /* PMAP: restore photon mapping defaults */
424 +        ray_defaults_pmap(rp);
425   }

Diff Legend

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