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

Comparing ray/src/rt/rcontrib.c (file contents):
Revision 2.7 by greg, Fri Jun 15 21:32:11 2012 UTC vs.
Revision 2.16 by greg, Tue Aug 14 22:25:12 2012 UTC

# Line 7 | Line 7 | static const char RCSid[] = "$Id$";
7   */
8  
9   #include "rcontrib.h"
10 – #include "source.h"
10   #include "otypes.h"
11 < #include "platform.h"
11 > #include "source.h"
12  
13   char    *shm_boundary = NULL;           /* boundary of shared memory */
14  
# Line 27 | Line 26 | int    do_irrad = 0;                   /* compute irradiance? */
26  
27   int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
28  
29 < double  dstrsrc = 0.0;                  /* square source distribution */
29 > double  dstrsrc = 0.9;                  /* square source distribution */
30   double  shadthresh = .03;               /* shadow threshold */
31   double  shadcert = .75;                 /* shadow certainty */
32   int     directrelay = 3;                /* number of source relays */
# Line 137 | Line 136 | addmodifier(char *modn, char *outf, char *binv, int bi
136   }
137  
138  
139 < /* add modifiers from a file list */
139 > /* Add modifiers from a file list */
140   void
141   addmodfile(char *fname, char *outf, char *binv, int bincnt)
142   {
# Line 159 | Line 158 | quit(                  /* quit program */
158   )
159   {
160          if (nchild > 0)         /* close children if any */
161 <                end_children();
161 >                end_children(code != 0);
162          exit(code);
163   }
164  
# Line 178 | Line 177 | rcinit()
177                                          /* set shared memory boundary */
178                  shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY");
179          }
180 <        if ((nproc > 1) & (accumulate <= 0))
181 <                put_zero_record(0);     /* prime our queue to accumulate */
183 <
180 >        for (i = 0; i < nsources; i++)  /* tracing to sources as well */
181 >                source[i].sflags |= SFOLLOW;
182          if (yres > 0) {                 /* set up flushing & ray counts */
183                  if (xres > 0)
184                          raysleft = (RNUMBER)xres*yres;
# Line 191 | Line 189 | rcinit()
189          if ((account = accumulate) > 1)
190                  raysleft *= accumulate;
191          waitflush = (yres > 0) & (xres > 1) ? 0 : xres;
194 –                                        /* tracing to sources as well */
195 –        for (i = 0; i < nsources; i++)
196 –                source[i].sflags |= SFOLLOW;
192  
193          if (nproc > 1 && in_rchild())   /* forked child? */
194                  return;                 /* return to main processing loop */
195  
196          if (recover) {                  /* recover previous output? */
197 <                if (accumulate <= 0) {
197 >                if (accumulate <= 0)
198                          reload_output();
199 <                        if (nproc > 1)
205 <                                queue_modifiers();
206 <                } else
199 >                else
200                          recover_output();
201          }
202          if (nproc == 1)                 /* single process? */
203                  return;
204 <
205 <        parental_loop();                /* else run controller */
204 >                                        /* else run appropriate controller */
205 >        if (accumulate <= 0)
206 >                feeder_loop();
207 >        else
208 >                parental_loop();
209          quit(0);                        /* parent musn't return! */
210   }
211  
# Line 230 | Line 226 | trace_contrib(RAY *r)
226  
227          if (mp == NULL)                         /* not in our list? */
228                  return;
229 +                                                /* shadow ray not on source? */
230 +        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
231 +                return;
232  
233 <        worldfunc(RCCONTEXT, r);                /* get bin number */
233 >        worldfunc(RCCONTEXT, r);                /* else get bin number */
234          bn = (int)(evalue(mp->binv) + .5);
235          if ((bn < 0) | (bn >= mp->nbins)) {
236                  error(WARNING, "bad bin number (ignored)");
237                  return;
238          }
239 <        raycontrib(contr, r, PRIMARY);
239 >        raycontrib(contr, r, PRIMARY);          /* compute coefficient */
240          if (contrib)
241 <                multcolor(contr, r->rcol);
241 >                multcolor(contr, r->rcol);      /* -> contribution */
242          addcolor(mp->cbin[bn], contr);
243   }
244  
# Line 256 | Line 255 | eval_irrad(FVECT org, FVECT dir)
255          thisray.rdir[2] = -dir[2];
256          thisray.rmax = 0.0;
257          rayorigin(&thisray, PRIMARY, NULL, NULL);
258 <        thisray.rot = 1e-5;             /* pretend we hit surface */
258 >                                        /* pretend we hit surface */
259 >        thisray.rt = thisray.rot = 1e-5;
260          thisray.rod = 1.0;
261 +        VCOPY(thisray.ron, dir);
262          VSUM(thisray.rop, org, dir, 1e-4);
263          samplendx++;                    /* compute result */
264          (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
# Line 346 | Line 347 | rcontrib()
347                  account = 1;            /* output accumulated totals */
348                  done_contrib();
349          }
350 +        lu_done(&ofiletab);             /* close output files */
351          if (raysleft)
352                  error(USER, "unexpected EOF on input");
351 –        lu_done(&ofiletab);             /* close output files */
353   }

Diff Legend

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