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.26 by greg, Wed Jul 15 23:38:37 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 102 | Line 103 | addmodifier(char *modn, char *outf, char *binv, int bi
103          }
104          if (nmods >= MAXMODLIST)
105                  error(INTERNAL, "too many modifiers");
106 +        if (!strcmp(modn, VOIDID)) {
107 +                sprintf(errmsg, "cannot track '%s' modifier", VOIDID);
108 +                error(USER, errmsg);
109 +        }
110          modname[nmods++] = modn;        /* XXX assumes static string */
111          lep->key = modn;                /* XXX assumes static string */
112          if (binv == NULL)
# Line 126 | Line 131 | addmodifier(char *modn, char *outf, char *binv, int bi
131                  error(SYSTEM, "out of memory in addmodifier");
132          mp->outspec = outf;             /* XXX assumes static string */
133          mp->modname = modn;             /* XXX assumes static string */
134 +        mp->params = prms;
135          mp->binv = ebinv;
136          mp->nbins = bincnt;
137          memset(mp->cbin, 0, sizeof(DCOLOR)*bincnt);
# Line 137 | Line 143 | addmodifier(char *modn, char *outf, char *binv, int bi
143   }
144  
145  
146 < /* add modifiers from a file list */
146 > /* Add modifiers from a file list */
147   void
148 < addmodfile(char *fname, char *outf, char *binv, int bincnt)
148 > addmodfile(char *fname, char *outf, char *prms, char *binv, int bincnt)
149   {
150          char    *mname[MAXMODLIST];
151          int     i;
# Line 149 | Line 155 | addmodfile(char *fname, char *outf, char *binv, int bi
155                  error(SYSTEM, errmsg);
156          }
157          for (i = 0; mname[i]; i++)      /* add each one */
158 <                addmodifier(mname[i], outf, binv, bincnt);
158 >                addmodifier(mname[i], outf, prms, binv, bincnt);
159   }
160  
161  
# Line 159 | Line 165 | quit(                  /* quit program */
165   )
166   {
167          if (nchild > 0)         /* close children if any */
168 <                end_children();
168 >                end_children(code != 0);
169          exit(code);
170   }
171  
# Line 178 | Line 184 | rcinit()
184                                          /* set shared memory boundary */
185                  shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY");
186          }
187 +        for (i = 0; i < nsources; i++)  /* tracing to sources as well */
188 +                source[i].sflags |= SFOLLOW;
189          if (yres > 0) {                 /* set up flushing & ray counts */
190                  if (xres > 0)
191                          raysleft = (RNUMBER)xres*yres;
# Line 188 | Line 196 | rcinit()
196          if ((account = accumulate) > 1)
197                  raysleft *= accumulate;
198          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;
199  
200          if (nproc > 1 && in_rchild())   /* forked child? */
201                  return;                 /* return to main processing loop */
# Line 218 | Line 223 | static void
223   trace_contrib(RAY *r)
224   {
225          MODCONT *mp;
226 +        double  bval;
227          int     bn;
228          RREAL   contr[3];
229  
230          if (r->ro == NULL || r->ro->omod == OVOID)
231                  return;
232 +                                                /* shadow ray not on source? */
233 +        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
234 +                return;
235  
236          mp = (MODCONT *)lu_find(&modconttab,objptr(r->ro->omod)->oname)->data;
237  
238          if (mp == NULL)                         /* not in our list? */
239                  return;
240  
241 <        worldfunc(RCCONTEXT, r);                /* get bin number */
242 <        bn = (int)(evalue(mp->binv) + .5);
243 <        if ((bn < 0) | (bn >= mp->nbins)) {
244 <                error(WARNING, "bad bin number (ignored)");
241 >        worldfunc(RCCONTEXT, r);                /* else set context */
242 >        set_eparams((char *)mp->params);
243 >        if ((bval = evalue(mp->binv)) <= -.5)   /* and get bin number */
244 >                return;                         /* silently ignore negatives */
245 >        if ((bn = (int)(bval + .5)) >= mp->nbins) {
246 >                sprintf(errmsg, "bad bin number (%d ignored)", bn);
247 >                error(WARNING, errmsg);
248                  return;
249          }
250 <        raycontrib(contr, r, PRIMARY);
250 >        raycontrib(contr, r, PRIMARY);          /* compute coefficient */
251          if (contrib)
252 <                multcolor(contr, r->rcol);
252 >                multcolor(contr, r->rcol);      /* -> contribution */
253          addcolor(mp->cbin[bn], contr);
254   }
255  
# Line 254 | Line 266 | eval_irrad(FVECT org, FVECT dir)
266          thisray.rdir[2] = -dir[2];
267          thisray.rmax = 0.0;
268          rayorigin(&thisray, PRIMARY, NULL, NULL);
269 <        thisray.rot = 1e-5;             /* pretend we hit surface */
269 >                                        /* pretend we hit surface */
270 >        thisray.rt = thisray.rot = 1e-5;
271          thisray.rod = 1.0;
272 +        VCOPY(thisray.ron, dir);
273          VSUM(thisray.rop, org, dir, 1e-4);
274          samplendx++;                    /* compute result */
275          (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
# Line 313 | Line 327 | rcontrib()
327   #endif
328          while (getvec(orig) == 0 && getvec(direc) == 0) {
329                  d = normalize(direc);
330 <                if (nchild != -1 && (d == 0.0) & (accumulate != 1)) {
330 >                if (nchild != -1 && (d == 0.0) & (accumulate == 0)) {
331                          if (!ignore_warning_given++)
332                                  error(WARNING,
333                                  "dummy ray(s) ignored during accumulation\n");
# Line 323 | Line 337 | rcontrib()
337                          lastray = lastdone = 0;
338                  ++lastray;
339                  if (d == 0.0) {                         /* zero ==> flush */
340 <                        if ((yres <= 0) | (xres <= 0))
341 <                                waitflush = 1;          /* flush right after */
342 <                        account = 1;
340 >                        if ((yres <= 0) | (xres <= 1))
341 >                                waitflush = 1;          /* flush after */
342 >                        if (nchild == -1)
343 >                                account = 1;
344                  } else if (imm_irrad) {                 /* else compute */
345                          eval_irrad(orig, direc);
346                  } else {
# Line 336 | Line 351 | rcontrib()
351                  if (raysleft && !--raysleft)
352                          break;          /* preemptive EOI */
353          }
354 <        if ((accumulate <= 0) | (account < accumulate)) {
354 >        if (nchild != -1 && (accumulate <= 0) | (account < accumulate)) {
355                  if (account < accumulate) {
356                          error(WARNING, "partial accumulation in final record");
357                          accumulate -= account;
# Line 344 | Line 359 | rcontrib()
359                  account = 1;            /* output accumulated totals */
360                  done_contrib();
361          }
362 +        lu_done(&ofiletab);             /* close output files */
363          if (raysleft)
364                  error(USER, "unexpected EOF on input");
349 –        lu_done(&ofiletab);             /* close output files */
365   }

Diff Legend

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