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

Comparing ray/src/rt/rxtmain.cpp (file contents):
Revision 2.2 by greg, Tue Mar 12 16:54:51 2024 UTC vs.
Revision 2.12 by greg, Wed Mar 5 18:56:28 2025 UTC

# Line 2 | Line 2
2   static const char       RCSid[] = "$Id$";
3   #endif
4   /*
5 < *  rtmain.c - main for rtrace per-ray calculation program
5 > *  rxtmain.cpp - main for per-ray calculation program
6   */
7  
8   #include "copyright.h"
# Line 43 | Line 43 | extern void  tranotify(OBJECT obj);
43   char  *tralist[MAXMODLIST];             /* list of modifers to trace (or no) */
44   int  traincl = -1;                      /* include == 1, exclude == 0 */
45  
46 + double  (*sens_curve)(const SCOLOR scol) = NULL;        /* spectral conversion for 1-channel */
47 + double  out_scalefactor = 1;            /* output calibration scale factor */
48 + RGBPRIMP  out_prims = stdprims;         /* output color primitives (NULL if spectral) */
49 + static RGBPRIMS  our_prims;             /* private output color primitives */
50 +
51   static void onsig(int  signo);
52   static void sigdie(int  signo, const char  *msg);
53   static void printdefaults(void);
54  
55 < #define RATRACE_FEATURES        "IrradianceCalc\nIrradianceCalc\nDistanceLimiting\n" \
56 <                        "HessianAmbientCache\nAmbientAveraging\n" \
57 <                        "AmbientValueSharing\nAdaptiveShadowTesting\n" \
58 <                        "Outputs=o,d,v,V,w,W,l,L,c,p,n,N,s,m,M,r,x,R,X,~\n"
55 > #define RXTRACE_FEATURES        "IrradianceCalc\nMultiprocessing\nDistanceLimiting\n" \
56 >                                "HessianAmbientCache\nAmbientAveraging\n" \
57 >                                "AmbientValueSharing\nAdaptiveShadowTesting\n" \
58 >                                "Outputs=o,d,v,V,w,W,l,L,c,p,n,N,s,m,M,r,x,R,X,~\n" \
59 >                                "OutputCS=RGB,XYZ,Y,S,M,prims,spec\n"
60  
61   int
62   main(int  argc, char  *argv[])
# Line 69 | Line 75 | main(int  argc, char  *argv[])
75          int  rval;
76          int  i;
77                                          /* global program name */
78 <        progname = argv[0] = fixargv0(argv[0]);
78 >        progname = argv[0];
79                                          /* feature check only? */
80 <        strcat(RFeatureList, RATRACE_FEATURES);
80 >        strcat(RFeatureList, RXTRACE_FEATURES);
81          if (argc > 1 && !strcmp(argv[1], "-features"))
82                  return feature_status(argc-2, argv+2);
83                                          /* add trace notify function */
# Line 136 | Line 142 | main(int  argc, char  *argv[])
142                          break;
143                  case 'I':                               /* immed. irradiance */
144                          rval = myRTmanager.rtFlags & RTimmIrrad;
145 <                        check_bool(3,rval);
145 >                        check_bool(2,rval);
146                          if (rval) myRTmanager.rtFlags |= RTimmIrrad;
147                          else myRTmanager.rtFlags &= ~RTimmIrrad;
148                          break;
# Line 221 | Line 227 | main(int  argc, char  *argv[])
227                                  goto badopt;
228                          }
229                          break;
230 +                case 'p':                               /* value output */
231 +                        switch (argv[i][2]) {
232 +                        case 'R':                       /* standard RGB output */
233 +                                if (strcmp(argv[i]+2, "RGB"))
234 +                                        goto badopt;
235 +                                out_prims = stdprims;
236 +                                out_scalefactor = 1;
237 +                                sens_curve = NULL;
238 +                                break;
239 +                        case 'X':                       /* XYZ output */
240 +                                if (strcmp(argv[i]+2, "XYZ"))
241 +                                        goto badopt;
242 +                                out_prims = xyzprims;
243 +                                out_scalefactor = WHTEFFICACY;
244 +                                sens_curve = NULL;
245 +                                break;
246 +                        case 'c': {
247 +                                int     j;
248 +                                check(3,"ffffffff");
249 +                                rval = 0;
250 +                                for (j = 0; j < 8; j++) {
251 +                                        our_prims[0][j] = atof(argv[++i]);
252 +                                        rval |= fabs(our_prims[0][j]-stdprims[0][j]) > .001;
253 +                                }
254 +                                if (rval) {
255 +                                        if (!colorprimsOK(our_prims))
256 +                                                error(USER, "illegal primary chromaticities");
257 +                                        out_prims = our_prims;
258 +                                } else
259 +                                        out_prims = stdprims;
260 +                                out_scalefactor = 1;
261 +                                sens_curve = NULL;
262 +                                } break;
263 +                        case 'Y':                       /* photopic response */
264 +                                if (argv[i][3])
265 +                                        goto badopt;
266 +                                sens_curve = scolor_photopic;
267 +                                out_scalefactor = WHTEFFICACY;
268 +                                break;
269 +                        case 'S':                       /* scotopic response */
270 +                                if (argv[i][3])
271 +                                        goto badopt;
272 +                                sens_curve = scolor_scotopic;
273 +                                out_scalefactor = WHTSCOTOPIC;
274 +                                break;
275 +                        case 'M':                       /* melanopic response */
276 +                                if (argv[i][3])
277 +                                        goto badopt;
278 +                                sens_curve = scolor_melanopic;
279 +                                out_scalefactor = WHTMELANOPIC;
280 +                                break;
281 +                        default:
282 +                                goto badopt;
283 +                        }
284 +                        break;
285 + #if MAXCSAMP>3
286 +                case 'c':                               /* output spectral results */
287 +                        if (argv[i][2] != 'o')
288 +                                goto badopt;
289 +                        rval = (out_prims == NULL) & (sens_curve == NULL);
290 +                        check_bool(3,rval);
291 +                        if (rval) {
292 +                                out_prims = NULL;
293 +                                sens_curve = NULL;
294 +                        } else if (out_prims == NULL)
295 +                                out_prims = stdprims;
296 +                        break;
297 + #endif
298                  default:
299                          goto badopt;
300                  }
301          }
302 +                                        /* set/check spectral sampling */
303 +        rval = setspectrsamp(CNDX, WLPART);
304 +        if (rval < 0)
305 +                error(USER, "unsupported spectral sampling");
306 +        if (sens_curve != NULL)
307 +                out_prims = NULL;
308 +        else if (out_prims != NULL) {
309 +                if (!rval)
310 +                        error(WARNING, "spectral range incompatible with color output");
311 +        } else if (NCSAMP == 3)
312 +                out_prims = stdprims;   /* 3 samples do not a spectrum make */
313                                          /* set up signal handling */
314          sigdie(SIGINT, "Interrupt");
315   #ifdef SIGHUP
# Line 268 | Line 353 | main(int  argc, char  *argv[])
353          if (outform != 'a')
354                  SET_FILE_BINARY(stdout);
355          if (doheader) {                 /* print header? */
356 <                static char     fmt[] = OCTFMT;
357 <                FILE *          octfp = fopen(argv[i], "rb");
273 <                if (checkheader(octfp, fmt, stdout) < 0)
274 <                        error(USER, "bad octree header");
275 <                fclose(octfp);
356 >                newheader("RADIANCE", stdout);
357 >                fputs(myRTmanager.GetHeadStr(), stdout);
358                  printargs(i, argv, stdout);
359                  printf("SOFTWARE= %s\n", VersionID);
360                  fputnow(stdout);
361                  if (rval > 0)           /* saved from setrtoutput() call */
362 <                        printf("NCOMP=%d\n", rval);
362 >                        fputncomp(rval, stdout);
363 >                if (NCSAMP > 3)
364 >                        fputwlsplit(WLPART, stdout);
365 >                if ((out_prims != stdprims) & (out_prims != NULL))
366 >                        fputprims(out_prims, stdout);
367                  if ((outform == 'f') | (outform == 'd'))
368                          fputendian(stdout);
369                  fputformat(formstr(outform), stdout);
370 <                putchar('\n');
370 >                fputc('\n', stdout);    /* end of header */
371          }
372          rtrace(NULL, nproc);            /* trace rays */
373          quit(0);                        /* clean up & exit */
# Line 307 | Line 393 | wputs(                         /* warning output function */
393  
394   void
395   eputs(                          /* put string to stderr */
396 <        char  *s
396 >        const char  *s
397   )
398   {
399          static int  midline = 0;
# Line 398 | Line 484 | printdefaults(void)                    /* print default values to stdou
484                  case '~': printf(" tilde"); break;
485                  }
486          putchar('\n');
487 +        if (sens_curve == scolor_photopic)
488 +                printf("-pY\t\t\t\t# photopic output\n");
489 +        else if (sens_curve == scolor_scotopic)
490 +                printf("-pS\t\t\t\t# scotopic output\n");
491 +        else if (sens_curve == scolor_melanopic)
492 +                printf("-pM\t\t\t\t# melanopic output\n");
493 +        else if (out_prims == stdprims)
494 +                printf("-pRGB\t\t\t\t# standard RGB color output\n");
495 +        else if (out_prims == xyzprims)
496 +                printf("-pXYZ\t\t\t\t# CIE XYZ color output\n");
497 +        else if (out_prims != NULL)
498 +                printf("-pc %.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f\t# output color primaries and white point\n",
499 +                                out_prims[RED][0], out_prims[RED][1],
500 +                                out_prims[GRN][0], out_prims[GRN][1],
501 +                                out_prims[BLU][0], out_prims[BLU][1],
502 +                                out_prims[WHT][0], out_prims[WHT][1]);
503 +        if ((sens_curve == NULL) & (NCSAMP > 3))
504 +                printf(out_prims != NULL ? "-co-\t\t\t\t# output tristimulus colors\n" :
505 +                                "-co+\t\t\t\t# output spectral values\n");
506          printf(erract[WARNING].pf != NULL ?
507                          "-w+\t\t\t\t# warning messages on\n" :
508                          "-w-\t\t\t\t# warning messages off\n");

Diff Legend

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