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 = 1; /* pure Monte Carlo sampling? */ |
142 |
|
int backvis = 1; /* back face visibility */ |
143 |
|
|
144 |
|
int maxdepth = -10; /* maximum recursion depth */ |
145 |
< |
double minweight = 2e-3; /* minimum ray weight */ |
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); |
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 |
< |
rand_samp = rp->rand_samp; |
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 |
|
|
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 = -10; |
406 |
< |
rp->minweight = 2e-3; |
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; |