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.30 by greg, Tue Apr 11 17:52:40 2017 UTC vs.
Revision 2.40 by greg, Tue Feb 16 20:06:06 2021 UTC

# Line 21 | 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 */
25 < 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 52 | 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 67 | 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 101 | Line 102 | addmodifier(char *modn, char *outf, char *prms, char *
102                  sprintf(errmsg, "duplicate modifier '%s'", modn);
103                  error(USER, errmsg);
104          }
104 –        if (nmods >= MAXMODLIST) {
105 –                sprintf(errmsg, "too many modifiers (%d limit)", MAXMODLIST);
106 –                error(INTERNAL, errmsg);
107 –        }
105          if (!strcmp(modn, VOIDID)) {
106                  sprintf(errmsg, "cannot track '%s' modifier", VOIDID);
107                  error(USER, errmsg);
108          }
109 +        if (nmods >= modasiz) {         /* need bigger modifier array */
110 +                modasiz += modasiz/3 + 64;
111 +                if (modname == NULL)
112 +                        modname = (char **)malloc(modasiz*sizeof(char *));
113 +                else
114 +                        modname = (char **)realloc(modname, modasiz*sizeof(char *));
115 +                if (modname == NULL)
116 +                        error(SYSTEM, "Out of memory in addmodifier()");
117 +        }
118          modname[nmods++] = modn;        /* XXX assumes static string */
119          lep->key = modn;                /* XXX assumes static string */
120          if (binv == NULL)
# Line 153 | Line 159 | addmodifier(char *modn, char *outf, char *prms, char *
159   void
160   addmodfile(char *fname, char *outf, char *prms, char *binv, int bincnt)
161   {
162 <        char    *mname[MAXMODLIST];
163 <        int     i;
164 <                                        /* find the file & store strings */
165 <        i = wordfile(mname, MAXMODLIST, getpath(fname, getrlibpath(), R_OK));
166 <        if (i < 0) {
167 <                sprintf(errmsg, "cannot find modifier file '%s'", fname);
162 >        char    *path = getpath(fname, getrlibpath(), R_OK);
163 >        char    mod[MAXSTR];
164 >        FILE    *fp;
165 >
166 >        if (path == NULL || (fp = fopen(path, "r")) == NULL) {
167 >                sprintf(errmsg, "cannot load modifier file '%s'",
168 >                                path != NULL ? path : fname);
169                  error(SYSTEM, errmsg);
170          }
171 <        if (i >= MAXMODLIST-1) {
172 <                sprintf(errmsg, "too many modifiers (%d limit) in file '%s'",
173 <                                MAXMODLIST-1, fname);
174 <                error(INTERNAL, errmsg);
171 >        while (fgetword(mod, sizeof(mod), fp) != NULL)
172 >                addmodifier(mod, outf, prms, binv, bincnt);
173 >        fclose(fp);
174 > }
175 >
176 >
177 > /* Check if we have any more rays left (and report progress) */
178 > int
179 > morays(void)
180 > {
181 >        static RNUMBER  total_rays;
182 >        static time_t   tstart, last_report;
183 >        time_t          tnow;
184 >
185 >        if (!raysleft)
186 >                return(1);      /* unknown total, so nothing to do or say */
187 >
188 >        if (report_intvl > 0 && (tnow = time(0)) >= last_report+report_intvl) {
189 >                if (!total_rays) {
190 >                        total_rays = raysleft;
191 >                        tstart = tnow;
192 >                } else {
193 >                        sprintf(errmsg, "%.2f%% done after %.3f hours\n",
194 >                                        100.-100.*raysleft/total_rays,
195 >                                        (1./3600.)*(tnow - tstart));
196 >                        eputs(errmsg);
197 >                }
198 >                last_report = tnow;
199          }
200 <        for (i = 0; mname[i]; i++)      /* add each one */
170 <                addmodifier(mname[i], outf, prms, binv, bincnt);
200 >        return(--raysleft);
201   }
202  
203  
204 + /* Quit program */
205   void
206 < quit(                   /* quit program */
206 > quit(
207          int  code
208   )
209   {
210          if (nchild > 0)         /* close children if any */
211                  end_children(code != 0);
212 +        else if (nchild < 0)
213 +                _exit(code);    /* avoid flush() in child */
214          exit(code);
215   }
216  
# Line 196 | Line 229 | rcinit(void)
229                                          /* set shared memory boundary */
230                  shm_boundary = strcpy((char *)malloc(16), "SHM_BOUNDARY");
231          }
232 +        trace = trace_contrib;          /* set up trace call-back */
233          for (i = 0; i < nsources; i++)  /* tracing to sources as well */
234                  source[i].sflags |= SFOLLOW;
235          if (yres > 0) {                 /* set up flushing & ray counts */
# Line 279 | Line 313 | eval_irrad(FVECT org, FVECT dir)
313          thisray.rmax = 0.0;
314          rayorigin(&thisray, PRIMARY, NULL, NULL);
315                                          /* pretend we hit surface */
316 <        thisray.rt = thisray.rot = 1e-5;
316 >        thisray.rxt = thisray.rot = 1e-5;
317          thisray.rod = 1.0;
318          VCOPY(thisray.ron, dir);
319          VSUM(thisray.rop, org, dir, 1e-4);
# Line 360 | Line 394 | rcontrib(void)
394                  }
395                  done_contrib();         /* accumulate/output */
396                  ++lastdone;
397 <                if (raysleft && !--raysleft)
397 >                if (!morays())
398                          break;          /* preemptive EOI */
399          }
400          if (nchild != -1 && (accumulate <= 0) | (account < accumulate)) {

Diff Legend

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