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.19 by greg, Sat Jul 19 00:07:03 2014 UTC vs.
Revision 2.31 by greg, Fri Feb 23 03:21:24 2018 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"
12   #include "otypes.h"
13   #include "source.h"
# Line 27 | Line 29 | int    do_irrad = 0;                   /* compute irradiance? */
29   int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
30  
31   double  dstrsrc = 0.9;                  /* square source distribution */
32 < double  shadthresh = .03;               /* shadow threshold */
32 > double  shadthresh = 0.;                /* shadow threshold */
33   double  shadcert = .75;                 /* shadow certainty */
34   int     directrelay = 3;                /* number of source relays */
35   int     vspretest = 512;                /* virtual source pretest density */
# Line 39 | Line 41 | COLOR  salbedo = BLKCOLOR;             /* global scattering albedo
41   double  seccg = 0.;                     /* global scattering eccentricity */
42   double  ssampdist = 0.;                 /* scatter sampling distance */
43  
44 < double  specthresh = .15;               /* specular sampling threshold */
44 > double  specthresh = .02;               /* specular sampling threshold */
45   double  specjitter = 1.;                /* specular sampling jitter */
46  
47   int     backvis = 1;                    /* back face visibility */
# Line 50 | Line 52 | double minweight = 2e-3;               /* minimum ray weight */
52   char    *ambfile = NULL;                /* ambient file name */
53   COLOR   ambval = BLKCOLOR;              /* ambient value */
54   int     ambvwt = 0;                     /* initial weight for ambient value */
55 < double  ambacc = 0;                     /* ambient accuracy */
55 > double  ambacc = 0.;                    /* ambient accuracy */
56   int     ambres = 256;                   /* ambient resolution */
57   int     ambdiv = 350;                   /* ambient divisions */
58   int     ambssamp = 0;                   /* ambient super-samples */
# Line 88 | 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 99 | Line 101 | addmodifier(char *modn, char *outf, char *binv, int bi
101                  sprintf(errmsg, "duplicate modifier '%s'", modn);
102                  error(USER, errmsg);
103          }
104 <        if (nmods >= MAXMODLIST)
105 <                error(INTERNAL, "too many modifiers");
104 >        if (nmods >= MAXMODLIST) {
105 >                sprintf(errmsg, "too many modifiers (%d limit)", MAXMODLIST);
106 >                error(INTERNAL, errmsg);
107 >        }
108 >        if (!strcmp(modn, VOIDID)) {
109 >                sprintf(errmsg, "cannot track '%s' modifier", VOIDID);
110 >                error(USER, errmsg);
111 >        }
112          modname[nmods++] = modn;        /* XXX assumes static string */
113          lep->key = modn;                /* XXX assumes static string */
114          if (binv == NULL)
# Line 125 | Line 133 | addmodifier(char *modn, char *outf, char *binv, int bi
133                  error(SYSTEM, "out of memory in addmodifier");
134          mp->outspec = outf;             /* XXX assumes static string */
135          mp->modname = modn;             /* XXX assumes static string */
136 +        mp->params = prms;              /* XXX assumes static string */
137          mp->binv = ebinv;
138 +        mp->bin0 = 0;
139          mp->nbins = bincnt;
140          memset(mp->cbin, 0, sizeof(DCOLOR)*bincnt);
141 <                                        /* allocate output streams */
142 <        for (i = bincnt; i-- > 0; )
143 <                getostream(mp->outspec, mp->modname, i, 1);
141 >                                        /* figure out starting bin */
142 >        while (!getostream(mp->outspec, mp->modname, mp->bin0, 1))
143 >                mp->bin0++;
144 >                                        /* allocate other output streams */
145 >        for (i = 0; ++i < mp->nbins; )
146 >                getostream(mp->outspec, mp->modname, mp->bin0+i, 1);
147          lep->data = (char *)mp;
148          return(mp);
149   }
# Line 138 | Line 151 | addmodifier(char *modn, char *outf, char *binv, int bi
151  
152   /* Add modifiers from a file list */
153   void
154 < addmodfile(char *fname, char *outf, char *binv, int bincnt)
154 > addmodfile(char *fname, char *outf, char *prms, char *binv, int bincnt)
155   {
156          char    *mname[MAXMODLIST];
157          int     i;
158                                          /* find the file & store strings */
159 <        if (wordfile(mname, getpath(fname, getrlibpath(), R_OK)) < 0) {
159 >        i = wordfile(mname, MAXMODLIST, getpath(fname, getrlibpath(), R_OK));
160 >        if (i < 0) {
161                  sprintf(errmsg, "cannot find modifier file '%s'", fname);
162                  error(SYSTEM, errmsg);
163          }
164 +        if (i >= MAXMODLIST-1) {
165 +                sprintf(errmsg, "too many modifiers (%d limit) in file '%s'",
166 +                                MAXMODLIST-1, fname);
167 +                error(INTERNAL, errmsg);
168 +        }
169          for (i = 0; mname[i]; i++)      /* add each one */
170 <                addmodifier(mname[i], outf, binv, bincnt);
170 >                addmodifier(mname[i], outf, prms, binv, bincnt);
171   }
172  
173  
# Line 165 | Line 184 | quit(                  /* quit program */
184  
185   /* Initialize our process(es) */
186   static void
187 < rcinit()
187 > rcinit(void)
188   {
189          int     i;
190  
# Line 222 | Line 241 | trace_contrib(RAY *r)
241  
242          if (r->ro == NULL || r->ro->omod == OVOID)
243                  return;
244 +                                                /* shadow ray not on source? */
245 +        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
246 +                return;
247  
248          mp = (MODCONT *)lu_find(&modconttab,objptr(r->ro->omod)->oname)->data;
249  
250          if (mp == NULL)                         /* not in our list? */
251                  return;
230 –                                                /* shadow ray not on source? */
231 –        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
232 –                return;
252  
253 <        worldfunc(RCCONTEXT, r);                /* else get bin number */
254 <        if ((bval = evalue(mp->binv)) <= -.5)
255 <                return;                         /* silently ignore */
253 >        worldfunc(RCCONTEXT, r);                /* else set context */
254 >        set_eparams((char *)mp->params);
255 >        if ((bval = evalue(mp->binv)) <= -.5)   /* and get bin number */
256 >                return;                         /* silently ignore negatives */
257          if ((bn = (int)(bval + .5)) >= mp->nbins) {
258 <                error(WARNING, "bad bin number (ignored)");
258 >                sprintf(errmsg, "bad bin number (%d ignored)", bn);
259 >                error(WARNING, errmsg);
260                  return;
261          }
262          raycontrib(contr, r, PRIMARY);          /* compute coefficient */
# Line 284 | Line 305 | eval_rad(FVECT org, FVECT dir, double dmax)
305  
306   /* Accumulate and/or output ray contributions (child or only process) */
307   static void
308 < done_contrib()
308 > done_contrib(void)
309   {
310          MODCONT *mp;
311          int     i;
# Line 305 | Line 326 | done_contrib()
326  
327   /* Principal calculation loop (called by main) */
328   void
329 < rcontrib()
329 > rcontrib(void)
330   {
331          static int      ignore_warning_given = 0;
332          FVECT           orig, direc;

Diff Legend

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