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.8 by greg, Sat Jun 16 17:09:49 2012 UTC vs.
Revision 2.25 by greg, Wed May 20 13:16:20 2015 UTC

# Line 6 | Line 6 | static const char RCSid[] = "$Id$";
6   * Initialization and calculation routines
7   */
8  
9 + #include "copyright.h"
10 +
11   #include "rcontrib.h"
10 – #include "source.h"
12   #include "otypes.h"
13 < #include "platform.h"
13 > #include "source.h"
14  
15   char    *shm_boundary = NULL;           /* boundary of shared memory */
16  
# Line 89 | Line 90 | formstr(                               /* return format identifier */
90  
91   /* Add modifier to our list to track */
92   MODCONT *
93 < addmodifier(char *modn, char *outf, char *binv, int bincnt)
93 > addmodifier(char *modn, char *outf, char *prms, char *binv, int bincnt)
94   {
95          LUENT   *lep = lu_find(&modconttab,modn);
96          MODCONT *mp;
# Line 126 | Line 127 | addmodifier(char *modn, char *outf, char *binv, int bi
127                  error(SYSTEM, "out of memory in addmodifier");
128          mp->outspec = outf;             /* XXX assumes static string */
129          mp->modname = modn;             /* XXX assumes static string */
130 +        mp->params = prms;
131          mp->binv = ebinv;
132          mp->nbins = bincnt;
133          memset(mp->cbin, 0, sizeof(DCOLOR)*bincnt);
# Line 137 | Line 139 | addmodifier(char *modn, char *outf, char *binv, int bi
139   }
140  
141  
142 < /* add modifiers from a file list */
142 > /* Add modifiers from a file list */
143   void
144 < addmodfile(char *fname, char *outf, char *binv, int bincnt)
144 > addmodfile(char *fname, char *outf, char *prms, char *binv, int bincnt)
145   {
146          char    *mname[MAXMODLIST];
147          int     i;
# Line 149 | Line 151 | addmodfile(char *fname, char *outf, char *binv, int bi
151                  error(SYSTEM, errmsg);
152          }
153          for (i = 0; mname[i]; i++)      /* add each one */
154 <                addmodifier(mname[i], outf, binv, bincnt);
154 >                addmodifier(mname[i], outf, prms, binv, bincnt);
155   }
156  
157  
# Line 159 | Line 161 | quit(                  /* quit program */
161   )
162   {
163          if (nchild > 0)         /* close children if any */
164 <                end_children();
164 >                end_children(code != 0);
165          exit(code);
166   }
167  
# Line 178 | Line 180 | rcinit()
180                                          /* set shared memory boundary */
181                  shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY");
182          }
183 +        for (i = 0; i < nsources; i++)  /* tracing to sources as well */
184 +                source[i].sflags |= SFOLLOW;
185          if (yres > 0) {                 /* set up flushing & ray counts */
186                  if (xres > 0)
187                          raysleft = (RNUMBER)xres*yres;
# Line 188 | Line 192 | rcinit()
192          if ((account = accumulate) > 1)
193                  raysleft *= accumulate;
194          waitflush = (yres > 0) & (xres > 1) ? 0 : xres;
191 –                                        /* tracing to sources as well */
192 –        for (i = 0; i < nsources; i++)
193 –                source[i].sflags |= SFOLLOW;
195  
196          if (nproc > 1 && in_rchild())   /* forked child? */
197                  return;                 /* return to main processing loop */
# Line 218 | Line 219 | static void
219   trace_contrib(RAY *r)
220   {
221          MODCONT *mp;
222 +        double  bval;
223          int     bn;
224          RREAL   contr[3];
225  
226          if (r->ro == NULL || r->ro->omod == OVOID)
227                  return;
228 +                                                /* shadow ray not on source? */
229 +        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
230 +                return;
231  
232          mp = (MODCONT *)lu_find(&modconttab,objptr(r->ro->omod)->oname)->data;
233  
234          if (mp == NULL)                         /* not in our list? */
235                  return;
236  
237 <        worldfunc(RCCONTEXT, r);                /* get bin number */
238 <        bn = (int)(evalue(mp->binv) + .5);
239 <        if ((bn < 0) | (bn >= mp->nbins)) {
240 <                error(WARNING, "bad bin number (ignored)");
237 >        worldfunc(RCCONTEXT, r);                /* else set context */
238 >        set_eparams((char *)mp->params);
239 >        if ((bval = evalue(mp->binv)) <= -.5)   /* and get bin number */
240 >                return;                         /* silently ignore negatives */
241 >        if ((bn = (int)(bval + .5)) >= mp->nbins) {
242 >                sprintf(errmsg, "bad bin number (%d ignored)", bn);
243 >                error(WARNING, errmsg);
244                  return;
245          }
246 <        raycontrib(contr, r, PRIMARY);
246 >        raycontrib(contr, r, PRIMARY);          /* compute coefficient */
247          if (contrib)
248 <                multcolor(contr, r->rcol);
248 >                multcolor(contr, r->rcol);      /* -> contribution */
249          addcolor(mp->cbin[bn], contr);
250   }
251  
# Line 254 | Line 262 | eval_irrad(FVECT org, FVECT dir)
262          thisray.rdir[2] = -dir[2];
263          thisray.rmax = 0.0;
264          rayorigin(&thisray, PRIMARY, NULL, NULL);
265 <        thisray.rot = 1e-5;             /* pretend we hit surface */
265 >                                        /* pretend we hit surface */
266 >        thisray.rt = thisray.rot = 1e-5;
267          thisray.rod = 1.0;
268 +        VCOPY(thisray.ron, dir);
269          VSUM(thisray.rop, org, dir, 1e-4);
270          samplendx++;                    /* compute result */
271          (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
# Line 313 | Line 323 | rcontrib()
323   #endif
324          while (getvec(orig) == 0 && getvec(direc) == 0) {
325                  d = normalize(direc);
326 <                if (nchild != -1 && (d == 0.0) & (accumulate != 1)) {
326 >                if (nchild != -1 && (d == 0.0) & (accumulate == 0)) {
327                          if (!ignore_warning_given++)
328                                  error(WARNING,
329                                  "dummy ray(s) ignored during accumulation\n");
# Line 323 | Line 333 | rcontrib()
333                          lastray = lastdone = 0;
334                  ++lastray;
335                  if (d == 0.0) {                         /* zero ==> flush */
336 <                        if ((yres <= 0) | (xres <= 0))
337 <                                waitflush = 1;          /* flush right after */
338 <                        account = 1;
336 >                        if ((yres <= 0) | (xres <= 1))
337 >                                waitflush = 1;          /* flush after */
338 >                        if (nchild == -1)
339 >                                account = 1;
340                  } else if (imm_irrad) {                 /* else compute */
341                          eval_irrad(orig, direc);
342                  } else {
# Line 336 | Line 347 | rcontrib()
347                  if (raysleft && !--raysleft)
348                          break;          /* preemptive EOI */
349          }
350 <        if ((accumulate <= 0) | (account < accumulate)) {
350 >        if (nchild != -1 && (accumulate <= 0) | (account < accumulate)) {
351                  if (account < accumulate) {
352                          error(WARNING, "partial accumulation in final record");
353                          accumulate -= account;
# Line 344 | Line 355 | rcontrib()
355                  account = 1;            /* output accumulated totals */
356                  done_contrib();
357          }
358 +        lu_done(&ofiletab);             /* close output files */
359          if (raysleft)
360                  error(USER, "unexpected EOF on input");
349 –        lu_done(&ofiletab);             /* close output files */
361   }

Diff Legend

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