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