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

Comparing ray/src/rt/rpmain.c (file contents):
Revision 2.23 by greg, Thu Oct 20 17:18:33 2022 UTC vs.
Revision 2.36 by greg, Fri Nov 15 20:47:42 2024 UTC

# Line 30 | Line 30 | static const char      RCSid[] = "$Id$";
30   char  *progname;                        /* argv[0] */
31   char  *octname;                         /* octree name */
32   char  *sigerr[NSIG];                    /* signal error messages */
33 – char  *shm_boundary = NULL;             /* boundary of shared memory */
33   char  *errfile = NULL;                  /* error output file */
34 <
34 > #ifdef SSKIPOPT
35 > int  *ssndx_map = NULL;                 /* loaded source index map */
36 > float  *scorr_map = NULL;               /* loaded source correction map */
37 > #endif
38   extern time_t  time();
39   extern time_t  tstart;                  /* start time */
40  
# Line 52 | Line 54 | extern double  mblur;                  /* motion blur parameter */
54  
55   extern double  dblur;                   /* depth-of-field blur parameter */
56  
57 + RGBPRIMP  out_prims = stdprims;         /* output color primitives */
58 + static RGBPRIMS  our_prims;             /* private output color primitives */
59 +
60   static void onsig(int signo);
61   static void sigdie(int  signo, char  *msg);
62   static void printdefaults(void);
63                                          /* rpict additional features */
64   #ifdef PERSIST
65   #define RPICT_FEATURES  "Persist\nParallelPersist\n" \
66 +                "ParticipatingMedia=Mist\n" \
67                  "Recovery\nIrradianceCalc\nViewTypes=v,l,a,h,s,c\n" \
68                  "HessianAmbientCache\nAmbientAveraging\nAmbientValueSharing\n" \
69                  "PixelJitter\nPixelSampling\nPixelMotion\nPixelDepthOfField\n" \
70                  "SmallSourceDrawing\nViewSequence\nProgressReporting\n" \
71 <                "AdaptiveShadowTesting\nOutputs=v,l\n"
71 >                "AdaptiveShadowTesting\nOutputs=v,l\n" \
72 >                "OutputCS=RGB,XYZ,prims\n"
73   #else
74   #define RPICT_FEATURES  "Recovery\nIrradianceCalc\nViewTypes=v,l,a,h,s,c\n" \
75 +                "ParticipatingMedia=Mist\n" \
76                  "HessianAmbientCache\nAmbientAveraging\nAmbientValueSharing\n" \
77                  "PixelJitter\nPixelSampling\nPixelMotion\nPixelDepthOfField\n" \
78                  "SmallSourceDrawing\nViewSequence\nProgressReporting\n" \
79 <                "AdaptiveShadowTesting\nOutputs=v,l\n"
79 >                "AdaptiveShadowTesting\nOutputs=v,l\n" \
80 >                "OutputCS=RGB,XYZ,prims\n"
81   #endif
82  
83  
# Line 89 | Line 98 | main(int  argc, char  *argv[])
98          char  *recover = NULL;
99          char  *outfile = NULL;
100          char  *zfile = NULL;
101 + #ifdef SSKIPOPT
102 +        char  *sskipspec = NULL;
103 +        char  *scorrfile = NULL;
104 + #endif
105          int  loadflags = ~IO_FILES;
106          int  seqstart = 0;
107          int  persist = 0;
# Line 101 | Line 114 | main(int  argc, char  *argv[])
114          progname = argv[0] = fixargv0(argv[0]);
115                                          /* feature check only? */
116          strcat(RFeatureList, RPICT_FEATURES);
117 <        if (!strcmp(argv[1], "-features"))
117 >        if (argc > 1 && !strcmp(argv[1], "-features"))
118                  return feature_status(argc-2, argv+2);
119                                          /* option city */
120          for (i = 1; i < argc; i++) {
# Line 178 | Line 191 | main(int  argc, char  *argv[])
191                                  check(3,"f");
192                                  dblur = atof(argv[++i]);
193                                  break;
194 +                        case 'R':                               /* standard RGB output */
195 +                                if (strcmp(argv[i]+2, "RGB"))
196 +                                        goto badopt;
197 +                                out_prims = stdprims;
198 +                                break;
199 +                        case 'X':                               /* XYZ output */
200 +                                if (strcmp(argv[i]+2, "XYZ"))
201 +                                        goto badopt;
202 +                                out_prims = xyzprims;
203 +                                break;
204 +                        case 'c': {                             /* chromaticities */
205 +                                int     j;
206 +                                check(3,"ffffffff");
207 +                                rval = 0;
208 +                                for (j = 0; j < 8; j++) {
209 +                                        our_prims[0][j] = atof(argv[++i]);
210 +                                        rval |= fabs(our_prims[0][j]-stdprims[0][j]) > .001;
211 +                                }
212 +                                if (rval) {
213 +                                        if (!colorprimsOK(our_prims))
214 +                                                error(USER, "illegal primary chromaticities");
215 +                                        out_prims = our_prims;
216 +                                } else
217 +                                        out_prims = stdprims;
218 +                                } break;
219                          default:
220                                  goto badopt;
221                          }
# Line 210 | Line 248 | main(int  argc, char  *argv[])
248                                  check(2,"s");
249                          recover = argv[++i];
250                          break;
251 <                case 't':                               /* timer */
252 <                        check(2,"i");
253 <                        ralrm = atoi(argv[++i]);
251 > #ifdef SSKIPOPT
252 >                case 'd':
253 >                        switch (argv[i][2]) {
254 >                        case 'S':                       /* direct skip */
255 >                                check(3,"s");
256 >                                sskipspec = argv[++i];
257 >                                break;
258 >                        case 'C':                       /* direct correction */
259 >                                check(3,"s");
260 >                                scorrfile = argv[++i];
261 >                                break;
262 >                        default:
263 >                                goto badopt;
264 >                        }
265                          break;
266 + #endif
267   #ifdef  PERSIST
268                  case 'P':                               /* persist file */
269                          if (argv[i][2] == 'P') {
# Line 226 | Line 276 | main(int  argc, char  *argv[])
276                          persistfile(argv[++i]);
277                          break;
278   #endif
279 +                case 't':                               /* timer */
280 +                        check(2,"i");
281 +                        ralrm = atoi(argv[++i]);
282 +                        break;
283                  case 'w':                               /* warnings */
284                          rval = erract[WARNING].pf != NULL;
285                          check_bool(2,rval);
# Line 240 | Line 294 | main(int  argc, char  *argv[])
294                          goto badopt;
295                  }
296          }
297 +                                        /* set/check spectral sampling */
298 +        if (setspectrsamp(CNDX, WLPART) <= 0)
299 +                error(USER, "unsupported spectral sampling");
300 +
301          err = setview(&ourview);        /* set viewing parameters */
302          if (err != NULL)
303                  error(USER, err);
# Line 319 | Line 377 | main(int  argc, char  *argv[])
377          ray_init_pmap();     /* PMAP: set up & load photon maps */
378  
379          marksources();                  /* find and mark sources */
380 <
380 > #ifdef SSKIPOPT
381 >        if (sskipspec != NULL) {        /* load source skip maps? */
382 >                if ((seqstart > 0) | (persist > 0))
383 >                        error(USER,
384 >                        "source skip map incompatible with -S/-P/-PP");
385 >                if (srcskip_open(sskipspec, scorrfile) <= 0)
386 >                        quit(1);
387 >                if ((ssndx_map = srcskip_ndxmap()) == NULL)
388 >                        quit(1);
389 >                if (scorrfile != NULL &&
390 >                                (scorr_map = srcskip_corrmap()) == NULL)
391 >                        quit(1);
392 >                srcskip_close();        /* leaves maps loaded */
393 >        }
394 > #endif
395          setambient();                   /* initialize ambient calculation */
396          
397          fflush(stdout);                 /* in case we're duplicating header */
# Line 331 | Line 403 | main(int  argc, char  *argv[])
403                          close(duped1);
404                  }
405                  if (persist == PARALLEL) {      /* multiprocessing */
406 <                        preload_objs();         /* preload scene */
335 <                        shm_boundary = (char *)malloc(16);
336 <                        strcpy(shm_boundary, "SHM_BOUNDARY");
406 >                        cow_memshare();         /* preloads scene */
407                          while ((rval=fork()) == 0) {    /* keep on forkin' */
408                                  pflock(1);
409                                  pfhold();
# Line 392 | Line 462 | badopt:
462   #undef  check_bool
463   }
464  
465 + #ifdef SSKIPOPT
466 + void                            /* source skip parameters for primary */
467 + sskip_ray(RAY *r, double h, double v)
468 + {
469 +        size_t  ssi;
470  
471 +        if (ssndx_map == NULL)
472 +                return;         /* nothing to do */
473 +
474 +                                /* rows counted from the top! */
475 +        ssi = (size_t)((1.-v)*sskip_dim[1])*sskip_dim[0] +
476 +                        (size_t)(h*sskip_dim[0]);
477 +
478 +        if (ssi >= sskip_dim[0]*sskip_dim[1])
479 +                error(CONSISTENCY, "bad index in sskip_ray()");
480 +
481 +        r->rsrc = ssndx_map[ssi];
482 +        if (scorr_map != NULL)
483 +                r->scorr = scorr_map[ssi];
484 + }
485 + #endif
486 +
487   void
488   wputs(                          /* warning output function */
489 <        char    *s
489 >        const char      *s
490   )
491   {
492          int  lasterrno = errno;
493 +        if (erract[WARNING].pf == NULL)
494 +                return;         /* called by calcomp or someone */
495          eputs(s);
496          errno = lasterrno;
497   }
# Line 406 | Line 499 | wputs(                         /* warning output function */
499  
500   void
501   eputs(                          /* put string to stderr */
502 <        char  *s
502 >        const char  *s
503   )
504   {
505          static int  midline = 0;
# Line 483 | Line 576 | printdefaults(void)                    /* print default values to stdou
576          printf("-vl %f\t\t\t# view lift\n", ourview.voff);
577          printf("-x  %-9d\t\t\t# x resolution\n", hresolu);
578          printf("-y  %-9d\t\t\t# y resolution\n", vresolu);
579 +        if (out_prims == stdprims)
580 +                printf("-pRGB\t\t\t\t# standard RGB color output\n");
581 +        else if (out_prims == xyzprims)
582 +                printf("-pXYZ\t\t\t\t# CIE XYZ color output\n");
583 +        else if (out_prims != NULL)
584 +                printf("-pc %.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f\t# output color primaries and white point\n",
585 +                                out_prims[RED][0], out_prims[RED][1],
586 +                                out_prims[GRN][0], out_prims[GRN][1],
587 +                                out_prims[BLU][0], out_prims[BLU][1],
588 +                                out_prims[WHT][0], out_prims[WHT][1]);
589          printf("-pa %f\t\t\t# pixel aspect ratio\n", pixaspect);
590          printf("-pj %f\t\t\t# pixel jitter\n", dstrpix);
591          printf("-pm %f\t\t\t# pixel motion\n", mblur);

Diff Legend

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