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.23 by greg, Mon Aug 25 13:18:30 2014 UTC vs.
Revision 2.36 by greg, Mon Jul 20 15:54:29 2020 UTC

# Line 8 | Line 8 | static const char RCSid[] = "$Id$";
8  
9   #include "copyright.h"
10  
11 – #include <ctype.h>
11   #include "rcontrib.h"
12   #include "otypes.h"
13   #include "source.h"
# Line 22 | Line 21 | int    dimlist[MAXDIM];                /* sampling dimensions */
21   int     ndims = 0;                      /* number of sampling dimensions */
22   int     samplendx = 0;                  /* index for this sample */
23  
24 < static void     trace_contrib(RAY *r);  /* our trace callback */
26 < void    (*trace)() = trace_contrib;
24 > void    (*trace)() = NULL;              /* trace call (NULL before rcinit) */
25  
26   int     do_irrad = 0;                   /* compute irradiance? */
27  
28   int     rand_samp = 1;                  /* pure Monte Carlo sampling? */
29  
30   double  dstrsrc = 0.9;                  /* square source distribution */
31 < double  shadthresh = .03;               /* shadow threshold */
31 > double  shadthresh = 0.;                /* shadow threshold */
32   double  shadcert = .75;                 /* shadow certainty */
33   int     directrelay = 3;                /* number of source relays */
34   int     vspretest = 512;                /* virtual source pretest density */
# Line 42 | Line 40 | COLOR  salbedo = BLKCOLOR;             /* global scattering albedo
40   double  seccg = 0.;                     /* global scattering eccentricity */
41   double  ssampdist = 0.;                 /* scatter sampling distance */
42  
43 < double  specthresh = .15;               /* specular sampling threshold */
43 > double  specthresh = .02;               /* specular sampling threshold */
44   double  specjitter = 1.;                /* specular sampling jitter */
45  
46   int     backvis = 1;                    /* back face visibility */
# Line 53 | Line 51 | double minweight = 2e-3;               /* minimum ray weight */
51   char    *ambfile = NULL;                /* ambient file name */
52   COLOR   ambval = BLKCOLOR;              /* ambient value */
53   int     ambvwt = 0;                     /* initial weight for ambient value */
54 < double  ambacc = 0;                     /* ambient accuracy */
54 > double  ambacc = 0.;                    /* ambient accuracy */
55   int     ambres = 256;                   /* ambient resolution */
56   int     ambdiv = 350;                   /* ambient divisions */
57   int     ambssamp = 0;                   /* ambient super-samples */
# Line 68 | Line 66 | long   waitflush;                      /* how long until next flush */
66   RNUMBER lastray = 0;                    /* last ray number sent */
67   RNUMBER lastdone = 0;                   /* last ray output */
68  
69 + static void     trace_contrib(RAY *r);  /* our trace callback */
70 +
71   static void mcfree(void *p) { epfree((*(MODCONT *)p).binv); free(p); }
72  
73   LUTAB   modconttab = LU_SINIT(NULL,mcfree);     /* modifier lookup table */
# Line 102 | Line 102 | addmodifier(char *modn, char *outf, char *prms, char *
102                  sprintf(errmsg, "duplicate modifier '%s'", modn);
103                  error(USER, errmsg);
104          }
105 <        if (nmods >= MAXMODLIST)
106 <                error(INTERNAL, "too many modifiers");
105 >        if (nmods >= MAXMODLIST) {
106 >                sprintf(errmsg, "too many modifiers (%d limit)", MAXMODLIST);
107 >                error(INTERNAL, errmsg);
108 >        }
109 >        if (!strcmp(modn, VOIDID)) {
110 >                sprintf(errmsg, "cannot track '%s' modifier", VOIDID);
111 >                error(USER, errmsg);
112 >        }
113          modname[nmods++] = modn;        /* XXX assumes static string */
114          lep->key = modn;                /* XXX assumes static string */
115          if (binv == NULL)
# Line 128 | Line 134 | addmodifier(char *modn, char *outf, char *prms, char *
134                  error(SYSTEM, "out of memory in addmodifier");
135          mp->outspec = outf;             /* XXX assumes static string */
136          mp->modname = modn;             /* XXX assumes static string */
137 <        mp->params = prms;
137 >        mp->params = prms;              /* XXX assumes static string */
138          mp->binv = ebinv;
139 +        mp->bin0 = 0;
140          mp->nbins = bincnt;
141          memset(mp->cbin, 0, sizeof(DCOLOR)*bincnt);
142 <                                        /* allocate output streams */
143 <        for (i = bincnt; i-- > 0; )
144 <                getostream(mp->outspec, mp->modname, i, 1);
142 >                                        /* figure out starting bin */
143 >        while (!getostream(mp->outspec, mp->modname, mp->bin0, 1))
144 >                mp->bin0++;
145 >                                        /* allocate other output streams */
146 >        for (i = 0; ++i < mp->nbins; )
147 >                getostream(mp->outspec, mp->modname, mp->bin0+i, 1);
148          lep->data = (char *)mp;
149          return(mp);
150   }
# Line 147 | Line 157 | addmodfile(char *fname, char *outf, char *prms, char *
157          char    *mname[MAXMODLIST];
158          int     i;
159                                          /* find the file & store strings */
160 <        if (wordfile(mname, getpath(fname, getrlibpath(), R_OK)) < 0) {
160 >        i = wordfile(mname, MAXMODLIST, getpath(fname, getrlibpath(), R_OK));
161 >        if (i < 0) {
162                  sprintf(errmsg, "cannot find modifier file '%s'", fname);
163                  error(SYSTEM, errmsg);
164          }
165 +        if (i >= MAXMODLIST-1) {
166 +                sprintf(errmsg, "too many modifiers (%d limit) in file '%s'",
167 +                                MAXMODLIST-1, fname);
168 +                error(INTERNAL, errmsg);
169 +        }
170          for (i = 0; mname[i]; i++)      /* add each one */
171                  addmodifier(mname[i], outf, prms, binv, bincnt);
172   }
# Line 163 | Line 179 | quit(                  /* quit program */
179   {
180          if (nchild > 0)         /* close children if any */
181                  end_children(code != 0);
182 +        else if (nchild < 0)
183 +                _exit(code);    /* avoid flush() in child */
184          exit(code);
185   }
186  
187  
188   /* Initialize our process(es) */
189   static void
190 < rcinit()
190 > rcinit(void)
191   {
192          int     i;
193  
# Line 181 | Line 199 | rcinit()
199                                          /* set shared memory boundary */
200                  shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY");
201          }
202 +        trace = trace_contrib;          /* set up trace call-back */
203          for (i = 0; i < nsources; i++)  /* tracing to sources as well */
204                  source[i].sflags |= SFOLLOW;
205          if (yres > 0) {                 /* set up flushing & ray counts */
# Line 215 | Line 234 | rcinit()
234  
235   /************************** MAIN CALCULATION PROCESS ***********************/
236  
218 – /* Set parameters for current bin evaluation */
219 – void
220 – set_eparams(char *prms)
221 – {
222 –        static char     *last_params = NULL;
223 –        char            vname[RMAXWORD];
224 –        double          value;
225 –        char            *cpd;
226 –                                        /* check if already set */
227 –        if ((prms == NULL) | (prms == last_params))
228 –                return;
229 –        if (last_params != NULL && !strcmp(prms, last_params))
230 –                return;
231 –        last_params = prms;             /* assign each variable */
232 –        while (*prms) {
233 –                if (isspace(*prms)) {
234 –                        ++prms; continue;
235 –                }
236 –                if (!isalpha(*prms))
237 –                        goto bad_params;
238 –                cpd = vname;
239 –                while (*prms && (*prms != '=') & !isspace(*prms)) {
240 –                        if (!isid(*prms))
241 –                                goto bad_params;
242 –                        *cpd++ = *prms++;
243 –                }
244 –                if (cpd == vname)
245 –                        goto bad_params;
246 –                *cpd = '\0';
247 –                while (isspace(*prms)) prms++;
248 –                if (*prms++ != '=')
249 –                        goto bad_params;
250 –                value = atof(prms);
251 –                if ((prms = fskip(prms)) == NULL)
252 –                        goto bad_params;
253 –                while (isspace(*prms)) prms++;
254 –                prms += (*prms == ',') | (*prms == ';');
255 –                varset(vname, '=', value);
256 –        }
257 –        return;
258 – bad_params:
259 –        sprintf(errmsg, "bad parameter list '%s'", last_params);
260 –        error(USER, errmsg);
261 – }
262 –
263 –
237   /* Our trace call to sum contributions */
238   static void
239   trace_contrib(RAY *r)
# Line 310 | Line 283 | eval_irrad(FVECT org, FVECT dir)
283          thisray.rmax = 0.0;
284          rayorigin(&thisray, PRIMARY, NULL, NULL);
285                                          /* pretend we hit surface */
286 <        thisray.rt = thisray.rot = 1e-5;
286 >        thisray.rxt = thisray.rot = 1e-5;
287          thisray.rod = 1.0;
288          VCOPY(thisray.ron, dir);
289          VSUM(thisray.rop, org, dir, 1e-4);
# Line 336 | Line 309 | eval_rad(FVECT org, FVECT dir, double dmax)
309  
310   /* Accumulate and/or output ray contributions (child or only process) */
311   static void
312 < done_contrib()
312 > done_contrib(void)
313   {
314          MODCONT *mp;
315          int     i;
# Line 357 | Line 330 | done_contrib()
330  
331   /* Principal calculation loop (called by main) */
332   void
333 < rcontrib()
333 > rcontrib(void)
334   {
335          static int      ignore_warning_given = 0;
336          FVECT           orig, direc;

Diff Legend

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