| 59 |
|
|
| 60 |
|
typedef void oputf_t(RAY *r); |
| 61 |
|
static oputf_t oputo, oputd, oputv, oputV, oputl, oputL, oputc, oputp, |
| 62 |
< |
oputn, oputN, oputs, oputw, oputW, oputm, oputM, oputtilde; |
| 62 |
> |
oputr, oputR, oputx, oputX, oputn, oputN, oputs, |
| 63 |
> |
oputw, oputW, oputm, oputM, oputtilde; |
| 64 |
|
|
| 65 |
|
static void setoutput(char *vs); |
| 66 |
|
extern void tranotify(OBJECT obj); |
| 73 |
|
static void tabin(RAY *r); |
| 74 |
|
static void ourtrace(RAY *r); |
| 75 |
|
|
| 76 |
< |
static oputf_t *ray_out[16], *every_out[16]; |
| 76 |
> |
static oputf_t *ray_out[32], *every_out[32]; |
| 77 |
|
static putf_t *putreal; |
| 78 |
|
|
| 79 |
|
|
| 118 |
|
unsigned long vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu |
| 119 |
|
: (unsigned long)vresolu; |
| 120 |
|
long nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : hresolu; |
| 121 |
+ |
int something2flush = 0; |
| 122 |
|
FILE *fp; |
| 123 |
|
double d; |
| 124 |
|
FVECT orig, direc; |
| 137 |
|
castonly = 0; |
| 138 |
|
else if (castonly) |
| 139 |
|
nproc = 1; /* don't bother multiprocessing */ |
| 140 |
+ |
if ((nextflush > 0) & (nproc > nextflush)) { |
| 141 |
+ |
error(WARNING, "reducing number of processes to match flush interval"); |
| 142 |
+ |
nproc = nextflush; |
| 143 |
+ |
} |
| 144 |
|
switch (outform) { |
| 145 |
|
case 'a': putreal = puta; break; |
| 146 |
|
case 'f': putreal = putf; break; |
| 159 |
|
if (hresolu > 0) { |
| 160 |
|
if (vresolu > 0) |
| 161 |
|
fprtresolu(hresolu, vresolu, stdout); |
| 162 |
< |
fflush(stdout); |
| 162 |
> |
else |
| 163 |
> |
fflush(stdout); |
| 164 |
|
} |
| 165 |
|
/* process file */ |
| 166 |
|
while (getvec(orig, inform, fp) == 0 && |
| 168 |
|
|
| 169 |
|
d = normalize(direc); |
| 170 |
|
if (d == 0.0) { /* zero ==> flush */ |
| 171 |
< |
if (--nextflush <= 0 || !vcount) { |
| 172 |
< |
if (nproc > 1 && ray_fifo_flush() < 0) |
| 171 |
> |
if ((--nextflush <= 0) | !vcount && something2flush) { |
| 172 |
> |
if (ray_pnprocs > 1 && ray_fifo_flush() < 0) |
| 173 |
|
error(USER, "child(ren) died"); |
| 174 |
|
bogusray(); |
| 175 |
|
fflush(stdout); |
| 176 |
+ |
something2flush = 0; |
| 177 |
|
nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : |
| 178 |
|
hresolu; |
| 179 |
|
} else |
| 182 |
|
rtcompute(orig, direc, lim_dist ? d : 0.0); |
| 183 |
|
/* flush if time */ |
| 184 |
|
if (!--nextflush) { |
| 185 |
< |
if (nproc > 1 && ray_fifo_flush() < 0) |
| 185 |
> |
if (ray_pnprocs > 1 && ray_fifo_flush() < 0) |
| 186 |
|
error(USER, "child(ren) died"); |
| 187 |
|
fflush(stdout); |
| 188 |
|
nextflush = hresolu; |
| 189 |
< |
} |
| 189 |
> |
} else |
| 190 |
> |
something2flush = 1; |
| 191 |
|
} |
| 192 |
|
if (ferror(stdout)) |
| 193 |
|
error(SYSTEM, "write error"); |
| 194 |
|
if (vcount && !--vcount) /* check for end */ |
| 195 |
|
break; |
| 196 |
|
} |
| 197 |
< |
if (nproc > 1) { /* clean up children */ |
| 197 |
> |
if (ray_pnprocs > 1) { /* clean up children */ |
| 198 |
|
if (ray_fifo_flush() < 0) |
| 199 |
|
error(USER, "unable to complete processing"); |
| 200 |
|
ray_pclose(0); |
| 245 |
|
case 'd': /* direction */ |
| 246 |
|
*table++ = oputd; |
| 247 |
|
break; |
| 248 |
+ |
case 'r': /* reflected contrib. */ |
| 249 |
+ |
*table++ = oputr; |
| 250 |
+ |
castonly = 0; |
| 251 |
+ |
break; |
| 252 |
+ |
case 'R': /* reflected distance */ |
| 253 |
+ |
*table++ = oputR; |
| 254 |
+ |
castonly = 0; |
| 255 |
+ |
break; |
| 256 |
+ |
case 'x': /* xmit contrib. */ |
| 257 |
+ |
*table++ = oputx; |
| 258 |
+ |
castonly = 0; |
| 259 |
+ |
break; |
| 260 |
+ |
case 'X': /* xmit distance */ |
| 261 |
+ |
*table++ = oputX; |
| 262 |
+ |
castonly = 0; |
| 263 |
+ |
break; |
| 264 |
|
case 'v': /* value */ |
| 265 |
|
*table++ = oputv; |
| 266 |
|
castonly = 0; |
| 267 |
|
break; |
| 268 |
|
case 'V': /* contribution */ |
| 269 |
|
*table++ = oputV; |
| 270 |
+ |
castonly = 0; |
| 271 |
|
if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0)) |
| 272 |
|
error(WARNING, |
| 273 |
|
"-otV accuracy depends on -aa 0 -as 0"); |
| 300 |
|
break; |
| 301 |
|
case 'W': /* coefficient */ |
| 302 |
|
*table++ = oputW; |
| 303 |
+ |
castonly = 0; |
| 304 |
|
if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0)) |
| 305 |
|
error(WARNING, |
| 306 |
|
"-otW accuracy depends on -aa 0 -as 0"); |
| 322 |
|
static void |
| 323 |
|
bogusray(void) /* print out empty record */ |
| 324 |
|
{ |
| 298 |
– |
thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] = |
| 299 |
– |
thisray.rdir[0] = thisray.rdir[1] = thisray.rdir[2] = 0.0; |
| 300 |
– |
thisray.rmax = 0.0; |
| 325 |
|
rayorigin(&thisray, PRIMARY, NULL, NULL); |
| 326 |
|
printvals(&thisray); |
| 327 |
|
} |
| 348 |
|
) |
| 349 |
|
{ |
| 350 |
|
void (*old_revf)(RAY *) = r->revf; |
| 351 |
< |
|
| 352 |
< |
r->rot = 1e-5; /* pretend we hit surface */ |
| 351 |
> |
/* pretend we hit surface */ |
| 352 |
> |
r->rxt = r->rot = 1e-5; |
| 353 |
|
VSUM(r->rop, r->rorg, r->rdir, r->rot); |
| 354 |
|
r->ron[0] = -r->rdir[0]; |
| 355 |
|
r->ron[1] = -r->rdir[1]; |
| 533 |
|
|
| 534 |
|
|
| 535 |
|
static void |
| 536 |
+ |
oputr( /* print mirrored contribution */ |
| 537 |
+ |
RAY *r |
| 538 |
+ |
) |
| 539 |
+ |
{ |
| 540 |
+ |
RREAL cval[3]; |
| 541 |
+ |
|
| 542 |
+ |
cval[0] = colval(r->mcol,RED); |
| 543 |
+ |
cval[1] = colval(r->mcol,GRN); |
| 544 |
+ |
cval[2] = colval(r->mcol,BLU); |
| 545 |
+ |
(*putreal)(cval, 3); |
| 546 |
+ |
} |
| 547 |
+ |
|
| 548 |
+ |
|
| 549 |
+ |
|
| 550 |
+ |
static void |
| 551 |
+ |
oputR( /* print mirrored distance */ |
| 552 |
+ |
RAY *r |
| 553 |
+ |
) |
| 554 |
+ |
{ |
| 555 |
+ |
(*putreal)(&r->rmt, 1); |
| 556 |
+ |
} |
| 557 |
+ |
|
| 558 |
+ |
|
| 559 |
+ |
static void |
| 560 |
+ |
oputx( /* print unmirrored contribution */ |
| 561 |
+ |
RAY *r |
| 562 |
+ |
) |
| 563 |
+ |
{ |
| 564 |
+ |
RREAL cval[3]; |
| 565 |
+ |
|
| 566 |
+ |
cval[0] = colval(r->rcol,RED) - colval(r->mcol,RED); |
| 567 |
+ |
cval[1] = colval(r->rcol,GRN) - colval(r->mcol,GRN); |
| 568 |
+ |
cval[2] = colval(r->rcol,BLU) - colval(r->mcol,BLU); |
| 569 |
+ |
(*putreal)(cval, 3); |
| 570 |
+ |
} |
| 571 |
+ |
|
| 572 |
+ |
|
| 573 |
+ |
static void |
| 574 |
+ |
oputX( /* print unmirrored distance */ |
| 575 |
+ |
RAY *r |
| 576 |
+ |
) |
| 577 |
+ |
{ |
| 578 |
+ |
(*putreal)(&r->rxt, 1); |
| 579 |
+ |
} |
| 580 |
+ |
|
| 581 |
+ |
|
| 582 |
+ |
static void |
| 583 |
|
oputv( /* print value */ |
| 584 |
|
RAY *r |
| 585 |
|
) |
| 619 |
|
RAY *r |
| 620 |
|
) |
| 621 |
|
{ |
| 622 |
< |
(*putreal)(&r->rt, 1); |
| 622 |
> |
RREAL d = raydistance(r); |
| 623 |
> |
|
| 624 |
> |
(*putreal)(&d, 1); |
| 625 |
|
} |
| 626 |
|
|
| 627 |
|
|