| 25 |
|
|
| 26 |
|
int do_irrad = 0; /* compute irradiance? */ |
| 27 |
|
|
| 28 |
< |
int rand_samp = 0; /* pure Monte Carlo sampling? */ |
| 28 |
> |
int rand_samp = 1; /* pure Monte Carlo sampling? */ |
| 29 |
|
|
| 30 |
|
double dstrsrc = 0.0; /* square source distribution */ |
| 31 |
|
double shadthresh = .03; /* shadow threshold */ |
| 45 |
|
|
| 46 |
|
int backvis = 1; /* back face visibility */ |
| 47 |
|
|
| 48 |
< |
int maxdepth = 8; /* maximum recursion depth */ |
| 49 |
< |
double minweight = 5e-4; /* minimum ray weight */ |
| 48 |
> |
int maxdepth = -10; /* maximum recursion depth */ |
| 49 |
> |
double minweight = 2e-3; /* minimum ray weight */ |
| 50 |
|
|
| 51 |
|
char *ambfile = NULL; /* ambient file name */ |
| 52 |
|
COLOR ambval = BLKCOLOR; /* ambient value */ |
| 63 |
|
RNUMBER raysleft; /* number of rays left to trace */ |
| 64 |
|
long waitflush; /* how long until next flush */ |
| 65 |
|
|
| 66 |
< |
int lastray = 0; /* last ray number sent */ |
| 67 |
< |
int lastdone = 0; /* last ray output */ |
| 66 |
> |
RNUMBER lastray = 0; /* last ray number sent */ |
| 67 |
> |
RNUMBER lastdone = 0; /* last ray output */ |
| 68 |
|
|
| 69 |
– |
/* Close output stream and free record */ |
| 70 |
– |
static void |
| 71 |
– |
closestream(void *p) |
| 72 |
– |
{ |
| 73 |
– |
STREAMOUT *sop = (STREAMOUT *)p; |
| 74 |
– |
int status = 0; |
| 75 |
– |
if (sop->outpipe) |
| 76 |
– |
status = pclose(sop->ofp); |
| 77 |
– |
else if (sop->ofp != stdout) |
| 78 |
– |
status = fclose(sop->ofp); |
| 79 |
– |
if (status) |
| 80 |
– |
error(SYSTEM, "error closing output stream"); |
| 81 |
– |
free(p); |
| 82 |
– |
} |
| 83 |
– |
|
| 84 |
– |
LUTAB ofiletab = LU_SINIT(free,closestream); /* output file table */ |
| 85 |
– |
|
| 86 |
– |
#define OF_MODIFIER 01 |
| 87 |
– |
#define OF_BIN 02 |
| 88 |
– |
|
| 69 |
|
static void mcfree(void *p) { epfree((*(MODCONT *)p).binv); free(p); } |
| 70 |
|
|
| 71 |
|
LUTAB modconttab = LU_SINIT(NULL,mcfree); /* modifier lookup table */ |
| 212 |
|
if ((nproc > 1) & (accumulate <= 0)) |
| 213 |
|
put_zero_record(0); /* prime our queue to accumulate */ |
| 214 |
|
|
| 235 |
– |
if (recover) { /* recover previous output? */ |
| 236 |
– |
if (accumulate <= 0) { |
| 237 |
– |
reload_output(); |
| 238 |
– |
if (nproc > 1) |
| 239 |
– |
queue_modifiers(); |
| 240 |
– |
} else |
| 241 |
– |
recover_output(); |
| 242 |
– |
} |
| 215 |
|
if (yres > 0) { /* set up flushing & ray counts */ |
| 216 |
|
if (xres > 0) |
| 217 |
|
raysleft = (RNUMBER)xres*yres; |
| 226 |
|
for (i = 0; i < nsources; i++) |
| 227 |
|
source[i].sflags |= SFOLLOW; |
| 228 |
|
|
| 229 |
< |
if (nproc == 1 || in_rchild()) /* single process or child */ |
| 229 |
> |
if (nproc > 1 && in_rchild()) /* forked child? */ |
| 230 |
|
return; /* return to main processing loop */ |
| 231 |
|
|
| 232 |
+ |
if (recover) { /* recover previous output? */ |
| 233 |
+ |
if (accumulate <= 0) { |
| 234 |
+ |
reload_output(); |
| 235 |
+ |
if (nproc > 1) |
| 236 |
+ |
queue_modifiers(); |
| 237 |
+ |
} else |
| 238 |
+ |
recover_output(); |
| 239 |
+ |
} |
| 240 |
+ |
if (nproc == 1) /* single process? */ |
| 241 |
+ |
return; |
| 242 |
+ |
|
| 243 |
|
parental_loop(); /* else run controller */ |
| 244 |
|
quit(0); /* parent musn't return! */ |
| 245 |
|
} |
| 353 |
|
#endif |
| 354 |
|
while (getvec(orig) == 0 && getvec(direc) == 0) { |
| 355 |
|
d = normalize(direc); |
| 356 |
< |
if ((d == 0.0) & (accumulate != 1)) { |
| 356 |
> |
if (nchild != -1 && (d == 0.0) & (accumulate != 1)) { |
| 357 |
|
if (!ignore_warning_given++) |
| 358 |
|
error(WARNING, |
| 359 |
|
"dummy ray(s) ignored during accumulation\n"); |
| 365 |
|
if (d == 0.0) { /* zero ==> flush */ |
| 366 |
|
if ((yres <= 0) | (xres <= 0)) |
| 367 |
|
waitflush = 1; /* flush right after */ |
| 368 |
+ |
account = 1; |
| 369 |
|
} else { /* else compute */ |
| 370 |
|
eval_ray(orig, direc, lim_dist ? d : 0.0); |
| 371 |
|
} |
| 374 |
|
if (raysleft && !--raysleft) |
| 375 |
|
break; /* preemptive EOI */ |
| 376 |
|
} |
| 377 |
< |
if (accumulate <= 0 || account < accumulate) { |
| 377 |
> |
if (nchild != -1 && (accumulate <= 0) | (account < accumulate)) { |
| 378 |
|
if (account < accumulate) { |
| 379 |
|
error(WARNING, "partial accumulation in final record"); |
| 380 |
|
accumulate -= account; |