--- ray/src/rt/rcontrib.c 2012/06/11 05:07:55 2.3 +++ ray/src/rt/rcontrib.c 2012/06/13 00:16:42 2.5 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: rcontrib.c,v 2.3 2012/06/11 05:07:55 greg Exp $"; +static const char RCSid[] = "$Id: rcontrib.c,v 2.5 2012/06/13 00:16:42 greg Exp $"; #endif /* * Accumulate ray contributions for a set of materials @@ -63,8 +63,8 @@ int account; /* current accumulation count */ RNUMBER raysleft; /* number of rays left to trace */ long waitflush; /* how long until next flush */ -int lastray = 0; /* last ray number sent */ -int lastdone = 0; /* last ray output */ +RNUMBER lastray = 0; /* last ray number sent */ +RNUMBER lastdone = 0; /* last ray output */ static void mcfree(void *p) { epfree((*(MODCONT *)p).binv); free(p); } @@ -212,14 +212,6 @@ rcinit() if ((nproc > 1) & (accumulate <= 0)) put_zero_record(0); /* prime our queue to accumulate */ - if (recover) { /* recover previous output? */ - if (accumulate <= 0) { - reload_output(); - if (nproc > 1) - queue_modifiers(); - } else - recover_output(); - } if (yres > 0) { /* set up flushing & ray counts */ if (xres > 0) raysleft = (RNUMBER)xres*yres; @@ -234,9 +226,20 @@ rcinit() for (i = 0; i < nsources; i++) source[i].sflags |= SFOLLOW; - if (nproc == 1 || in_rchild()) /* single process or child */ + if (nproc > 1 && in_rchild()) /* forked child? */ return; /* return to main processing loop */ + if (recover) { /* recover previous output? */ + if (accumulate <= 0) { + reload_output(); + if (nproc > 1) + queue_modifiers(); + } else + recover_output(); + } + if (nproc == 1) /* single process? */ + return; + parental_loop(); /* else run controller */ quit(0); /* parent musn't return! */ } @@ -350,7 +353,7 @@ rcontrib() #endif while (getvec(orig) == 0 && getvec(direc) == 0) { d = normalize(direc); - if ((d == 0.0) & (accumulate != 1)) { + if (nchild != -1 && (d == 0.0) & (accumulate != 1)) { if (!ignore_warning_given++) error(WARNING, "dummy ray(s) ignored during accumulation\n"); @@ -362,6 +365,7 @@ rcontrib() if (d == 0.0) { /* zero ==> flush */ if ((yres <= 0) | (xres <= 0)) waitflush = 1; /* flush right after */ + account = 1; } else { /* else compute */ eval_ray(orig, direc, lim_dist ? d : 0.0); } @@ -370,7 +374,7 @@ rcontrib() if (raysleft && !--raysleft) break; /* preemptive EOI */ } - if (accumulate <= 0 || account < accumulate) { + if (nchild != -1 && (accumulate <= 0) | (account < accumulate)) { if (account < accumulate) { error(WARNING, "partial accumulation in final record"); accumulate -= account;