ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rcmain.c
(Generate patch)

Comparing ray/src/rt/rcmain.c (file contents):
Revision 2.12 by greg, Sat Jul 19 18:19:33 2014 UTC vs.
Revision 2.32 by greg, Wed Nov 15 18:02:53 2023 UTC

# Line 13 | Line 13 | static const char      RCSid[] = "$Id$";
13   #include "random.h"
14   #include "source.h"
15   #include "ambient.h"
16 + #include "pmapray.h"
17 + #include "pmapcontrib.h"
18  
19   int     gargc;                          /* global argc */
20   char    **gargv;                        /* global argv */
# Line 41 | Line 43 | int    using_stdout = 0;               /* are we using stdout? */
43   int     imm_irrad = 0;                  /* compute immediate irradiance? */
44   int     lim_dist = 0;                   /* limit distance? */
45  
46 < const char      *modname[MAXMODLIST];   /* ordered modifier name list */
45 < int             nmods = 0;              /* number of modifiers */
46 > int     report_intvl = 0;               /* reporting interval (seconds) */
47  
48 + char    **modname = NULL;               /* ordered modifier name list */
49 + int     nmods = 0;                      /* number of modifiers */
50 + int     modasiz = 0;                    /* allocated modifier array size */
51 +
52   void    (*addobjnotify[8])() = {ambnotify, NULL};
53  
54 < char    RCCONTEXT[] = "RC";             /* our special evaluation context */
54 > char    RCCONTEXT[] = "RC.";            /* our special evaluation context */
55  
56 + #if defined(_WIN32) || defined(_WIN64)
57 + #define RCONTRIB_FEATURES       "Accumulation\nSummation\nRecovery\n" \
58 +                                "Hyperspectral\nImmediateIrradiance\n" \
59 +                                "ProgressReporting\nDistanceLimiting\n" \
60 +                                "InputFormats=a,f,d\nOutputFormats=a,f,d,c\n" \
61 +                                "Outputs=V,W\n"
62 + #else
63 + #define RCONTRIB_FEATURES       "Multiprocessing\n" \
64 +                                "Accumulation\nSummation\nRecovery\n" \
65 +                                "Hyperspectral\nImmediateIrradiance\n" \
66 +                                "ProgressReporting\nDistanceLimiting\n" \
67 +                                "InputFormats=a,f,d\nOutputFormats=a,f,d,c\n" \
68 +                                "Outputs=V,W\n"
69 + #endif
70  
71   static void
72   printdefaults(void)                     /* print default values to stdout */
73   {
74          printf("-c %-5d\t\t\t# accumulated rays per record\n", accumulate);
75 +        if (NCSAMP > 3) {
76 +                printf("-cs %-2d\t\t\t\t# number of spectral bins\n", NCSAMP);
77 +                printf("-cw %3.0f %3.0f\t\t\t# wavelength limits (nm)\n",
78 +                                WLPART[3], WLPART[0]);
79 +        }
80          printf("-V%c\t\t\t\t# output %s\n", contrib ? '+' : '-',
81                          contrib ? "contributions" : "coefficients");
82          if (imm_irrad)
# Line 63 | Line 87 | printdefaults(void)                    /* print default values to stdou
87          printf("-y %-9d\t\t\t# y resolution\n", yres);
88          printf(lim_dist ? "-ld+\t\t\t\t# limit distance on\n" :
89                          "-ld-\t\t\t\t# limit distance off\n");
90 <        printf("-h%c\t\t\t\t# %s header\n", header ? '+' : '-',
91 <                        header ? "output" : "no");
90 >        printf(header ? "-h+\t\t\t\t# output header\n" :
91 >                        "-h-\t\t\t\t# no header\n");
92          printf("-f%c%c\t\t\t\t# format input/output = %s/%s\n",
93                          inpfmt, outfmt, formstr(inpfmt), formstr(outfmt));
94          printf(erract[WARNING].pf != NULL ?
# Line 161 | Line 185 | main(int argc, char *argv[])
185   #define  check(ol,al)           if (argv[i][ol] || \
186                                  badarg(argc-i-1,argv+i+1,al)) \
187                                  goto badopt
188 < #define  bool(olen,var)         switch (argv[i][olen]) { \
188 > #define  check_bool(olen,var)           switch (argv[i][olen]) { \
189                                  case '\0': var = !var; break; \
190                                  case 'y': case 'Y': case 't': case 'T': \
191                                  case '+': case '1': var = 1; break; \
# Line 178 | Line 202 | main(int argc, char *argv[])
202          progname = argv[0] = fixargv0(argv[0]);
203          gargv = argv;
204          gargc = argc;
205 +                                        /* feature check only? */
206 +        strcat(RFeatureList, RCONTRIB_FEATURES);
207 +        if (argc > 1 && !strcmp(argv[1], "-features"))
208 +                return feature_status(argc-2, argv+2);
209 + #if defined(_WIN32) || defined(_WIN64)  /* increase file limit to maximum */
210 +        for (i = 8192; i > _IOB_ENTRIES; i >>= 1)
211 +                if (_setmaxstdio(i) == i)
212 +                        break;
213 + #endif
214                                          /* initialize calcomp routines early */
215          initfunc();
216 <        setcontext(RCCONTEXT);
216 >        calcontext(RCCONTEXT);
217                                          /* option city */
218          for (i = 1; i < argc; i++) {
219                                                  /* expand arguments */
# Line 215 | Line 248 | main(int argc, char *argv[])
248                                  error(USER, "bad number of processes");
249                          break;
250                  case 'V':                       /* output contributions */
251 <                        bool(2,contrib);
251 >                        check_bool(2,contrib);
252                          break;
253                  case 'x':                       /* x resolution */
254                          check(2,"i");
# Line 227 | Line 260 | main(int argc, char *argv[])
260                          break;
261                  case 'w':                       /* warnings */
262                          rval = (erract[WARNING].pf != NULL);
263 <                        bool(2,rval);
263 >                        check_bool(2,rval);
264                          if (rval) erract[WARNING].pf = wputs;
265                          else erract[WARNING].pf = NULL;
266                          break;
# Line 238 | Line 271 | main(int argc, char *argv[])
271                  case 'l':                       /* limit distance */
272                          if (argv[i][2] != 'd')
273                                  goto badopt;
274 <                        bool(3,lim_dist);
274 >                        check_bool(3,lim_dist);
275                          break;
276                  case 'I':                       /* immed. irradiance */
277 <                        bool(2,imm_irrad);
277 >                        check_bool(2,imm_irrad);
278                          break;
279                  case 'f':                       /* file or force or format */
280                          if (!argv[i][2]) {
# Line 250 | Line 283 | main(int argc, char *argv[])
283                                  break;
284                          }
285                          if (argv[i][2] == 'o') {
286 <                                bool(3,force_open);
286 >                                check_bool(3,force_open);
287                                  break;
288                          }
289                          setformat(argv[i]+2);
# Line 259 | Line 292 | main(int argc, char *argv[])
292                          check(2,"s");
293                          curout = argv[++i];
294                          break;
262 –                case 'c':                       /* input rays per output */
263 –                        check(2,"i");
264 –                        accumulate = atoi(argv[++i]);
265 –                        break;
295                  case 'r':                       /* recover output */
296 <                        bool(2,recover);
296 >                        check_bool(2,recover);
297                          break;
298                  case 'h':                       /* header output */
299 <                        bool(2,header);
299 >                        check_bool(2,header);
300                          break;
301                  case 'p':                       /* parameter setting(s) */
302                          check(2,"s");
303 <                        prms = argv[++i];
303 >                        set_eparams(prms = argv[++i]);
304                          break;
305 +                case 'c':                               /* spectral sampling or count */
306 +                        switch (argv[i][2]) {
307 + #if MAXCSAMP>3
308 +                        case 's':                       /* spectral bin count */
309 +                                check(3,"i");
310 +                                NCSAMP = atoi(argv[++i]);
311 +                                break;
312 +                        case 'w':                       /* wavelength extrema */
313 +                                check(3,"ff");
314 +                                WLPART[0] = atof(argv[++i]);
315 +                                WLPART[3] = atof(argv[++i]);
316 +                                break;
317 + #endif
318 +                        case '\0':                      /* sample count */
319 +                                check(2,"i");
320 +                                accumulate = atoi(argv[++i]);
321 +                                break;
322 +                        default:
323 +                                goto badopt;
324 +                        }
325 +                        break;
326                  case 'b':                       /* bin expression/count */
327                          if (argv[i][2] == 'n') {
328                                  check(3,"s");
# Line 290 | Line 340 | main(int argc, char *argv[])
340                          check(2,"s");
341                          addmodfile(argv[++i], curout, prms, binval, bincnt);
342                          break;
343 +                case 't':                       /* reporting interval */
344 +                        check(2,"i");
345 +                        report_intvl = atoi(argv[++i]);
346 +                        break;
347                  default:
348                          goto badopt;
349                  }
# Line 340 | Line 394 | main(int argc, char *argv[])
394          readoct(octname, ~(IO_FILES|IO_INFO), &thescene, NULL);
395          nsceneobjs = nobjects;
396  
397 +        /* PMAP: set up & load photon maps */
398 +        ray_init_pmap();    
399 +        
400          marksources();                  /* find and mark sources */
401 +        
402 +        /* PMAP: init photon map for light source contributions */
403 +        initPmapContrib(&modconttab, nmods);
404  
405          setambient();                   /* initialize ambient calculation */
406 <
406 >        
407          rcontrib();                     /* trace ray contributions (loop) */
408  
409          ambsync();                      /* flush ambient file */
410  
411 +        /* PMAP: free photon maps */
412 +        ray_done_pmap();    
413 +        
414          quit(0);        /* exit clean */
415  
416   badopt:
417          fprintf(stderr,
418 < "Usage: %s [-n nprocs][-V][-r][-e expr][-f source][-o ospec][-p p1=V1,p2=V2][-b binv][-bn N] {-m mod | -M file} [rtrace options] octree\n",
418 > "Usage: %s [-n nprocs][-V][-c count][-r][-e expr][-f source][-o ospec][-p p1=V1,p2=V2][-b binv][-bn N] {-m mod | -M file} [rtrace options] octree\n",
419                          progname);
420          sprintf(errmsg, "command line error at '%s'", argv[i]);
421          error(USER, errmsg);
422          return(1);      /* pro forma return */
423  
424   #undef  check
425 < #undef  bool
425 > #undef  check_bool
426   }
427  
428  
429   void
430   wputs(                          /* warning output function */
431 <        char    *s
431 >        const char      *s
432   )
433   {
434          int  lasterrno = errno;
# Line 376 | Line 439 | wputs(                         /* warning output function */
439  
440   void
441   eputs(                          /* put string to stderr */
442 <        char  *s
442 >        const char  *s
443   )
444   {
445          static int  midline = 0;

Diff Legend

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