ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rtmain.c
Revision: 2.37
Committed: Tue Oct 18 18:07:50 2022 UTC (18 months, 3 weeks ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.36: +2 -2 lines
Log Message:
fix(rtrace): corrected defaults reporting with tilde

File Contents

# User Rev Content
1 greg 2.1 #ifndef lint
2 greg 2.37 static const char RCSid[] = "$Id: rtmain.c,v 2.36 2022/04/06 15:28:18 greg Exp $";
3 greg 2.1 #endif
4     /*
5     * rtmain.c - main for rtrace per-ray calculation program
6     */
7    
8 greg 2.2 #include "copyright.h"
9 greg 2.1
10 schorsch 2.4 #include <signal.h>
11    
12     #include "platform.h"
13 schorsch 2.7 #include "rtprocess.h" /* getpid() */
14 schorsch 2.9 #include "resolu.h"
15 greg 2.1 #include "ray.h"
16     #include "source.h"
17     #include "ambient.h"
18     #include "random.h"
19     #include "paths.h"
20 greg 2.23 #include "pmapray.h"
21 greg 2.1
22 greg 2.18 extern char *progname; /* global argv[0] */
23    
24     extern char *shm_boundary; /* boundary of shared memory */
25    
26 greg 2.1 /* persistent processes define */
27     #ifdef F_SETLKW
28     #define PERSIST 1 /* normal persist */
29     #define PARALLEL 2 /* parallel persist */
30     #define PCHILD 3 /* child of normal persist */
31     #endif
32    
33     char *sigerr[NSIG]; /* signal error messages */
34     char *errfile = NULL; /* error output file */
35    
36 greg 2.18 int nproc = 1; /* number of processes */
37    
38 greg 2.31 extern char *formstr(int f); /* string from format */
39     extern int setrtoutput(void); /* set output values */
40 greg 2.32
41 greg 2.18 int inform = 'a'; /* input format */
42     int outform = 'a'; /* output format */
43     char *outvals = "v"; /* output specification */
44 greg 2.1
45 greg 2.18 int hresolu = 0; /* horizontal (scan) size */
46     int vresolu = 0; /* vertical resolution */
47 greg 2.1
48 greg 2.33 extern int castonly; /* only doing ray-casting? */
49 greg 2.32
50 greg 2.18 int imm_irrad = 0; /* compute immediate irradiance? */
51     int lim_dist = 0; /* limit distance? */
52    
53     #ifndef MAXMODLIST
54     #define MAXMODLIST 1024 /* maximum modifiers we'll track */
55     #endif
56 greg 2.1
57 greg 2.18 extern void (*addobjnotify[])(); /* object notification calls */
58     extern void tranotify(OBJECT obj);
59    
60     char *tralist[MAXMODLIST]; /* list of modifers to trace (or no) */
61     int traincl = -1; /* include == 1, exclude == 0 */
62 greg 2.1
63 greg 2.10 static int loadflags = ~IO_FILES; /* what to load from octree */
64    
65 schorsch 2.9 static void onsig(int signo);
66     static void sigdie(int signo, char *msg);
67     static void printdefaults(void);
68 greg 2.1
69    
70     int
71 schorsch 2.9 main(int argc, char *argv[])
72 greg 2.1 {
73     #define check(ol,al) if (argv[i][ol] || \
74     badarg(argc-i-1,argv+i+1,al)) \
75     goto badopt
76 schorsch 2.26 #define check_bool(olen,var) switch (argv[i][olen]) { \
77 greg 2.1 case '\0': var = !var; break; \
78     case 'y': case 'Y': case 't': case 'T': \
79     case '+': case '1': var = 1; break; \
80     case 'n': case 'N': case 'f': case 'F': \
81     case '-': case '0': var = 0; break; \
82     default: goto badopt; }
83 greg 2.25 extern char *octname;
84 greg 2.1 int persist = 0;
85 greg 2.18 char *octnm = NULL;
86 greg 2.22 char **tralp = NULL;
87     int duped1 = -1;
88 greg 2.1 int rval;
89     int i;
90     /* global program name */
91     progname = argv[0] = fixargv0(argv[0]);
92 greg 2.18 /* add trace notify function */
93     for (i = 0; addobjnotify[i] != NULL; i++)
94     ;
95     addobjnotify[i] = tranotify;
96 greg 2.1 /* option city */
97     for (i = 1; i < argc; i++) {
98     /* expand arguments */
99     while ((rval = expandarg(&argc, &argv, i)) > 0)
100     ;
101     if (rval < 0) {
102     sprintf(errmsg, "cannot expand '%s'", argv[i]);
103     error(SYSTEM, errmsg);
104     }
105     if (argv[i] == NULL || argv[i][0] != '-')
106     break; /* break from options */
107     if (!strcmp(argv[i], "-version")) {
108     puts(VersionID);
109     quit(0);
110     }
111     if (!strcmp(argv[i], "-defaults") ||
112     !strcmp(argv[i], "-help")) {
113     printdefaults();
114     quit(0);
115     }
116     rval = getrenderopt(argc-i, argv+i);
117     if (rval >= 0) {
118     i += rval;
119     continue;
120     }
121     switch (argv[i][1]) {
122 greg 2.18 case 'n': /* number of cores */
123     check(2,"i");
124     nproc = atoi(argv[++i]);
125     if (nproc <= 0)
126     error(USER, "bad number of processes");
127     break;
128 greg 2.1 case 'x': /* x resolution */
129     check(2,"i");
130     hresolu = atoi(argv[++i]);
131     break;
132     case 'y': /* y resolution */
133     check(2,"i");
134     vresolu = atoi(argv[++i]);
135     break;
136     case 'w': /* warnings */
137     rval = erract[WARNING].pf != NULL;
138 schorsch 2.26 check_bool(2,rval);
139 greg 2.1 if (rval) erract[WARNING].pf = wputs;
140     else erract[WARNING].pf = NULL;
141     break;
142     case 'e': /* error file */
143     check(2,"s");
144     errfile = argv[++i];
145     break;
146     case 'l': /* limit distance */
147     if (argv[i][2] != 'd')
148     goto badopt;
149 schorsch 2.26 check_bool(3,lim_dist);
150 greg 2.1 break;
151     case 'I': /* immed. irradiance */
152 schorsch 2.26 check_bool(2,imm_irrad);
153 greg 2.1 break;
154     case 'f': /* format i/o */
155     switch (argv[i][2]) {
156     case 'a': /* ascii */
157     case 'f': /* float */
158     case 'd': /* double */
159     inform = argv[i][2];
160     break;
161     default:
162     goto badopt;
163     }
164     switch (argv[i][3]) {
165     case '\0':
166     outform = inform;
167     break;
168     case 'a': /* ascii */
169     case 'f': /* float */
170     case 'd': /* double */
171     case 'c': /* color */
172     check(4,"");
173     outform = argv[i][3];
174     break;
175     default:
176     goto badopt;
177     }
178     break;
179     case 'o': /* output */
180     outvals = argv[i]+2;
181     break;
182     case 'h': /* header output */
183     rval = loadflags & IO_INFO;
184 schorsch 2.26 check_bool(2,rval);
185 greg 2.1 loadflags = rval ? loadflags | IO_INFO :
186     loadflags & ~IO_INFO;
187     break;
188     case 't': /* trace */
189     switch (argv[i][2]) {
190     case 'i': /* include */
191     case 'I':
192     check(3,"s");
193     if (traincl != 1) {
194     traincl = 1;
195     tralp = tralist;
196     }
197     if (argv[i][2] == 'I') { /* file */
198 greg 2.27 rval = wordfile(tralp, MAXMODLIST-(tralp-tralist),
199 greg 2.3 getpath(argv[++i],getrlibpath(),R_OK));
200 greg 2.1 if (rval < 0) {
201     sprintf(errmsg,
202     "cannot open trace include file \"%s\"",
203     argv[i]);
204     error(SYSTEM, errmsg);
205     }
206     tralp += rval;
207     } else {
208     *tralp++ = argv[++i];
209     *tralp = NULL;
210     }
211     break;
212     case 'e': /* exclude */
213     case 'E':
214     check(3,"s");
215     if (traincl != 0) {
216     traincl = 0;
217     tralp = tralist;
218     }
219     if (argv[i][2] == 'E') { /* file */
220 greg 2.27 rval = wordfile(tralp, MAXMODLIST-(tralp-tralist),
221 greg 2.3 getpath(argv[++i],getrlibpath(),R_OK));
222 greg 2.1 if (rval < 0) {
223     sprintf(errmsg,
224     "cannot open trace exclude file \"%s\"",
225     argv[i]);
226     error(SYSTEM, errmsg);
227     }
228     tralp += rval;
229     } else {
230     *tralp++ = argv[++i];
231     *tralp = NULL;
232     }
233     break;
234     default:
235     goto badopt;
236     }
237     break;
238     #ifdef PERSIST
239     case 'P': /* persist file */
240     if (argv[i][2] == 'P') {
241     check(3,"s");
242     persist = PARALLEL;
243     } else {
244     check(2,"s");
245     persist = PERSIST;
246     }
247     persistfile(argv[++i]);
248     break;
249     #endif
250     default:
251     goto badopt;
252     }
253     }
254 greg 2.34 if (nproc > 1 && persist)
255     error(USER, "multiprocessing incompatible with persist file");
256 greg 2.1 /* initialize object types */
257     initotypes();
258     /* initialize urand */
259 greg 2.14 if (rand_samp) {
260     srandom((long)time(0));
261     initurand(0);
262     } else {
263     srandom(0L);
264     initurand(2048);
265     }
266 greg 2.1 /* set up signal handling */
267     sigdie(SIGINT, "Interrupt");
268 schorsch 2.5 #ifdef SIGHUP
269 greg 2.1 sigdie(SIGHUP, "Hangup");
270 schorsch 2.5 #endif
271 greg 2.1 sigdie(SIGTERM, "Terminate");
272 schorsch 2.5 #ifdef SIGPIPE
273 greg 2.1 sigdie(SIGPIPE, "Broken pipe");
274 schorsch 2.5 #endif
275     #ifdef SIGALRM
276 greg 2.1 sigdie(SIGALRM, "Alarm clock");
277 schorsch 2.5 #endif
278 greg 2.1 #ifdef SIGXCPU
279     sigdie(SIGXCPU, "CPU limit exceeded");
280     sigdie(SIGXFSZ, "File size exceeded");
281     #endif
282     /* open error file */
283     if (errfile != NULL) {
284     if (freopen(errfile, "a", stderr) == NULL)
285     quit(2);
286     fprintf(stderr, "**************\n*** PID %5d: ",
287     getpid());
288     printargs(argc, argv, stderr);
289     putc('\n', stderr);
290     fflush(stderr);
291     }
292     #ifdef NICE
293     nice(NICE); /* lower priority */
294     #endif
295     /* get octree */
296     if (i == argc)
297 greg 2.18 octnm = NULL;
298 greg 2.1 else if (i == argc-1)
299 greg 2.18 octnm = argv[i];
300 greg 2.1 else
301     goto badopt;
302 greg 2.18 if (octnm == NULL)
303 greg 2.1 error(USER, "missing octree argument");
304     /* set up output */
305     #ifdef PERSIST
306     if (persist) {
307     duped1 = dup(fileno(stdout)); /* don't lose our output */
308     openheader();
309     }
310     #endif
311     if (outform != 'a')
312 schorsch 2.4 SET_FILE_BINARY(stdout);
313 greg 2.31 rval = setrtoutput();
314 greg 2.24 readoct(octname = octnm, loadflags, &thescene, NULL);
315 greg 2.1 nsceneobjs = nobjects;
316    
317     if (loadflags & IO_INFO) { /* print header */
318     printargs(i, argv, stdout);
319     printf("SOFTWARE= %s\n", VersionID);
320     fputnow(stdout);
321 greg 2.31 if (rval > 0) /* saved from setrtoutput() call */
322     printf("NCOMP=%d\n", rval);
323 greg 2.30 if ((outform == 'f') | (outform == 'd'))
324     fputendian(stdout);
325 greg 2.1 fputformat(formstr(outform), stdout);
326     putchar('\n');
327     }
328    
329 greg 2.32 if (!castonly) { /* any actual ray traversal to do? */
330    
331     ray_init_pmap(); /* PMAP: set up & load photon maps */
332    
333     marksources(); /* find and mark sources */
334    
335     setambient(); /* initialize ambient calculation */
336 greg 2.33 } else
337     distantsources(); /* else mark only distant sources */
338 greg 2.35
339 greg 2.36 fflush(stdout); /* in case we're duplicating header */
340 greg 2.35
341 greg 2.1 #ifdef PERSIST
342     if (persist) {
343     /* reconnect stdout */
344     dup2(duped1, fileno(stdout));
345     close(duped1);
346     if (persist == PARALLEL) { /* multiprocessing */
347     preload_objs(); /* preload scene */
348     shm_boundary = (char *)malloc(16);
349     strcpy(shm_boundary, "SHM_BOUNDARY");
350     while ((rval=fork()) == 0) { /* keep on forkin' */
351     pflock(1);
352     pfhold();
353 greg 2.15 ambsync(); /* load new values */
354 greg 2.1 }
355     if (rval < 0)
356     error(SYSTEM, "cannot fork child for persist function");
357 greg 2.15 pfdetach(); /* parent will run then exit */
358 greg 2.1 }
359     }
360     runagain:
361     if (persist)
362     dupheader(); /* send header to stdout */
363     #endif
364     /* trace rays */
365 greg 2.19 rtrace(NULL, nproc);
366 greg 2.1 /* flush ambient file */
367     ambsync();
368     #ifdef PERSIST
369     if (persist == PERSIST) { /* first run-through */
370     if ((rval=fork()) == 0) { /* child loops until killed */
371     pflock(1);
372     persist = PCHILD;
373     } else { /* original process exits */
374     if (rval < 0)
375     error(SYSTEM, "cannot fork child for persist function");
376     pfdetach(); /* parent exits */
377     }
378     }
379     if (persist == PCHILD) { /* wait for a signal then go again */
380     pfhold();
381     raynum = nrays = 0; /* reinitialize */
382     goto runagain;
383     }
384     #endif
385 greg 2.23
386     ray_done_pmap(); /* PMAP: free photon maps */
387    
388 greg 2.1 quit(0);
389    
390     badopt:
391     sprintf(errmsg, "command line error at '%s'", argv[i]);
392     error(USER, errmsg);
393 schorsch 2.9 return 1; /* pro forma return */
394 greg 2.1
395     #undef check
396 schorsch 2.26 #undef check_bool
397 greg 2.1 }
398    
399    
400     void
401 schorsch 2.9 wputs( /* warning output function */
402     char *s
403     )
404 greg 2.1 {
405     int lasterrno = errno;
406     eputs(s);
407     errno = lasterrno;
408     }
409    
410    
411     void
412 schorsch 2.9 eputs( /* put string to stderr */
413     register char *s
414     )
415 greg 2.1 {
416     static int midline = 0;
417    
418     if (!*s)
419     return;
420     if (!midline++) {
421     fputs(progname, stderr);
422     fputs(": ", stderr);
423     }
424     fputs(s, stderr);
425     if (s[strlen(s)-1] == '\n') {
426     fflush(stderr);
427     midline = 0;
428     }
429     }
430    
431    
432 schorsch 2.9 static void
433     onsig( /* fatal signal */
434     int signo
435     )
436 greg 2.1 {
437     static int gotsig = 0;
438    
439     if (gotsig++) /* two signals and we're gone! */
440     _exit(signo);
441    
442 schorsch 2.5 #ifdef SIGALRM
443 greg 2.1 alarm(15); /* allow 15 seconds to clean up */
444     signal(SIGALRM, SIG_DFL); /* make certain we do die */
445 schorsch 2.5 #endif
446 greg 2.1 eputs("signal - ");
447     eputs(sigerr[signo]);
448     eputs("\n");
449     quit(3);
450     }
451    
452    
453 schorsch 2.9 static void
454     sigdie( /* set fatal signal */
455     int signo,
456     char *msg
457     )
458 greg 2.1 {
459     if (signal(signo, onsig) == SIG_IGN)
460     signal(signo, SIG_IGN);
461     sigerr[signo] = msg;
462     }
463    
464    
465 schorsch 2.9 static void
466 schorsch 2.8 printdefaults(void) /* print default values to stdout */
467 greg 2.1 {
468     register char *cp;
469    
470     if (imm_irrad)
471     printf("-I+\t\t\t\t# immediate irradiance on\n");
472 greg 2.18 printf("-n %-2d\t\t\t\t# number of rendering processes\n", nproc);
473 greg 2.21 printf("-x %-9d\t\t\t# %s\n", hresolu,
474     vresolu && hresolu ? "x resolution" : "flush interval");
475 greg 2.16 printf("-y %-9d\t\t\t# y resolution\n", vresolu);
476 greg 2.1 printf(lim_dist ? "-ld+\t\t\t\t# limit distance on\n" :
477     "-ld-\t\t\t\t# limit distance off\n");
478 greg 2.10 printf("-h%c\t\t\t\t# %s header\n", loadflags & IO_INFO ? '+' : '-',
479     loadflags & IO_INFO ? "output" : "no");
480 greg 2.1 printf("-f%c%c\t\t\t\t# format input/output = %s/%s\n",
481     inform, outform, formstr(inform), formstr(outform));
482 greg 2.10 printf("-o%-9s\t\t\t# output", outvals);
483 greg 2.1 for (cp = outvals; *cp; cp++)
484     switch (*cp) {
485 greg 2.10 case 't': case 'T': printf(" trace"); break;
486 greg 2.1 case 'o': printf(" origin"); break;
487     case 'd': printf(" direction"); break;
488 greg 2.28 case 'r': printf(" reflect_contrib"); break;
489     case 'R': printf(" reflect_length"); break;
490     case 'x': printf(" unreflect_contrib"); break;
491     case 'X': printf(" unreflect_length"); break;
492 greg 2.1 case 'v': printf(" value"); break;
493 greg 2.13 case 'V': printf(" contribution"); break;
494 greg 2.1 case 'l': printf(" length"); break;
495     case 'L': printf(" first_length"); break;
496     case 'p': printf(" point"); break;
497     case 'n': printf(" normal"); break;
498     case 'N': printf(" unperturbed_normal"); break;
499     case 's': printf(" surface"); break;
500     case 'w': printf(" weight"); break;
501 greg 2.13 case 'W': printf(" coefficient"); break;
502 greg 2.1 case 'm': printf(" modifier"); break;
503 greg 2.10 case 'M': printf(" material"); break;
504 greg 2.37 case '~': printf(" tilde"); break;
505 greg 2.1 }
506     putchar('\n');
507     printf(erract[WARNING].pf != NULL ?
508     "-w+\t\t\t\t# warning messages on\n" :
509     "-w-\t\t\t\t# warning messages off\n");
510     print_rdefaults();
511     }