ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rtrace.c
(Generate patch)

Comparing ray/src/rt/rtrace.c (file contents):
Revision 2.63 by greg, Wed Dec 16 03:30:50 2009 UTC vs.
Revision 2.72 by greg, Thu Nov 8 00:54:07 2018 UTC

# Line 28 | Line 28 | static const char      RCSid[] = "$Id$";
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  
# Line 53 | Line 54 | OBJECT traset[MAXTSET+1]={0};          /* trace include/exclud
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);
# Line 114 | Line 115 | rtrace(                                /* trace rays from file */
115   )
116   {
117          unsigned long  vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu
118 <                                              : vresolu;
119 <        long  nextflush = hresolu;
118 >                                              : (unsigned long)vresolu;
119 >        long  nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : hresolu;
120          FILE  *fp;
121          double  d;
122          FVECT  orig, direc;
# Line 162 | Line 163 | rtrace(                                /* trace rays from file */
163                  if (d == 0.0) {                         /* zero ==> flush */
164                          if (--nextflush <= 0 || !vcount) {
165                                  if (nproc > 1 && ray_fifo_flush() < 0)
166 <                                        error(USER, "lost children");
166 >                                        error(USER, "child(ren) died");
167                                  bogusray();
168                                  fflush(stdout);
169 <                                nextflush = hresolu;
169 >                                nextflush = (vresolu > 0) & (hresolu > 1) ? 0 :
170 >                                                                hresolu;
171                          } else
172                                  bogusray();
173                  } else {                                /* compute and print */
# Line 173 | Line 175 | rtrace(                                /* trace rays from file */
175                                                          /* flush if time */
176                          if (!--nextflush) {
177                                  if (nproc > 1 && ray_fifo_flush() < 0)
178 <                                        error(USER, "lost children");
178 >                                        error(USER, "child(ren) died");
179                                  fflush(stdout);
180                                  nextflush = hresolu;
181                          }
# Line 409 | Line 411 | getvec(                /* get a vector from fp */
411                  }
412                  break;
413          case 'f':                                       /* binary float */
414 <                if (fread((char *)vf, sizeof(float), 3, fp) != 3)
414 >                if (getbinary(vf, sizeof(float), 3, fp) != 3)
415                          return(-1);
416                  VCOPY(vec, vf);
417                  break;
418          case 'd':                                       /* binary double */
419 <                if (fread((char *)vd, sizeof(double), 3, fp) != 3)
419 >                if (getbinary(vd, sizeof(double), 3, fp) != 3)
420                          return(-1);
421                  VCOPY(vec, vd);
422                  break;
# Line 493 | Line 495 | oputo(                         /* print origin */
495          RAY  *r
496   )
497   {
498 <        (*putreal)(r->rorg[0]);
497 <        (*putreal)(r->rorg[1]);
498 <        (*putreal)(r->rorg[2]);
498 >        (*putreal)(r->rorg, 3);
499   }
500  
501  
# Line 504 | Line 504 | oputd(                         /* print direction */
504          RAY  *r
505   )
506   {
507 <        (*putreal)(r->rdir[0]);
508 <        (*putreal)(r->rdir[1]);
509 <        (*putreal)(r->rdir[2]);
507 >        (*putreal)(r->rdir, 3);
508   }
509  
510  
# Line 515 | Line 513 | oputv(                         /* print value */
513          RAY  *r
514   )
515   {
516 +        RREAL   cval[3];
517 +
518          if (outform == 'c') {
519                  COLR  cout;
520                  setcolr(cout,   colval(r->rcol,RED),
521                                  colval(r->rcol,GRN),
522                                  colval(r->rcol,BLU));
523 <                fwrite((char *)cout, sizeof(cout), 1, stdout);
523 >                putbinary(cout, sizeof(cout), 1, stdout);
524                  return;
525          }
526 <        (*putreal)(colval(r->rcol,RED));
527 <        (*putreal)(colval(r->rcol,GRN));
528 <        (*putreal)(colval(r->rcol,BLU));
526 >        cval[0] = colval(r->rcol,RED);
527 >        cval[1] = colval(r->rcol,GRN);
528 >        cval[2] = colval(r->rcol,BLU);
529 >        (*putreal)(cval, 3);
530   }
531  
532  
# Line 534 | Line 535 | oputV(                         /* print value contribution */
535          RAY *r
536   )
537   {
538 <        double  contr[3];
538 >        RREAL   contr[3];
539  
540          raycontrib(contr, r, PRIMARY);
541          multcolor(contr, r->rcol);
542 <        (*putreal)(contr[RED]);
542 <        (*putreal)(contr[GRN]);
543 <        (*putreal)(contr[BLU]);
542 >        (*putreal)(contr, 3);
543   }
544  
545  
# Line 549 | Line 548 | oputl(                         /* print effective distance */
548          RAY  *r
549   )
550   {
551 <        (*putreal)(r->rt);
551 >        (*putreal)(&r->rt, 1);
552   }
553  
554  
# Line 558 | Line 557 | oputL(                         /* print single ray length */
557          RAY  *r
558   )
559   {
560 <        (*putreal)(r->rot);
560 >        (*putreal)(&r->rot, 1);
561   }
562  
563  
# Line 567 | Line 566 | oputc(                         /* print local coordinates */
566          RAY  *r
567   )
568   {
569 <        (*putreal)(r->uv[0]);
571 <        (*putreal)(r->uv[1]);
569 >        (*putreal)(r->uv, 2);
570   }
571  
572  
573 + static RREAL    vdummy[3] = {0.0, 0.0, 0.0};
574 +
575 +
576   static void
577   oputp(                          /* print point */
578          RAY  *r
579   )
580   {
581 <        if (r->rot < FHUGE) {
582 <                (*putreal)(r->rop[0]);
583 <                (*putreal)(r->rop[1]);
584 <                (*putreal)(r->rop[2]);
584 <        } else {
585 <                (*putreal)(0.0);
586 <                (*putreal)(0.0);
587 <                (*putreal)(0.0);
588 <        }
581 >        if (r->rot < FHUGE)
582 >                (*putreal)(r->rop, 3);
583 >        else
584 >                (*putreal)(vdummy, 3);
585   }
586  
587  
# Line 594 | Line 590 | oputN(                         /* print unperturbed normal */
590          RAY  *r
591   )
592   {
593 <        if (r->rot < FHUGE) {
594 <                (*putreal)(r->ron[0]);
595 <                (*putreal)(r->ron[1]);
596 <                (*putreal)(r->ron[2]);
601 <        } else {
602 <                (*putreal)(0.0);
603 <                (*putreal)(0.0);
604 <                (*putreal)(0.0);
605 <        }
593 >        if (r->rot < FHUGE)
594 >                (*putreal)(r->ron, 3);
595 >        else
596 >                (*putreal)(vdummy, 3);
597   }
598  
599  
# Line 614 | Line 605 | oputn(                         /* print perturbed normal */
605          FVECT  pnorm;
606  
607          if (r->rot >= FHUGE) {
608 <                (*putreal)(0.0);
618 <                (*putreal)(0.0);
619 <                (*putreal)(0.0);
608 >                (*putreal)(vdummy, 3);
609                  return;
610          }
611          raynormal(pnorm, r);
612 <        (*putreal)(pnorm[0]);
624 <        (*putreal)(pnorm[1]);
625 <        (*putreal)(pnorm[2]);
612 >        (*putreal)(pnorm, 3);
613   }
614  
615  
# Line 644 | Line 631 | oputw(                         /* print weight */
631          RAY  *r
632   )
633   {
634 <        (*putreal)(r->rweight);
634 >        RREAL   rwt = r->rweight;
635 >        
636 >        (*putreal)(&rwt, 1);
637   }
638  
639  
# Line 653 | Line 642 | oputW(                         /* print coefficient */
642          RAY  *r
643   )
644   {
645 <        double  contr[3];
645 >        RREAL   contr[3];
646 >                                /* shadow ray not on source? */
647 >        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
648 >                setcolor(contr, 0.0, 0.0, 0.0);
649 >        else
650 >                raycontrib(contr, r, PRIMARY);
651  
652 <        raycontrib(contr, r, PRIMARY);
659 <        (*putreal)(contr[RED]);
660 <        (*putreal)(contr[GRN]);
661 <        (*putreal)(contr[BLU]);
652 >        (*putreal)(contr, 3);
653   }
654  
655  
# Line 706 | Line 697 | oputtilde(                     /* output tilde (spacer) */
697  
698  
699   static void
700 < puta(                           /* print ascii value */
701 <        double  v
700 > puta(                           /* print ascii value(s) */
701 >        RREAL *v, int n
702   )
703   {
704 <        printf("%e\t", v);
704 >        if (n == 3) {
705 >                printf("%e\t%e\t%e\t", v[0], v[1], v[2]);
706 >                return;
707 >        }
708 >        while (n--)
709 >                printf("%e\t", *v++);
710   }
711  
712  
713   static void
714 < putd(v)                         /* print binary double */
719 < double  v;
714 > putd(RREAL *v, int n)           /* print binary double(s) */
715   {
716 <        fwrite((char *)&v, sizeof(v), 1, stdout);
716 > #ifdef  SMLFLT
717 >        double  da[3];
718 >        int     i;
719 >
720 >        if (n > 3)
721 >                error(INTERNAL, "code error in putd()");
722 >        for (i = n; i--; )
723 >                da[i] = v[i];
724 >        putbinary(da, sizeof(double), n, stdout);
725 > #else
726 >        putbinary(v, sizeof(RREAL), n, stdout);
727 > #endif
728   }
729  
730  
731   static void
732 < putf(v)                         /* print binary float */
727 < double  v;
732 > putf(RREAL *v, int n)           /* print binary float(s) */
733   {
734 <        float f = v;
734 > #ifndef SMLFLT
735 >        float   fa[3];
736 >        int     i;
737  
738 <        fwrite((char *)&f, sizeof(f), 1, stdout);
738 >        if (n > 3)
739 >                error(INTERNAL, "code error in putf()");
740 >        for (i = n; i--; )
741 >                fa[i] = v[i];
742 >        putbinary(fa, sizeof(float), n, stdout);
743 > #else
744 >        putbinary(v, sizeof(RREAL), n, stdout);
745 > #endif
746   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines