--- ray/src/rt/rtmain.c 2003/02/22 02:07:29 2.1
+++ ray/src/rt/rtmain.c 2023/12/11 18:33:53 2.50
@@ -1,80 +1,28 @@
#ifndef lint
-static const char RCSid[] = "$Id: rtmain.c,v 2.1 2003/02/22 02:07:29 greg Exp $";
+static const char RCSid[] = "$Id: rtmain.c,v 2.50 2023/12/11 18:33:53 greg Exp $";
#endif
/*
* rtmain.c - main for rtrace per-ray calculation program
*/
-/* ====================================================================
- * The Radiance Software License, Version 1.0
- *
- * Copyright (c) 1990 - 2002 The Regents of the University of California,
- * through Lawrence Berkeley National Laboratory. All rights reserved.
- *
- * Redistribution and use in source and binary forms, with or without
- * modification, are permitted provided that the following conditions
- * are met:
- *
- * 1. Redistributions of source code must retain the above copyright
- * notice, this list of conditions and the following disclaimer.
- *
- * 2. Redistributions in binary form must reproduce the above copyright
- * notice, this list of conditions and the following disclaimer in
- * the documentation and/or other materials provided with the
- * distribution.
- *
- * 3. The end-user documentation included with the redistribution,
- * if any, must include the following acknowledgment:
- * "This product includes Radiance software
- * (http://radsite.lbl.gov/)
- * developed by the Lawrence Berkeley National Laboratory
- * (http://www.lbl.gov/)."
- * Alternately, this acknowledgment may appear in the software itself,
- * if and wherever such third-party acknowledgments normally appear.
- *
- * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
- * and "The Regents of the University of California" must
- * not be used to endorse or promote products derived from this
- * software without prior written permission. For written
- * permission, please contact radiance@radsite.lbl.gov.
- *
- * 5. Products derived from this software may not be called "Radiance",
- * nor may "Radiance" appear in their name, without prior written
- * permission of Lawrence Berkeley National Laboratory.
- *
- * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
- * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
- * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
- * DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
- * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
- * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
- * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
- * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
- * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
- * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
- * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
- * SUCH DAMAGE.
- * ====================================================================
- *
- * This software consists of voluntary contributions made by many
- * individuals on behalf of Lawrence Berkeley National Laboratory. For more
- * information on Lawrence Berkeley National Laboratory, please see
- * .
- */
+#include "copyright.h"
-#include "ray.h"
+#include
+#include "platform.h"
+#include "rtprocess.h" /* getpid() */
+#include "resolu.h"
+#include "ray.h"
#include "source.h"
-
#include "ambient.h"
-
#include "random.h"
-
#include "paths.h"
+#include "pmapray.h"
-#include
+extern char *progname; /* global argv[0] */
-#include
+extern char *shm_boundary; /* boundary of shared memory */
+
/* persistent processes define */
#ifdef F_SETLKW
#define PERSIST 1 /* normal persist */
@@ -82,58 +30,94 @@ static const char RCSid[] = "$Id: rtmain.c,v 2.1 2003/
#define PCHILD 3 /* child of normal persist */
#endif
-char *progname; /* argv[0] */
+char *sigerr[NSIG]; /* signal error messages */
+char *errfile = NULL; /* error output file */
-char *octname; /* octree name */
+int nproc = 1; /* number of processes */
-char *sigerr[NSIG]; /* signal error messages */
+extern int setrtoutput(void); /* set output values */
-char *shm_boundary = NULL; /* boundary of shared memory */
+int inform = 'a'; /* input format */
+int outform = 'a'; /* output format */
+char *outvals = "v"; /* output specification */
-char *errfile = NULL; /* error output file */
+int hresolu = 0; /* horizontal (scan) size */
+int vresolu = 0; /* vertical resolution */
-extern char *formstr(); /* string from format */
-extern int inform; /* input format */
-extern int outform; /* output format */
-extern char *outvals; /* output values */
+extern int castonly; /* only doing ray-casting? */
-extern int hresolu; /* horizontal resolution */
-extern int vresolu; /* vertical resolution */
+int imm_irrad = 0; /* compute immediate irradiance? */
+int lim_dist = 0; /* limit distance? */
-extern int imm_irrad; /* compute immediate irradiance? */
-extern int lim_dist; /* limit distance? */
+#ifndef MAXMODLIST
+#define MAXMODLIST 1024 /* maximum modifiers we'll track */
+#endif
-extern char *tralist[]; /* list of modifers to trace (or no) */
-extern int traincl; /* include == 1, exclude == 0 */
+extern void (*addobjnotify[])(); /* object notification calls */
+extern void tranotify(OBJECT obj);
-void onsig();
-void sigdie();
-void printdefaults();
+char *tralist[MAXMODLIST]; /* list of modifers to trace (or no) */
+int traincl = -1; /* include == 1, exclude == 0 */
+double (*sens_curve)(SCOLOR scol) = NULL; /* spectral conversion for 1-channel */
+double out_scalefactor = 1; /* output calibration scale factor */
+RGBPRIMP out_prims = stdprims; /* output color primitives (NULL if spectral) */
+static RGBPRIMS our_prims; /* private output color primitives */
+static int loadflags = ~IO_FILES; /* what to load from octree */
+
+static void onsig(int signo);
+static void sigdie(int signo, char *msg);
+static void printdefaults(void);
+
+#ifdef PERSIST
+#define RTRACE_FEATURES "Persist\nParallelPersist\nMultiprocessing\n" \
+ "IrradianceCalc\nImmediateIrradiance\nDistanceLimiting\n" \
+ "ParticipatingMedia=Mist\n" \
+ "HessianAmbientCache\nAmbientAveraging\n" \
+ "AmbientValueSharing\nAdaptiveShadowTesting\n" \
+ "InputFormats=a,f,d\nOutputFormats=a,f,d,c\n" \
+ "Outputs=o,d,v,V,w,W,l,L,c,p,n,N,s,m,M,r,x,R,X,~\n"
+#else
+#define RTRACE_FEATURES "IrradianceCalc\nIrradianceCalc\nDistanceLimiting\n" \
+ "ParticipatingMedia=Mist\n" \
+ "HessianAmbientCache\nAmbientAveraging\n" \
+ "AmbientValueSharing\nAdaptiveShadowTesting\n" \
+ "InputFormats=a,f,d\nOutputFormats=a,f,d,c\n" \
+ "Outputs=o,d,v,V,w,W,l,L,c,p,n,N,s,m,M,r,x,R,X,~\n"
+#endif
+
+
int
-main(argc, argv)
-int argc;
-char *argv[];
+main(int argc, char *argv[])
{
#define check(ol,al) if (argv[i][ol] || \
badarg(argc-i-1,argv+i+1,al)) \
goto badopt
-#define bool(olen,var) switch (argv[i][olen]) { \
+#define check_bool(olen,var) switch (argv[i][olen]) { \
case '\0': var = !var; break; \
case 'y': case 'Y': case 't': case 'T': \
case '+': case '1': var = 1; break; \
case 'n': case 'N': case 'f': case 'F': \
case '-': case '0': var = 0; break; \
default: goto badopt; }
- int loadflags = ~IO_FILES;
+ extern char *octname;
int persist = 0;
- char **tralp;
- int duped1;
+ char *octnm = NULL;
+ char **tralp = NULL;
+ int duped1 = -1;
int rval;
int i;
/* global program name */
progname = argv[0] = fixargv0(argv[0]);
+ /* feature check only? */
+ strcat(RFeatureList, RTRACE_FEATURES);
+ if (argc > 1 && !strcmp(argv[1], "-features"))
+ return feature_status(argc-2, argv+2);
+ /* add trace notify function */
+ for (i = 0; addobjnotify[i] != NULL; i++)
+ ;
+ addobjnotify[i] = tranotify;
/* option city */
for (i = 1; i < argc; i++) {
/* expand arguments */
@@ -160,6 +144,12 @@ char *argv[];
continue;
}
switch (argv[i][1]) {
+ case 'n': /* number of cores */
+ check(2,"i");
+ nproc = atoi(argv[++i]);
+ if (nproc <= 0)
+ error(USER, "bad number of processes");
+ break;
case 'x': /* x resolution */
check(2,"i");
hresolu = atoi(argv[++i]);
@@ -168,9 +158,9 @@ char *argv[];
check(2,"i");
vresolu = atoi(argv[++i]);
break;
- case 'w': /* warnings */
+ case 'w': /* warnings & spectral */
rval = erract[WARNING].pf != NULL;
- bool(2,rval);
+ check_bool(2,rval);
if (rval) erract[WARNING].pf = wputs;
else erract[WARNING].pf = NULL;
break;
@@ -181,10 +171,10 @@ char *argv[];
case 'l': /* limit distance */
if (argv[i][2] != 'd')
goto badopt;
- bool(3,lim_dist);
+ check_bool(3,lim_dist);
break;
case 'I': /* immed. irradiance */
- bool(2,imm_irrad);
+ check_bool(2,imm_irrad);
break;
case 'f': /* format i/o */
switch (argv[i][2]) {
@@ -216,7 +206,7 @@ char *argv[];
break;
case 'h': /* header output */
rval = loadflags & IO_INFO;
- bool(2,rval);
+ check_bool(2,rval);
loadflags = rval ? loadflags | IO_INFO :
loadflags & ~IO_INFO;
break;
@@ -230,8 +220,8 @@ char *argv[];
tralp = tralist;
}
if (argv[i][2] == 'I') { /* file */
- rval = wordfile(tralp,
- getpath(argv[++i],getlibpath(),R_OK));
+ rval = wordfile(tralp, MAXMODLIST-(tralp-tralist),
+ getpath(argv[++i],getrlibpath(),R_OK));
if (rval < 0) {
sprintf(errmsg,
"cannot open trace include file \"%s\"",
@@ -252,8 +242,8 @@ char *argv[];
tralp = tralist;
}
if (argv[i][2] == 'E') { /* file */
- rval = wordfile(tralp,
- getpath(argv[++i],getlibpath(),R_OK));
+ rval = wordfile(tralp, MAXMODLIST-(tralp-tralist),
+ getpath(argv[++i],getrlibpath(),R_OK));
if (rval < 0) {
sprintf(errmsg,
"cannot open trace exclude file \"%s\"",
@@ -270,6 +260,74 @@ char *argv[];
goto badopt;
}
break;
+ case 'p': /* value output */
+ switch (argv[i][2]) {
+ case 'R': /* standard RGB output */
+ if (strcmp(argv[i]+2, "RGB"))
+ goto badopt;
+ out_prims = stdprims;
+ out_scalefactor = 1;
+ sens_curve = NULL;
+ break;
+ case 'X': /* XYZ output */
+ if (strcmp(argv[i]+2, "XYZ"))
+ goto badopt;
+ out_prims = xyzprims;
+ out_scalefactor = WHTEFFICACY;
+ sens_curve = NULL;
+ break;
+ case 'c': {
+ int j;
+ check(3,"ffffffff");
+ rval = 0;
+ for (j = 0; j < 8; j++) {
+ our_prims[0][j] = atof(argv[++i]);
+ rval |= fabs(our_prims[0][j]-stdprims[0][j]) > .001;
+ }
+ if (rval) {
+ if (!colorprimsOK(our_prims))
+ error(USER, "illegal primary chromaticities");
+ out_prims = our_prims;
+ } else
+ out_prims = stdprims;
+ out_scalefactor = 1;
+ sens_curve = NULL;
+ } break;
+ case 'Y': /* photopic response */
+ if (argv[i][3])
+ goto badopt;
+ sens_curve = scolor_photopic;
+ out_scalefactor = WHTEFFICACY;
+ break;
+ case 'S': /* scotopic response */
+ if (argv[i][3])
+ goto badopt;
+ sens_curve = scolor_scotopic;
+ out_scalefactor = WHTSCOTOPIC;
+ break;
+ case 'M': /* melanopic response */
+ if (argv[i][3])
+ goto badopt;
+ sens_curve = scolor_melanopic;
+ out_scalefactor = WHTMELANOPIC;
+ break;
+ default:
+ goto badopt;
+ }
+ break;
+#if MAXCSAMP>3
+ case 'c': /* output spectral results */
+ if (argv[i][2] != 'o')
+ goto badopt;
+ rval = (out_prims == NULL) & (sens_curve == NULL);
+ check_bool(3,rval);
+ if (rval) {
+ out_prims = NULL;
+ sens_curve = NULL;
+ } else if (out_prims == NULL)
+ out_prims = stdprims;
+ break;
+#endif
#ifdef PERSIST
case 'P': /* persist file */
if (argv[i][2] == 'P') {
@@ -286,16 +344,39 @@ char *argv[];
goto badopt;
}
}
+ /* set/check spectral sampling */
+ rval = setspectrsamp(CNDX, WLPART);
+ if (rval < 0)
+ error(USER, "unsupported spectral sampling");
+ if (out_prims != NULL) {
+ if (!rval)
+ error(WARNING, "spectral range incompatible with color output");
+ } else if (NCSAMP == 3)
+ out_prims = stdprims; /* 3 samples do not a spectrum make */
+ if (nproc > 1 && persist)
+ error(USER, "multiprocessing incompatible with persist file");
/* initialize object types */
initotypes();
/* initialize urand */
- initurand(2048);
+ if (rand_samp) {
+ srandom((long)time(0));
+ initurand(0);
+ } else {
+ srandom(0L);
+ initurand(2048);
+ }
/* set up signal handling */
sigdie(SIGINT, "Interrupt");
+#ifdef SIGHUP
sigdie(SIGHUP, "Hangup");
+#endif
sigdie(SIGTERM, "Terminate");
+#ifdef SIGPIPE
sigdie(SIGPIPE, "Broken pipe");
+#endif
+#ifdef SIGALRM
sigdie(SIGALRM, "Alarm clock");
+#endif
#ifdef SIGXCPU
sigdie(SIGXCPU, "CPU limit exceeded");
sigdie(SIGXFSZ, "File size exceeded");
@@ -315,12 +396,12 @@ char *argv[];
#endif
/* get octree */
if (i == argc)
- octname = NULL;
+ octnm = NULL;
else if (i == argc-1)
- octname = argv[i];
+ octnm = argv[i];
else
goto badopt;
- if (octname == NULL)
+ if (octnm == NULL)
error(USER, "missing octree argument");
/* set up output */
#ifdef PERSIST
@@ -329,12 +410,10 @@ char *argv[];
openheader();
}
#endif
-#ifdef MSDOS
if (outform != 'a')
- setmode(fileno(stdout), O_BINARY);
- if (octname == NULL)
- setmode(fileno(stdin), O_BINARY);
-#endif
+ SET_FILE_BINARY(stdout);
+ rval = setrtoutput();
+ octname = savqstr(octnm);
readoct(octname, loadflags, &thescene, NULL);
nsceneobjs = nobjects;
@@ -342,17 +421,32 @@ char *argv[];
printargs(i, argv, stdout);
printf("SOFTWARE= %s\n", VersionID);
fputnow(stdout);
+ if (rval > 0) /* saved from setrtoutput() call */
+ fputncomp(rval, stdout);
+ if (NCSAMP > 3)
+ fputwlsplit(WLPART, stdout);
+ if ((out_prims != stdprims) & (out_prims != NULL))
+ fputprims(out_prims, stdout);
+ if ((outform == 'f') | (outform == 'd'))
+ fputendian(stdout);
fputformat(formstr(outform), stdout);
- putchar('\n');
+ fputc('\n', stdout); /* end of header */
}
- marksources(); /* find and mark sources */
+ if (!castonly) { /* any actual ray traversal to do? */
- setambient(); /* initialize ambient calculation */
+ ray_init_pmap(); /* PMAP: set up & load photon maps */
+
+ marksources(); /* find and mark sources */
+ setambient(); /* initialize ambient calculation */
+ } else
+ distantsources(); /* else mark only distant sources */
+
+ fflush(stdout); /* in case we're duplicating header */
+
#ifdef PERSIST
if (persist) {
- fflush(stdout);
/* reconnect stdout */
dup2(duped1, fileno(stdout));
close(duped1);
@@ -363,10 +457,11 @@ char *argv[];
while ((rval=fork()) == 0) { /* keep on forkin' */
pflock(1);
pfhold();
+ ambsync(); /* load new values */
}
if (rval < 0)
error(SYSTEM, "cannot fork child for persist function");
- pfdetach(); /* parent exits */
+ pfdetach(); /* parent will run then exit */
}
}
runagain:
@@ -374,7 +469,7 @@ runagain:
dupheader(); /* send header to stdout */
#endif
/* trace rays */
- rtrace(NULL);
+ rtrace(NULL, nproc);
/* flush ambient file */
ambsync();
#ifdef PERSIST
@@ -389,26 +484,30 @@ runagain:
}
}
if (persist == PCHILD) { /* wait for a signal then go again */
- close(duped1); /* release output handle */
pfhold();
raynum = nrays = 0; /* reinitialize */
goto runagain;
}
#endif
+
+ ray_done_pmap(); /* PMAP: free photon maps */
+
quit(0);
badopt:
sprintf(errmsg, "command line error at '%s'", argv[i]);
error(USER, errmsg);
+ return 1; /* pro forma return */
#undef check
-#undef bool
+#undef check_bool
}
void
-wputs(s) /* warning output function */
-char *s;
+wputs( /* warning output function */
+ const char *s
+)
{
int lasterrno = errno;
eputs(s);
@@ -417,8 +516,9 @@ char *s;
void
-eputs(s) /* put string to stderr */
-register char *s;
+eputs( /* put string to stderr */
+ const char *s
+)
{
static int midline = 0;
@@ -436,17 +536,20 @@ register char *s;
}
-void
-onsig(signo) /* fatal signal */
-int signo;
+static void
+onsig( /* fatal signal */
+ int signo
+)
{
static int gotsig = 0;
if (gotsig++) /* two signals and we're gone! */
_exit(signo);
+#ifdef SIGALRM
alarm(15); /* allow 15 seconds to clean up */
signal(SIGALRM, SIG_DFL); /* make certain we do die */
+#endif
eputs("signal - ");
eputs(sigerr[signo]);
eputs("\n");
@@ -454,10 +557,11 @@ int signo;
}
-void
-sigdie(signo, msg) /* set fatal signal */
-int signo;
-char *msg;
+static void
+sigdie( /* set fatal signal */
+ int signo,
+ char *msg
+)
{
if (signal(signo, onsig) == SIG_IGN)
signal(signo, SIG_IGN);
@@ -465,26 +569,38 @@ char *msg;
}
-void
-printdefaults() /* print default values to stdout */
+static void
+printdefaults(void) /* print default values to stdout */
{
- register char *cp;
+ char *cp;
+ printf(erract[WARNING].pf != NULL ?
+ "-w+\t\t\t\t# warning messages on\n" :
+ "-w-\t\t\t\t# warning messages off\n");
if (imm_irrad)
printf("-I+\t\t\t\t# immediate irradiance on\n");
- printf("-x %-9d\t\t\t# x resolution\n", hresolu);
- printf("-y %-9d\t\t\t# y resolution\n", vresolu);
+ printf("-n %-2d\t\t\t\t# number of rendering processes\n", nproc);
+ printf("-x %-9d\t\t\t# %s\n", hresolu,
+ vresolu && hresolu ? "x resolution" : "flush interval");
+ printf("-y %-9d\t\t\t# y resolution\n", vresolu);
printf(lim_dist ? "-ld+\t\t\t\t# limit distance on\n" :
"-ld-\t\t\t\t# limit distance off\n");
+ printf(loadflags & IO_INFO ? "-h+\t\t\t\t# output header\n" :
+ "-h-\t\t\t\t# no header\n");
printf("-f%c%c\t\t\t\t# format input/output = %s/%s\n",
inform, outform, formstr(inform), formstr(outform));
- printf("-o%s\t\t\t\t# output", outvals);
+ printf("-o%-9s\t\t\t# output", outvals);
for (cp = outvals; *cp; cp++)
switch (*cp) {
- case 't': printf(" trace"); break;
+ case 't': case 'T': printf(" trace"); break;
case 'o': printf(" origin"); break;
case 'd': printf(" direction"); break;
+ case 'r': printf(" reflect_contrib"); break;
+ case 'R': printf(" reflect_length"); break;
+ case 'x': printf(" unreflect_contrib"); break;
+ case 'X': printf(" unreflect_length"); break;
case 'v': printf(" value"); break;
+ case 'V': printf(" contribution"); break;
case 'l': printf(" length"); break;
case 'L': printf(" first_length"); break;
case 'p': printf(" point"); break;
@@ -492,11 +608,30 @@ printdefaults() /* print default values to stdout */
case 'N': printf(" unperturbed_normal"); break;
case 's': printf(" surface"); break;
case 'w': printf(" weight"); break;
+ case 'W': printf(" coefficient"); break;
case 'm': printf(" modifier"); break;
+ case 'M': printf(" material"); break;
+ case '~': printf(" tilde"); break;
}
- putchar('\n');
- printf(erract[WARNING].pf != NULL ?
- "-w+\t\t\t\t# warning messages on\n" :
- "-w-\t\t\t\t# warning messages off\n");
+ fputc('\n', stdout);
+ if (sens_curve == scolor_photopic)
+ printf("-pY\t\t\t\t# photopic output\n");
+ else if (sens_curve == scolor_scotopic)
+ printf("-pS\t\t\t\t# scotopic output\n");
+ else if (sens_curve == scolor_melanopic)
+ printf("-pM\t\t\t\t# melanopic output\n");
+ else if (out_prims == stdprims)
+ printf("-pRGB\t\t\t\t# standard RGB color output\n");
+ else if (out_prims == xyzprims)
+ printf("-pXYZ\t\t\t\t# CIE XYZ color output\n");
+ else if (out_prims != NULL)
+ printf("-pc %.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f\t# output color primaries and white point\n",
+ out_prims[RED][0], out_prims[RED][1],
+ out_prims[GRN][0], out_prims[GRN][1],
+ out_prims[BLU][0], out_prims[BLU][1],
+ out_prims[WHT][0], out_prims[WHT][1]);
+ if (NCSAMP > 3)
+ printf(out_prims != NULL ? "-co-\t\t\t\t# output tristimulus colors\n" :
+ "-co+\t\t\t\t# output spectral values\n");
print_rdefaults();
}