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 = 1; /* pure Monte Carlo sampling? */ |
144 |
|
int backvis = 1; /* back face visibility */ |
145 |
|
|
146 |
|
int maxdepth = -10; /* maximum recursion depth */ |
147 |
< |
double minweight = 2e-3; /* minimum ray weight */ |
147 |
> |
double minweight = 1e-4; /* minimum ray weight */ |
148 |
|
|
149 |
|
char *ambfile = NULL; /* ambient file name */ |
150 |
|
COLOR ambval = BLKCOLOR; /* ambient value */ |
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 |
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 { |
174 |
< |
srandom(0L); |
175 |
< |
initurand(2048); |
176 |
< |
} |
177 |
< |
/* read scene octree */ |
178 |
< |
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 |
< |
/* PMAP: Init & load photon maps */ |
191 |
< |
ray_init_pmap(); |
192 |
< |
/* find and mark sources */ |
193 |
< |
marksources(); |
194 |
< |
/* initialize ambient calculation */ |
195 |
< |
setambient(); |
196 |
< |
/* ready to go... (almost) */ |
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 |
+ |
|
203 |
|
void |
204 |
|
ray_trace( /* trace a primary ray */ |
205 |
|
RAY *r |
304 |
|
RAYPARAMS *rp |
305 |
|
) |
306 |
|
{ |
307 |
< |
register int i; |
307 |
> |
int i; |
308 |
|
|
309 |
|
if (rp == NULL) { /* restore defaults */ |
310 |
|
RAYPARAMS dflt; |
314 |
|
} |
315 |
|
/* restore saved settings */ |
316 |
|
do_irrad = rp->do_irrad; |
317 |
< |
rand_samp = rp->rand_samp; |
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; |
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) |
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; |
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 = -10; |
405 |
< |
rp->minweight = 2e-3; |
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; |