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.39 by greg, Thu Apr 14 18:04:12 2005 UTC vs.
Revision 2.50 by greg, Fri Oct 7 03:45:14 2005 UTC

# Line 29 | Line 29 | static const char      RCSid[] = "$Id$";
29   #include  "source.h"
30   #include  "otypes.h"
31   #include  "resolu.h"
32 + #include  "random.h"
33  
34   CUBE  thescene;                         /* our scene */
35   OBJECT  nsceneobjs;                     /* number of objects in our scene */
# Line 46 | Line 47 | char  *outvals = "v";                  /* output specification */
47  
48   int  do_irrad = 0;                      /* compute irradiance? */
49  
50 + int  rand_samp = 0;                     /* pure Monte Carlo sampling? */
51 +
52   void  (*trace)() = NULL;                /* trace call */
53  
54   char  *tralist[128];                    /* list of modifers to trace (or no) */
55   int  traincl = -1;                      /* include == 1, exclude == 0 */
56 < #define  MAXTSET        511             /* maximum number in trace set */
56 > #ifndef  MAXTSET
57 > #define  MAXTSET        1024            /* maximum number in trace set */
58 > #endif
59   OBJECT  traset[MAXTSET+1]={0};          /* trace include/exclude set */
60  
61   int  hresolu = 0;                       /* horizontal (scan) size */
# Line 74 | Line 79 | double  specjitter = 1.;               /* specular sampling jitter
79  
80   int  backvis = 1;                       /* back face visibility */
81  
82 < int  maxdepth = 8;                      /* maximum recursion depth */
82 > int  maxdepth = -10;                    /* maximum recursion depth */
83   double  minweight = 2e-3;               /* minimum ray weight */
84  
85   char  *ambfile = NULL;                  /* ambient file name */
# Line 96 | Line 101 | typedef void putf_t(double v);
101   static putf_t puta, putd, putf;
102  
103   typedef void oputf_t(RAY *r);
104 < static oputf_t  oputo, oputd, oputv, oputl, oputL, oputc,
105 <                oputp, oputn, oputN, oputs, oputw, oputm, oputM;
104 > static oputf_t  oputo, oputd, oputv, oputl, oputL, oputc, oputp,
105 >                oputn, oputN, oputs, oputw, oputW, oputm, oputM, oputtilde;
106  
107   static void setoutput(char *vs);
108   static void tranotify(OBJECT obj);
# Line 148 | Line 153 | rtrace(                                /* trace rays from file */
153          char  *fname
154   )
155   {
156 <        long  vcount = hresolu>1 ? hresolu*vresolu : vresolu;
156 >        unsigned long  vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu
157 >                                              : vresolu;
158          long  nextflush = hresolu;
159          FILE  *fp;
160          double  d;
# Line 160 | Line 166 | rtrace(                                /* trace rays from file */
166                  sprintf(errmsg, "cannot open input file \"%s\"", fname);
167                  error(SYSTEM, errmsg);
168          }
163 – #ifdef _WIN32
169          if (inform != 'a')
170                  SET_FILE_BINARY(fp);
166 – #endif
171                                          /* set up output */
172          setoutput(outvals);
173          switch (outform) {
# Line 189 | Line 193 | rtrace(                                /* trace rays from file */
193                  d = normalize(direc);
194                  if (d == 0.0) {                         /* zero ==> flush */
195                          bogusray();
196 <                        if (--nextflush <= 0 || vcount <= 0) {
196 >                        if (--nextflush <= 0 || !vcount) {
197                                  fflush(stdout);
198                                  nextflush = hresolu;
199                          }
200                  } else {
201 <                        samplendx++;
201 >                        samplendx = rand_samp ? random() : samplendx+1;
202                                                          /* compute and print */
203                          if (imm_irrad)
204                                  irrad(orig, direc);
205                          else
206                                  rad(orig, direc, lim_dist ? d : 0.0);
207                                                          /* flush if time */
208 <                        if (--nextflush == 0) {
208 >                        if (!--nextflush) {
209                                  fflush(stdout);
210                                  nextflush = hresolu;
211                          }
212                  }
213                  if (ferror(stdout))
214                          error(SYSTEM, "write error");
215 <                if (--vcount == 0)                      /* check for end */
215 >                if (vcount && !--vcount)                /* check for end */
216                          break;
217          }
218 <        fflush(stdout);
219 <        if (vcount > 0)
220 <                error(USER, "read error");
218 >        if (fflush(stdout) < 0)
219 >                error(SYSTEM, "write error");
220 >        if (vcount)
221 >                error(USER, "unexpected EOF on input");
222          if (fname != NULL)
223                  fclose(fp);
224   }
225  
226  
227   static void
228 + trace_sources(void)                     /* trace rays to light sources, also */
229 + {
230 +        int     sn;
231 +        
232 +        for (sn = 0; sn < nsources; sn++)
233 +                source[sn].sflags |= SFOLLOW;
234 + }
235 +
236 +
237 + static void
238   setoutput(                              /* set up output tables */
239          register char  *vs
240   )
# Line 229 | Line 244 | setoutput(                             /* set up output tables */
244          castonly = 1;
245          while (*vs)
246                  switch (*vs++) {
247 +                case 'T':                               /* trace sources */
248 +                        if (!*vs) break;
249 +                        trace_sources();
250 +                        /* fall through */
251                  case 't':                               /* trace */
252 +                        if (!*vs) break;
253                          *table = NULL;
254                          table = every_out;
255                          trace = ourtrace;
# Line 271 | Line 291 | setoutput(                             /* set up output tables */
291                  case 'w':                               /* weight */
292                          *table++ = oputw;
293                          break;
294 +                case 'W':                               /* coefficient */
295 +                        *table++ = oputW;
296 +                        if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0))
297 +                                error(WARNING,
298 +                                        "-otW accuracy depends on -aa 0 -as 0");
299 +                        break;
300                  case 'm':                               /* modifier */
301                          *table++ = oputm;
302                          break;
303                  case 'M':                               /* material */
304                          *table++ = oputM;
305                          break;
306 +                case '~':                               /* tilde */
307 +                        *table++ = oputtilde;
308 +                        break;
309                  }
310          *table = NULL;
311   }
# Line 287 | Line 316 | bogusray(void)                 /* print out empty record */
316   {
317          thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] =
318          thisray.rdir[0] = thisray.rdir[1] = thisray.rdir[2] = 0.0;
319 <        rayorigin(&thisray, NULL, PRIMARY, 1.0);
319 >        thisray.rmax = 0.0;
320 >        rayorigin(&thisray, PRIMARY, NULL, NULL);
321          printvals(&thisray);
322   }
323  
# Line 302 | Line 332 | rad(           /* compute and print ray value(s) */
332          VCOPY(thisray.rorg, org);
333          VCOPY(thisray.rdir, dir);
334          thisray.rmax = dmax;
335 <        rayorigin(&thisray, NULL, PRIMARY, 1.0);
335 >        rayorigin(&thisray, PRIMARY, NULL, NULL);
336          if (castonly) {
337                  if (!localhit(&thisray, &thescene)) {
338                          if (thisray.ro == &Aftplane) {  /* clipped */
# Line 329 | Line 359 | irrad(                 /* compute immediate irradiance value */
359                  thisray.rorg[i] = org[i] + dir[i];
360                  thisray.rdir[i] = -dir[i];
361          }
362 <        rayorigin(&thisray, NULL, PRIMARY, 1.0);
362 >        thisray.rmax = 0.0;
363 >        rayorigin(&thisray, PRIMARY, NULL, NULL);
364                                          /* pretend we hit surface */
365          thisray.rot = 1.0-1e-4;
366          thisray.rod = 1.0;
# Line 442 | Line 473 | ourtrace(                              /* print ray values */
473          tabin(r);
474          for (tp = every_out; *tp != NULL; tp++)
475                  (**tp)(r);
476 <        putchar('\n');
476 >        if (outform == 'a')
477 >                putchar('\n');
478   }
479  
480  
# Line 451 | Line 483 | tabin(                         /* tab in appropriate amount */
483          RAY  *r
484   )
485   {
486 <        register RAY  *rp;
486 >        const RAY  *rp;
487  
488          for (rp = r->parent; rp != NULL; rp = rp->parent)
489                  putchar('\t');
# Line 485 | Line 517 | oputv(                         /* print value */
517          RAY  *r
518   )
519   {
488 –        COLR  cout;
489 –        
520          if (outform == 'c') {
521 +                COLR  cout;
522                  setcolr(cout,   colval(r->rcol,RED),
523                                  colval(r->rcol,GRN),
524                                  colval(r->rcol,BLU));
# Line 605 | Line 636 | oputw(                         /* print weight */
636  
637  
638   static void
639 + oputW(                          /* print contribution */
640 +        RAY  *r
641 + )
642 + {
643 +        double  contr[3];
644 +
645 +        raycontrib(contr, r, PRIMARY);
646 +        (*putreal)(contr[RED]);
647 +        (*putreal)(contr[GRN]);
648 +        (*putreal)(contr[BLU]);
649 + }
650 +
651 +
652 + static void
653   oputm(                          /* print modifier */
654          RAY  *r
655   )
# Line 635 | Line 680 | oputM(                         /* print material */
680          } else
681                  putchar('*');
682          putchar('\t');
683 + }
684 +
685 +
686 + static void
687 + oputtilde(                      /* output tilde (spacer) */
688 +        RAY  *r
689 + )
690 + {
691 +        fputs("~\t", stdout);
692   }
693  
694  

Diff Legend

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