11 |
|
#include "otypes.h" |
12 |
|
#include "platform.h" |
13 |
|
|
14 |
+ |
char *shm_boundary = NULL; /* boundary of shared memory */ |
15 |
+ |
|
16 |
|
CUBE thescene; /* our scene */ |
17 |
|
OBJECT nsceneobjs; /* number of objects in our scene */ |
18 |
|
|
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 |
|
|
67 |
– |
/* Close output stream and free record */ |
68 |
– |
static void |
69 |
– |
closestream(void *p) |
70 |
– |
{ |
71 |
– |
STREAMOUT *sop = (STREAMOUT *)p; |
72 |
– |
int status = 0; |
73 |
– |
if (sop->outpipe) |
74 |
– |
status = pclose(sop->ofp); |
75 |
– |
else if (sop->ofp != stdout) |
76 |
– |
status = fclose(sop->ofp); |
77 |
– |
if (status) |
78 |
– |
error(SYSTEM, "error closing output stream"); |
79 |
– |
free(p); |
80 |
– |
} |
81 |
– |
|
82 |
– |
LUTAB ofiletab = LU_SINIT(free,closestream); /* output file table */ |
83 |
– |
|
84 |
– |
#define OF_MODIFIER 01 |
85 |
– |
#define OF_BIN 02 |
86 |
– |
|
69 |
|
static void mcfree(void *p) { epfree((*(MODCONT *)p).binv); free(p); } |
70 |
|
|
71 |
|
LUTAB modconttab = LU_SINIT(NULL,mcfree); /* modifier lookup table */ |
204 |
|
if (nproc > MAXPROCESS) |
205 |
|
sprintf(errmsg, "too many processes requested -- reducing to %d", |
206 |
|
nproc = MAXPROCESS); |
207 |
< |
|
207 |
> |
if (nproc > 1) { |
208 |
> |
preload_objs(); /* preload auxiliary data */ |
209 |
> |
/* set shared memory boundary */ |
210 |
> |
shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY"); |
211 |
> |
} |
212 |
|
if ((nproc > 1) & (accumulate <= 0)) |
213 |
< |
zero_record(0); /* prime our queue to accumulate */ |
213 |
> |
put_zero_record(0); /* prime our queue to accumulate */ |
214 |
|
|
229 |
– |
if (recover) { /* recover previous output? */ |
230 |
– |
if (accumulate <= 0) { |
231 |
– |
reload_output(); |
232 |
– |
if ((nproc > 1) & (accumulate <= 0)) |
233 |
– |
queue_modifiers(); |
234 |
– |
} else |
235 |
– |
recover_output(); |
236 |
– |
} |
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 |
|
} |
345 |
|
static int ignore_warning_given = 0; |
346 |
|
FVECT orig, direc; |
347 |
|
double d; |
348 |
< |
/* initialize & fork */ |
348 |
> |
/* initialize (& fork more of us) */ |
349 |
|
rcinit(); |
350 |
|
/* load rays from stdin & process */ |
351 |
|
#ifdef getc_unlocked |
352 |
< |
flockfile(stdin); /* avoid lock/unlock overhead */ |
352 |
> |
flockfile(stdin); /* avoid mutex overhead */ |
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; |