ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rtmain.c
Revision: 2.31
Committed: Sun Apr 5 15:47:02 2020 UTC (4 years, 1 month ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.30: +6 -2 lines
Log Message:
Added NCOMP= to header and better checking of -o* options

File Contents

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