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.58 by greg, Sat Dec 12 23:08:13 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 63 | Line 64 | static oputf_t  oputo, oputd, oputv, oputV, oputl, opu
64   static void setoutput(char *vs);
65   extern void tranotify(OBJECT obj);
66   static void bogusray(void);
67 < static void rad(FVECT org, FVECT dir, double dmax);
68 < static void irrad(FVECT org, FVECT dir);
67 > static void raycast(RAY *r);
68 > static void rayirrad(RAY *r);
69 > static void rtcompute(FVECT org, FVECT dir, double dmax);
70   static int printvals(RAY *r);
71   static int getvec(FVECT vec, int fmt, FILE *fp);
72   static void tabin(RAY *r);
# Line 108 | Line 110 | formstr(                               /* return format identifier */
110  
111   extern void
112   rtrace(                         /* trace rays from file */
113 <        char  *fname
113 >        char  *fname,
114 >        int  nproc
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 128 | Line 131 | rtrace(                                /* trace rays from file */
131                  SET_FILE_BINARY(fp);
132                                          /* set up output */
133          setoutput(outvals);
134 +        if (imm_irrad)
135 +                castonly = 0;
136 +        else if (castonly)
137 +                nproc = 1;              /* don't bother multiprocessing */
138          switch (outform) {
139          case 'a': putreal = puta; break;
140          case 'f': putreal = putf; break;
# Line 139 | Line 146 | rtrace(                                /* trace rays from file */
146          default:
147                  error(CONSISTENCY, "botched output format");
148          }
149 <        if (ray_pnprocs > 1)
149 >        if (nproc > 1) {                /* start multiprocessing */
150 >                ray_popen(nproc);
151                  ray_fifo_out = printvals;
152 +        }
153          if (hresolu > 0) {
154                  if (vresolu > 0)
155                          fprtresolu(hresolu, vresolu, stdout);
# Line 152 | Line 161 | rtrace(                                /* trace rays from file */
161  
162                  d = normalize(direc);
163                  if (d == 0.0) {                         /* zero ==> flush */
155 –                        if (ray_pnprocs > 1 && ray_fifo_flush() < 0)
156 –                                error(USER, "lost children");
157 –                        bogusray();
164                          if (--nextflush <= 0 || !vcount) {
165 +                                if (nproc > 1 && ray_fifo_flush() < 0)
166 +                                        error(USER, "child(ren) died");
167 +                                bogusray();
168                                  fflush(stdout);
169 <                                nextflush = hresolu;
170 <                        }
169 >                                nextflush = (vresolu > 0) & (hresolu > 1) ? 0 :
170 >                                                                hresolu;
171 >                        } else
172 >                                bogusray();
173                  } else {                                /* compute and print */
174 <                        if (imm_irrad)
164 <                                irrad(orig, direc);
165 <                        else
166 <                                rad(orig, direc, lim_dist ? d : 0.0);
174 >                        rtcompute(orig, direc, lim_dist ? d : 0.0);
175                                                          /* flush if time */
176                          if (!--nextflush) {
177 <                                if (ray_pnprocs > 1 && ray_fifo_flush() < 0)
178 <                                        error(USER, "lost children");
177 >                                if (nproc > 1 && ray_fifo_flush() < 0)
178 >                                        error(USER, "child(ren) died");
179                                  fflush(stdout);
180                                  nextflush = hresolu;
181                          }
# Line 177 | Line 185 | rtrace(                                /* trace rays from file */
185                  if (vcount && !--vcount)                /* check for end */
186                          break;
187          }
188 <        if (ray_pnprocs > 1 && ray_fifo_flush() < 0)
189 <                error(USER, "unable to complete processing");
188 >        if (nproc > 1) {                                /* clean up children */
189 >                if (ray_fifo_flush() < 0)
190 >                        error(USER, "unable to complete processing");
191 >                ray_pclose(0);
192 >        }
193          if (fflush(stdout) < 0)
194                  error(SYSTEM, "write error");
195          if (vcount)
# Line 293 | Line 304 | bogusray(void)                 /* print out empty record */
304  
305  
306   static void
307 < rad(            /* compute and print ray value(s) */
308 <        FVECT  org,
298 <        FVECT  dir,
299 <        double  dmax
307 > raycast(                        /* compute first ray intersection only */
308 >        RAY *r
309   )
310   {
311 <        VCOPY(thisray.rorg, org);
312 <        VCOPY(thisray.rdir, dir);
313 <        thisray.rmax = dmax;
314 <        if (!castonly && ray_pnprocs > 1) {
315 <                if (ray_fifo_in(&thisray) < 0)
316 <                        error(USER, "lost children");
308 <                return;
311 >        if (!localhit(r, &thescene)) {
312 >                if (r->ro == &Aftplane) {       /* clipped */
313 >                        r->ro = NULL;
314 >                        r->rot = FHUGE;
315 >                } else
316 >                        sourcehit(r);
317          }
310 –        rayorigin(&thisray, PRIMARY, NULL, NULL);
311 –        if (castonly) {
312 –                if (!localhit(&thisray, &thescene)) {
313 –                        if (thisray.ro == &Aftplane) {  /* clipped */
314 –                                thisray.ro = NULL;
315 –                                thisray.rot = FHUGE;
316 –                        } else
317 –                                sourcehit(&thisray);
318 –                }
319 –        } else
320 –                ray_trace(&thisray);
321 –        printvals(&thisray);
318   }
319  
320  
321   static void
322 < irrad(                  /* compute immediate irradiance value */
322 > rayirrad(                       /* compute irradiance rather than radiance */
323 >        RAY *r
324 > )
325 > {
326 >        void    (*old_revf)(RAY *) = r->revf;
327 >
328 >        r->rot = 1e-5;                  /* pretend we hit surface */
329 >        VSUM(r->rop, r->rorg, r->rdir, r->rot);
330 >        r->ron[0] = -r->rdir[0];
331 >        r->ron[1] = -r->rdir[1];
332 >        r->ron[2] = -r->rdir[2];
333 >        r->rod = 1.0;
334 >                                        /* compute result */
335 >        r->revf = raytrace;
336 >        (*ofun[Lamb.otype].funp)(&Lamb, r);
337 >        r->revf = old_revf;
338 > }
339 >
340 >
341 > static void
342 > rtcompute(                      /* compute and print ray value(s) */
343          FVECT  org,
344 <        FVECT  dir
344 >        FVECT  dir,
345 >        double  dmax
346   )
347   {
348 <        VSUM(thisray.rorg, org, dir, 1.1e-4);
332 <        thisray.rdir[0] = -dir[0];
333 <        thisray.rdir[1] = -dir[1];
334 <        thisray.rdir[2] = -dir[2];
335 <        thisray.rmax = 0.0;
348 >                                        /* set up ray */
349          rayorigin(&thisray, PRIMARY, NULL, NULL);
350 <                                        /* pretend we hit surface */
351 <        thisray.rot = 1e-5;
352 <        thisray.rod = 1.0;
353 <        VCOPY(thisray.ron, dir);
354 <        VSUM(thisray.rop, org, dir, 1e-4);
355 <                                        /* compute and print */
356 <        (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
350 >        if (imm_irrad) {
351 >                VSUM(thisray.rorg, org, dir, 1.1e-4);
352 >                thisray.rdir[0] = -dir[0];
353 >                thisray.rdir[1] = -dir[1];
354 >                thisray.rdir[2] = -dir[2];
355 >                thisray.rmax = 0.0;
356 >                thisray.revf = rayirrad;
357 >        } else {
358 >                VCOPY(thisray.rorg, org);
359 >                VCOPY(thisray.rdir, dir);
360 >                thisray.rmax = dmax;
361 >                if (castonly)
362 >                        thisray.revf = raycast;
363 >        }
364 >        if (ray_pnprocs > 1) {          /* multiprocessing FIFO? */
365 >                if (ray_fifo_in(&thisray) < 0)
366 >                        error(USER, "lost children");
367 >                return;
368 >        }
369 >        samplendx++;                    /* else do it ourselves */
370 >        rayvalue(&thisray);
371          printvals(&thisray);
372   }
373  
# Line 384 | 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 468 | Line 495 | oputo(                         /* print origin */
495          RAY  *r
496   )
497   {
498 <        (*putreal)(r->rorg[0]);
472 <        (*putreal)(r->rorg[1]);
473 <        (*putreal)(r->rorg[2]);
498 >        (*putreal)(r->rorg, 3);
499   }
500  
501  
# Line 479 | Line 504 | oputd(                         /* print direction */
504          RAY  *r
505   )
506   {
507 <        (*putreal)(r->rdir[0]);
483 <        (*putreal)(r->rdir[1]);
484 <        (*putreal)(r->rdir[2]);
507 >        (*putreal)(r->rdir, 3);
508   }
509  
510  
# Line 490 | 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 509 | 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]);
517 <        (*putreal)(contr[GRN]);
518 <        (*putreal)(contr[BLU]);
542 >        (*putreal)(contr, 3);
543   }
544  
545  
# Line 524 | Line 548 | oputl(                         /* print effective distance */
548          RAY  *r
549   )
550   {
551 <        (*putreal)(r->rt);
551 >        (*putreal)(&r->rt, 1);
552   }
553  
554  
# Line 533 | 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 542 | Line 566 | oputc(                         /* print local coordinates */
566          RAY  *r
567   )
568   {
569 <        (*putreal)(r->uv[0]);
546 <        (*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]);
559 <        } else {
560 <                (*putreal)(0.0);
561 <                (*putreal)(0.0);
562 <                (*putreal)(0.0);
563 <        }
581 >        if (r->rot < FHUGE)
582 >                (*putreal)(r->rop, 3);
583 >        else
584 >                (*putreal)(vdummy, 3);
585   }
586  
587  
# Line 569 | 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]);
576 <        } else {
577 <                (*putreal)(0.0);
578 <                (*putreal)(0.0);
579 <                (*putreal)(0.0);
580 <        }
593 >        if (r->rot < FHUGE)
594 >                (*putreal)(r->ron, 3);
595 >        else
596 >                (*putreal)(vdummy, 3);
597   }
598  
599  
# Line 589 | Line 605 | oputn(                         /* print perturbed normal */
605          FVECT  pnorm;
606  
607          if (r->rot >= FHUGE) {
608 <                (*putreal)(0.0);
593 <                (*putreal)(0.0);
594 <                (*putreal)(0.0);
608 >                (*putreal)(vdummy, 3);
609                  return;
610          }
611          raynormal(pnorm, r);
612 <        (*putreal)(pnorm[0]);
599 <        (*putreal)(pnorm[1]);
600 <        (*putreal)(pnorm[2]);
612 >        (*putreal)(pnorm, 3);
613   }
614  
615  
# Line 619 | 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 628 | 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);
634 <        (*putreal)(contr[RED]);
635 <        (*putreal)(contr[GRN]);
636 <        (*putreal)(contr[BLU]);
652 >        (*putreal)(contr, 3);
653   }
654  
655  
# Line 681 | 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 */
694 < 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 */
702 < 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