106 |
|
|
107 |
|
char *octname; /* octree name we are given */ |
108 |
|
|
109 |
– |
char *shm_boundary = NULL; /* boundary of shared memory */ |
110 |
– |
|
109 |
|
CUBE thescene; /* our scene */ |
110 |
|
OBJECT nsceneobjs; /* number of objects in our scene */ |
111 |
|
|
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 */ |
141 |
|
|
142 |
|
int backvis = 1; /* back face visibility */ |
143 |
|
|
144 |
< |
int maxdepth = 8; /* maximum recursion depth */ |
145 |
< |
double minweight = 5e-4; /* 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 */ |
156 |
|
int ambincl = -1; /* include == 1, exclude == 0 */ |
157 |
|
|
158 |
|
|
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 |
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 |
< |
} |
177 |
< |
/* read scene octree */ |
178 |
< |
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 |
< |
/* PMAP: Init & load photon maps */ |
192 |
< |
ray_init_pmap(); |
193 |
< |
/* find and mark sources */ |
194 |
< |
marksources(); |
195 |
< |
/* initialize ambient calculation */ |
196 |
< |
setambient(); |
197 |
< |
/* ready to go... (almost) */ |
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 |
+ |
|
204 |
|
void |
205 |
|
ray_trace( /* trace a primary ray */ |
206 |
|
RAY *r |
225 |
|
donesets(); |
226 |
|
octdone(); |
227 |
|
thescene.cutree = EMPTY; |
228 |
< |
octname = NULL; |
228 |
> |
freeqstr(octname); octname = NULL; |
229 |
|
retainfonts = 0; |
230 |
|
if (freall) { |
231 |
|
freefont(NULL); |
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; |
270 |
|
rp->backvis = backvis; |
271 |
|
rp->maxdepth = maxdepth; |
272 |
|
rp->minweight = minweight; |
257 |
– |
copycolor(rp->ambval, ambval); |
258 |
– |
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; |
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) |
305 |
|
RAYPARAMS *rp |
306 |
|
) |
307 |
|
{ |
308 |
< |
register int i; |
308 |
> |
int i; |
309 |
|
|
310 |
|
if (rp == NULL) { /* restore defaults */ |
311 |
|
RAYPARAMS dflt; |
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; |
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) |
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 |
|
|
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; |
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; |
383 |
< |
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; |