28 |
|
#include "ambient.h" |
29 |
|
#include "source.h" |
30 |
|
#include "otypes.h" |
31 |
+ |
#include "otspecial.h" |
32 |
|
#include "resolu.h" |
33 |
|
#include "random.h" |
34 |
|
|
54 |
|
|
55 |
|
static RAY thisray; /* for our convenience */ |
56 |
|
|
57 |
< |
typedef void putf_t(double v); |
57 |
> |
typedef void putf_t(RREAL *v, int n); |
58 |
|
static putf_t puta, putd, putf; |
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 |
|
|
116 |
|
) |
117 |
|
{ |
118 |
|
unsigned long vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu |
119 |
< |
: vresolu; |
120 |
< |
long nextflush = hresolu; |
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 (nproc > 1 && ray_fifo_flush() < 0) |
172 |
< |
error(USER, "lost children"); |
173 |
< |
bogusray(); |
174 |
< |
if (--nextflush <= 0 || !vcount) { |
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 |
< |
nextflush = hresolu; |
177 |
< |
} |
176 |
> |
something2flush = 0; |
177 |
> |
nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : |
178 |
> |
hresolu; |
179 |
> |
} else |
180 |
> |
bogusray(); |
181 |
|
} else { /* compute and print */ |
182 |
|
rtcompute(orig, direc, lim_dist ? d : 0.0); |
183 |
|
/* flush if time */ |
184 |
|
if (!--nextflush) { |
185 |
< |
if (nproc > 1 && ray_fifo_flush() < 0) |
186 |
< |
error(USER, "lost children"); |
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 |
|
{ |
295 |
– |
thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] = |
296 |
– |
thisray.rdir[0] = thisray.rdir[1] = thisray.rdir[2] = 0.0; |
297 |
– |
thisray.rmax = 0.0; |
325 |
|
rayorigin(&thisray, PRIMARY, NULL, NULL); |
326 |
|
printvals(&thisray); |
327 |
|
} |
347 |
|
RAY *r |
348 |
|
) |
349 |
|
{ |
350 |
< |
r->rot = 1e-5; /* pretend we hit surface */ |
350 |
> |
void (*old_revf)(RAY *) = r->revf; |
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]; |
356 |
|
r->ron[2] = -r->rdir[2]; |
357 |
|
r->rod = 1.0; |
358 |
|
/* compute result */ |
359 |
+ |
r->revf = raytrace; |
360 |
|
(*ofun[Lamb.otype].funp)(&Lamb, r); |
361 |
+ |
r->revf = old_revf; |
362 |
|
} |
363 |
|
|
364 |
|
|
369 |
|
double dmax |
370 |
|
) |
371 |
|
{ |
372 |
< |
if (imm_irrad) { /* set up ray */ |
372 |
> |
/* set up ray */ |
373 |
> |
rayorigin(&thisray, PRIMARY, NULL, NULL); |
374 |
> |
if (imm_irrad) { |
375 |
|
VSUM(thisray.rorg, org, dir, 1.1e-4); |
376 |
|
thisray.rdir[0] = -dir[0]; |
377 |
|
thisray.rdir[1] = -dir[1]; |
378 |
|
thisray.rdir[2] = -dir[2]; |
379 |
|
thisray.rmax = 0.0; |
380 |
+ |
thisray.revf = rayirrad; |
381 |
|
} else { |
382 |
|
VCOPY(thisray.rorg, org); |
383 |
|
VCOPY(thisray.rdir, dir); |
384 |
|
thisray.rmax = dmax; |
385 |
+ |
if (castonly) |
386 |
+ |
thisray.revf = raycast; |
387 |
|
} |
352 |
– |
rayorigin(&thisray, PRIMARY, NULL, NULL); |
353 |
– |
/* special case evaluators */ |
354 |
– |
if (castonly) |
355 |
– |
thisray.revf = raycast; |
356 |
– |
else if (imm_irrad) |
357 |
– |
thisray.revf = rayirrad; |
358 |
– |
|
388 |
|
if (ray_pnprocs > 1) { /* multiprocessing FIFO? */ |
389 |
|
if (ray_fifo_in(&thisray) < 0) |
390 |
|
error(USER, "lost children"); |
435 |
|
} |
436 |
|
break; |
437 |
|
case 'f': /* binary float */ |
438 |
< |
if (fread((char *)vf, sizeof(float), 3, fp) != 3) |
438 |
> |
if (getbinary(vf, sizeof(float), 3, fp) != 3) |
439 |
|
return(-1); |
440 |
|
VCOPY(vec, vf); |
441 |
|
break; |
442 |
|
case 'd': /* binary double */ |
443 |
< |
if (fread((char *)vd, sizeof(double), 3, fp) != 3) |
443 |
> |
if (getbinary(vd, sizeof(double), 3, fp) != 3) |
444 |
|
return(-1); |
445 |
|
VCOPY(vec, vd); |
446 |
|
break; |
519 |
|
RAY *r |
520 |
|
) |
521 |
|
{ |
522 |
< |
(*putreal)(r->rorg[0]); |
494 |
< |
(*putreal)(r->rorg[1]); |
495 |
< |
(*putreal)(r->rorg[2]); |
522 |
> |
(*putreal)(r->rorg, 3); |
523 |
|
} |
524 |
|
|
525 |
|
|
528 |
|
RAY *r |
529 |
|
) |
530 |
|
{ |
531 |
< |
(*putreal)(r->rdir[0]); |
505 |
< |
(*putreal)(r->rdir[1]); |
506 |
< |
(*putreal)(r->rdir[2]); |
531 |
> |
(*putreal)(r->rdir, 3); |
532 |
|
} |
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 |
|
) |
586 |
|
{ |
587 |
+ |
RREAL cval[3]; |
588 |
+ |
|
589 |
|
if (outform == 'c') { |
590 |
|
COLR cout; |
591 |
|
setcolr(cout, colval(r->rcol,RED), |
592 |
|
colval(r->rcol,GRN), |
593 |
|
colval(r->rcol,BLU)); |
594 |
< |
fwrite((char *)cout, sizeof(cout), 1, stdout); |
594 |
> |
putbinary(cout, sizeof(cout), 1, stdout); |
595 |
|
return; |
596 |
|
} |
597 |
< |
(*putreal)(colval(r->rcol,RED)); |
598 |
< |
(*putreal)(colval(r->rcol,GRN)); |
599 |
< |
(*putreal)(colval(r->rcol,BLU)); |
597 |
> |
cval[0] = colval(r->rcol,RED); |
598 |
> |
cval[1] = colval(r->rcol,GRN); |
599 |
> |
cval[2] = colval(r->rcol,BLU); |
600 |
> |
(*putreal)(cval, 3); |
601 |
|
} |
602 |
|
|
603 |
|
|
606 |
|
RAY *r |
607 |
|
) |
608 |
|
{ |
609 |
< |
double contr[3]; |
609 |
> |
RREAL contr[3]; |
610 |
|
|
611 |
|
raycontrib(contr, r, PRIMARY); |
612 |
|
multcolor(contr, r->rcol); |
613 |
< |
(*putreal)(contr[RED]); |
539 |
< |
(*putreal)(contr[GRN]); |
540 |
< |
(*putreal)(contr[BLU]); |
613 |
> |
(*putreal)(contr, 3); |
614 |
|
} |
615 |
|
|
616 |
|
|
619 |
|
RAY *r |
620 |
|
) |
621 |
|
{ |
622 |
< |
(*putreal)(r->rt); |
622 |
> |
RREAL d = raydistance(r); |
623 |
> |
|
624 |
> |
(*putreal)(&d, 1); |
625 |
|
} |
626 |
|
|
627 |
|
|
630 |
|
RAY *r |
631 |
|
) |
632 |
|
{ |
633 |
< |
(*putreal)(r->rot); |
633 |
> |
(*putreal)(&r->rot, 1); |
634 |
|
} |
635 |
|
|
636 |
|
|
639 |
|
RAY *r |
640 |
|
) |
641 |
|
{ |
642 |
< |
(*putreal)(r->uv[0]); |
568 |
< |
(*putreal)(r->uv[1]); |
642 |
> |
(*putreal)(r->uv, 2); |
643 |
|
} |
644 |
|
|
645 |
|
|
646 |
+ |
static RREAL vdummy[3] = {0.0, 0.0, 0.0}; |
647 |
+ |
|
648 |
+ |
|
649 |
|
static void |
650 |
|
oputp( /* print point */ |
651 |
|
RAY *r |
652 |
|
) |
653 |
|
{ |
654 |
< |
if (r->rot < FHUGE) { |
655 |
< |
(*putreal)(r->rop[0]); |
656 |
< |
(*putreal)(r->rop[1]); |
657 |
< |
(*putreal)(r->rop[2]); |
581 |
< |
} else { |
582 |
< |
(*putreal)(0.0); |
583 |
< |
(*putreal)(0.0); |
584 |
< |
(*putreal)(0.0); |
585 |
< |
} |
654 |
> |
if (r->rot < FHUGE) |
655 |
> |
(*putreal)(r->rop, 3); |
656 |
> |
else |
657 |
> |
(*putreal)(vdummy, 3); |
658 |
|
} |
659 |
|
|
660 |
|
|
663 |
|
RAY *r |
664 |
|
) |
665 |
|
{ |
666 |
< |
if (r->rot < FHUGE) { |
667 |
< |
(*putreal)(r->ron[0]); |
668 |
< |
(*putreal)(r->ron[1]); |
669 |
< |
(*putreal)(r->ron[2]); |
598 |
< |
} else { |
599 |
< |
(*putreal)(0.0); |
600 |
< |
(*putreal)(0.0); |
601 |
< |
(*putreal)(0.0); |
602 |
< |
} |
666 |
> |
if (r->rot < FHUGE) |
667 |
> |
(*putreal)(r->ron, 3); |
668 |
> |
else |
669 |
> |
(*putreal)(vdummy, 3); |
670 |
|
} |
671 |
|
|
672 |
|
|
678 |
|
FVECT pnorm; |
679 |
|
|
680 |
|
if (r->rot >= FHUGE) { |
681 |
< |
(*putreal)(0.0); |
615 |
< |
(*putreal)(0.0); |
616 |
< |
(*putreal)(0.0); |
681 |
> |
(*putreal)(vdummy, 3); |
682 |
|
return; |
683 |
|
} |
684 |
|
raynormal(pnorm, r); |
685 |
< |
(*putreal)(pnorm[0]); |
621 |
< |
(*putreal)(pnorm[1]); |
622 |
< |
(*putreal)(pnorm[2]); |
685 |
> |
(*putreal)(pnorm, 3); |
686 |
|
} |
687 |
|
|
688 |
|
|
704 |
|
RAY *r |
705 |
|
) |
706 |
|
{ |
707 |
< |
(*putreal)(r->rweight); |
707 |
> |
RREAL rwt = r->rweight; |
708 |
> |
|
709 |
> |
(*putreal)(&rwt, 1); |
710 |
|
} |
711 |
|
|
712 |
|
|
715 |
|
RAY *r |
716 |
|
) |
717 |
|
{ |
718 |
< |
double contr[3]; |
718 |
> |
RREAL contr[3]; |
719 |
> |
/* shadow ray not on source? */ |
720 |
> |
if (r->rsrc >= 0 && source[r->rsrc].so != r->ro) |
721 |
> |
setcolor(contr, 0.0, 0.0, 0.0); |
722 |
> |
else |
723 |
> |
raycontrib(contr, r, PRIMARY); |
724 |
|
|
725 |
< |
raycontrib(contr, r, PRIMARY); |
656 |
< |
(*putreal)(contr[RED]); |
657 |
< |
(*putreal)(contr[GRN]); |
658 |
< |
(*putreal)(contr[BLU]); |
725 |
> |
(*putreal)(contr, 3); |
726 |
|
} |
727 |
|
|
728 |
|
|
770 |
|
|
771 |
|
|
772 |
|
static void |
773 |
< |
puta( /* print ascii value */ |
774 |
< |
double v |
773 |
> |
puta( /* print ascii value(s) */ |
774 |
> |
RREAL *v, int n |
775 |
|
) |
776 |
|
{ |
777 |
< |
printf("%e\t", v); |
777 |
> |
if (n == 3) { |
778 |
> |
printf("%e\t%e\t%e\t", v[0], v[1], v[2]); |
779 |
> |
return; |
780 |
> |
} |
781 |
> |
while (n--) |
782 |
> |
printf("%e\t", *v++); |
783 |
|
} |
784 |
|
|
785 |
|
|
786 |
|
static void |
787 |
< |
putd(v) /* print binary double */ |
716 |
< |
double v; |
787 |
> |
putd(RREAL *v, int n) /* print binary double(s) */ |
788 |
|
{ |
789 |
< |
fwrite((char *)&v, sizeof(v), 1, stdout); |
789 |
> |
#ifdef SMLFLT |
790 |
> |
double da[3]; |
791 |
> |
int i; |
792 |
> |
|
793 |
> |
if (n > 3) |
794 |
> |
error(INTERNAL, "code error in putd()"); |
795 |
> |
for (i = n; i--; ) |
796 |
> |
da[i] = v[i]; |
797 |
> |
putbinary(da, sizeof(double), n, stdout); |
798 |
> |
#else |
799 |
> |
putbinary(v, sizeof(RREAL), n, stdout); |
800 |
> |
#endif |
801 |
|
} |
802 |
|
|
803 |
|
|
804 |
|
static void |
805 |
< |
putf(v) /* print binary float */ |
724 |
< |
double v; |
805 |
> |
putf(RREAL *v, int n) /* print binary float(s) */ |
806 |
|
{ |
807 |
< |
float f = v; |
807 |
> |
#ifndef SMLFLT |
808 |
> |
float fa[3]; |
809 |
> |
int i; |
810 |
|
|
811 |
< |
fwrite((char *)&f, sizeof(f), 1, stdout); |
811 |
> |
if (n > 3) |
812 |
> |
error(INTERNAL, "code error in putf()"); |
813 |
> |
for (i = n; i--; ) |
814 |
> |
fa[i] = v[i]; |
815 |
> |
putbinary(fa, sizeof(float), n, stdout); |
816 |
> |
#else |
817 |
> |
putbinary(v, sizeof(RREAL), n, stdout); |
818 |
> |
#endif |
819 |
|
} |