ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/mkillum.c
(Generate patch)

Comparing ray/src/gen/mkillum.c (file contents):
Revision 2.27 by schorsch, Wed Jun 7 17:52:04 2006 UTC vs.
Revision 2.31 by greg, Fri Sep 21 05:53:21 2007 UTC

# Line 23 | Line 23 | static const char RCSid[] = "$Id$";
23   #define  S_COMPL        2               /* select all but element */
24   #define  S_ALL          3               /* select all */
25  
26 –                                /* rtrace command and defaults */
27 – char  *rtargv[64] = { "rtrace", "-dj", ".25", "-dr", "3", "-dv-",
28 –                "-ab", "2", "-ad", "1024", "-as", "512", "-aa", ".1", };
29 – int  rtargc = 14;
30 –                                /* overriding rtrace options */
31 – char  *myrtopts[] = { "-I-", "-i-", "-ld-", "-ov", "-h-",
32 –                        "-fff", "-y", "0", NULL };
33 –
34 – struct rtproc   rt0;            /* head of rtrace process list */
35 –
26   struct illum_args  thisillum = {        /* our illum and default values */
27                  0,
28 +                UDzpos,
29 +                0.,
30                  DFLMAT,
31                  DFLDAT,
32                  0,
33                  VOIDID,
34                  SAMPDENS,
35                  NSAMPS,
36 +                NULL,
37                  0.,
38          };
39  
40   char    matcheck[MAXSTR];       /* current material to include or exclude */
41   int     matselect = S_ALL;      /* selection criterion */
42  
50 – FUN     ofun[NUMOTYPE] = INIT_OTYPE;    /* object types */
51 –
52 – char    persistfn[] = "pfXXXXXX";       /* persist file name */
53 –
43   int     gargc;                  /* global argc */
44   char    **gargv;                /* global argv */
56 – #define  progname       gargv[0]
45  
46   int     doneheader = 0;         /* printed header yet? */
47   #define  checkhead()    if (!doneheader++) printhead(gargc,gargv)
48  
49   int     warnings = 1;           /* print warnings? */
50  
51 < int done_rprocs(struct rtproc *rtp);
64 < void init(int np);
51 > void init(char *octnm, int np);
52   void filter(register FILE       *infp, char     *name);
53   void xoptions(char      *s, char        *nm);
54   void printopts(void);
# Line 76 | Line 63 | main(          /* compute illum distributions using rtrace */
63   )
64   {
65          int     nprocs = 1;
79 –        char    *rtpath;
66          FILE    *fp;
67 +        int     rval;
68          register int    i;
69                                  /* set global arguments */
70          gargv = argv;
71 <                                /* check for -n option */
72 <        if (!strcmp(argv[1], "-n")) {
73 <                nprocs = atoi(argv[2]);
74 <                if (nprocs <= 0)
75 <                        error(USER, "illegal number of processes");
76 <                i = 3;
77 <        } else
78 <                i = 1;
79 <                                /* set up rtrace command */
80 <        for ( ; i < argc; i++) {
81 <                if (argv[i][0] == '<' && argv[i][1] == '\0')
71 >        progname = gargv[0];
72 >                                /* set up rendering defaults */
73 >        dstrsrc = 0.25;
74 >        directrelay = 3;
75 >        directvis = 0;
76 >        ambounce = 2;
77 >                                /* get options from command line */
78 >        for (i = 1; i < argc; i++) {
79 >                while ((rval = expandarg(&argc, &argv, i)) > 0)
80 >                        ;
81 >                if (rval < 0) {
82 >                        sprintf(errmsg, "cannot expand '%s'", argv[i]);
83 >                        error(SYSTEM, errmsg);
84 >                }
85 >                if (argv[i][0] != '-')
86                          break;
87 <                rtargv[rtargc++] = argv[i];
88 <                if (argv[i][0] == '-' && argv[i][1] == 'w')
89 <                        switch (argv[i][2]) {
99 <                        case '\0':
100 <                                warnings = !warnings;
101 <                                break;
102 <                        case '+':
103 <                        case 'T': case 't':
104 <                        case 'Y': case 'y':
105 <                        case '1':
106 <                                warnings = 1;
107 <                                break;
108 <                        case '-':
109 <                        case 'F': case 'f':
110 <                        case 'N': case 'n':
111 <                        case '0':
112 <                                warnings = 0;
113 <                                break;
114 <                        }
115 <        }
116 <        gargc = i;
117 <        if (!strcmp(rtargv[--rtargc], "-defaults"))
118 <                nprocs = 0;
119 <        if (nprocs > 1) {       /* add persist file if parallel invocation */
120 <                rtargv[rtargc++] = "-PP";
121 <                rtargv[rtargc++] = mktemp(persistfn);
122 <        }
123 <                                /* add "mandatory" rtrace options */
124 <        for (i = 0; myrtopts[i] != NULL; i++)
125 <                rtargv[rtargc++] = myrtopts[i];
126 <                                /* finally, put back final argument */
127 <        rtargv[rtargc++] = argv[gargc-1];
128 <        rtargv[rtargc] = NULL;
129 <        if (!nprocs) {          /* just asking for defaults? */
130 <                printopts(); fflush(stdout);
131 <                rtpath = getpath(rtargv[0], getenv("PATH"), X_OK);
132 <                if (rtpath == NULL) {
133 <                        eputs(rtargv[0]);
134 <                        eputs(": command not found\n");
135 <                        exit(1);
87 >                if (!strcmp(argv[i], "-w")) {
88 >                        warnings = 0;
89 >                        continue;
90                  }
91 <                execv(rtpath, rtargv);
92 <                perror(rtpath); /* execv() should not return */
93 <                exit(1);
91 >                if (!strcmp(argv[i], "-n")) {
92 >                        nprocs = atoi(argv[++i]);
93 >                        if (nprocs <= 0)
94 >                                error(USER, "illegal number of processes");
95 >                        continue;
96 >                }
97 >                if (!strcmp(argv[i], "-defaults")) {
98 >                        printopts();
99 >                        print_rdefaults();
100 >                        quit(0);
101 >                }
102 >                rval = getrenderopt(argc-i, argv+i);
103 >                if (rval < 0) {
104 >                        sprintf(errmsg, "bad render option at '%s'", argv[i]);
105 >                        error(USER, errmsg);
106 >                }
107 >                i += rval;
108          }
109 <        if (gargc < 2 || argv[gargc-1][0] == '-')
109 >        gargc = ++i;
110 >                                /* add "mandatory" render options */
111 >        do_irrad = 0;
112 >        if (gargc > argc || argv[gargc-1][0] == '-')
113                  error(USER, "missing octree argument");
114                                  /* else initialize and run our calculation */
115 <        init(nprocs);
116 <        if (gargc+1 < argc)
115 >        init(argv[gargc-1], nprocs);
116 >        if (gargc < argc) {
117 >                if (gargc == argc-1 || argv[gargc][0] != '<' || argv[gargc][1])
118 >                        error(USER, "use '< file1 file2 ..' for multiple inputs");
119                  for (i = gargc+1; i < argc; i++) {
120                          if ((fp = fopen(argv[i], "r")) == NULL) {
121                                  sprintf(errmsg,
# Line 152 | Line 125 | main(          /* compute illum distributions using rtrace */
125                          filter(fp, argv[i]);
126                          fclose(fp);
127                  }
128 <        else
128 >        } else
129                  filter(stdin, "standard input");
130          quit(0);
131          return 0; /* pro forma return */
132   }
133  
134  
162 – #ifndef SIGALRM
163 – #define SIGALRM SIGTERM
164 – #endif
165 – static void
166 – killpersist(void)                       /* kill persistent rtrace process */
167 – {
168 –        FILE    *fp = fopen(persistfn, "r");
169 –        int     pid;
170 –
171 –        if (fp == NULL)
172 –                return;
173 –        if (fscanf(fp, "%*s %d", &pid) != 1 || kill(pid, SIGALRM) < 0)
174 –                unlink(persistfn);
175 –        fclose(fp);
176 – }
177 –
178 –
179 – int
180 – done_rprocs(struct rtproc *rtp)
181 – {
182 –        int     st0, st1 = 0;
183 –
184 –        if (rtp->next != NULL) {        /* close last opened first! */
185 –                st1 = done_rprocs(rtp->next);
186 –                free((void *)rtp->next);
187 –                rtp->next = NULL;
188 –        }
189 –        st0 = close_process(&rtp->pd);
190 –        if (st0 < 0)
191 –                error(WARNING, "unknown return status from rtrace process");
192 –        else if (st0 > 0)
193 –                return(st0);
194 –        return(st1);
195 – }
196 –
135   void
136 < quit(int status)                        /* exit with status */
136 > init(char *octnm, int np)               /* start rendering process(es) */
137   {
200 –        int     rtstat;
201 –
202 –        if (rt0.next != NULL)           /* terminate persistent rtrace */
203 –                killpersist();
204 –                                        /* clean up rtrace process(es) */
205 –        rtstat = done_rprocs(&rt0);
206 –        if (status == 0)
207 –                status = rtstat;
208 –        exit(status);
209 – }
210 –
211 – void
212 – init(int np)                            /* start rtrace and set up buffers */
213 – {
214 –        struct rtproc   *rtp;
215 –        int     i;
216 –        int     maxbytes;
217 –                                        /* set up object functions */
218 –        ofun[OBJ_FACE].funp = o_face;
219 –        ofun[OBJ_SPHERE].funp = o_sphere;
220 –        ofun[OBJ_RING].funp = o_ring;
138                                          /* set up signal handling */
139          signal(SIGINT, quit);
140   #ifdef SIGHUP
# Line 229 | Line 146 | init(int np)                           /* start rtrace and set up buffers */
146   #ifdef SIGPIPE
147          signal(SIGPIPE, quit);
148   #endif
149 <        rtp = &rt0;                     /* start rtrace process(es) */
150 <        for (i = 0; i++ < np; ) {
234 <                errno = 0;
235 <                maxbytes = open_process(&rtp->pd, rtargv);
236 <                if (maxbytes == 0) {
237 <                        eputs(rtargv[0]);
238 <                        eputs(": command not found\n");
239 <                        exit(1);
240 <                }
241 <                if (maxbytes < 0)
242 <                        error(SYSTEM, "cannot start rtrace process");
243 <                rtp->bsiz = maxbytes/(6*sizeof(float));
244 <                rtp->buf = (float *)malloc(6*sizeof(float)*rtp->bsiz--);
245 <                rtp->dest = (float **)calloc(rtp->bsiz, sizeof(float *));
246 <                if (rtp->buf == NULL || rtp->dest == NULL)
247 <                        error(SYSTEM, "out of memory in init");
248 <                rtp->nrays = 0;
249 <                if (i == np)            /* last process? */
250 <                        break;
251 <                if (np > 1)
252 <                        sleep(2);       /* wait for persist file */
253 <                rtp->next = (struct rtproc *)malloc(sizeof(struct rtproc));
254 <                if (rtp->next == NULL)
255 <                        error(SYSTEM, "out of memory in init");
256 <                rtp = rtp->next;
257 <        }
258 <        rtp->next = NULL;
259 <                                        /* set up urand */
260 <        initurand(16384);
149 >                                        /* start rendering process(es) */
150 >        ray_pinit(octnm, np);
151   }
152  
153  
# Line 329 | Line 219 | xoptions(                      /* process options in string s */
219   {
220          extern FILE     *freopen();
221          char    buf[64];
222 +        int     negax;
223          int     nerrs = 0;
224          register char   *cp;
225  
# Line 408 | Line 299 | xoptions(                      /* process options in string s */
299                          }
300                          cp = sskip(cp);
301                          continue;
302 <                case 'd':                       /* point sample density */
302 >                case 'd':                       /* sample density / BSDF data */
303                          if (*++cp != '=')
304                                  break;
305 <                        if (!isintd(++cp, " \t\n\r"))
306 <                                break;
307 <                        thisillum.sampdens = atoi(cp);
305 >                        if (thisillum.sd != NULL) {
306 >                                free_BSDF(thisillum.sd);
307 >                                thisillum.sd = NULL;
308 >                        }
309 >                        if (!*++cp || isspace(*cp))
310 >                                continue;
311 >                        if (isintd(++cp, " \t\n\r")) {
312 >                                thisillum.sampdens = atoi(cp);
313 >                        } else {
314 >                                atos(buf, sizeof(buf), cp);
315 >                                thisillum.sd = load_BSDF(buf);
316 >                                if (thisillum.sd == NULL)
317 >                                        break;
318 >                        }
319                          cp = sskip(cp);
320                          continue;
321 <                case 's':                       /* point super-samples */
321 >                case 's':                       /* surface super-samples */
322                          if (*++cp != '=')
323                                  break;
324                          if (!isintd(++cp, " \t\n\r"))
# Line 458 | Line 360 | xoptions(                      /* process options in string s */
360                          }
361                          doneheader = 0;
362                          continue;
363 +                case 'u':                       /* up direction */
364 +                        if (*++cp != '=')
365 +                                break;
366 +                        if (!*++cp || isspace(*cp)) {
367 +                                thisillum.udir = UDunknown;
368 +                                continue;
369 +                        }
370 +                        negax = 0;
371 +                        if (*cp == '+')
372 +                                cp++;
373 +                        else if (*cp == '-') {
374 +                                negax++;
375 +                                cp++;
376 +                        }
377 +                        switch (*cp++) {
378 +                        case 'x':
379 +                        case 'X':
380 +                                thisillum.udir = negax ? UDxneg : UDxpos;
381 +                                break;
382 +                        case 'y':
383 +                        case 'Y':
384 +                                thisillum.udir = negax ? UDyneg : UDypos;
385 +                                break;
386 +                        case 'z':
387 +                        case 'Z':
388 +                                thisillum.udir = negax ? UDxneg : UDxpos;
389 +                                break;
390 +                        default:
391 +                                thisillum.udir = UDunknown;
392 +                                break;
393 +                        }
394 +                        if (thisillum.udir == UDunknown || !isspace(*cp))
395 +                                break;
396 +                        continue;
397 +                case 't':                       /* object thickness */
398 +                        if (*++cp != '=')
399 +                                break;
400 +                        if (!isfltd(++cp, " \t\n\r"))
401 +                                break;
402 +                        thisillum.thick = atof(cp);
403 +                        if (thisillum.thick < .0)
404 +                                thisillum.thick = .0;
405 +                        cp = sskip(cp);
406 +                        continue;
407                  case '!':                       /* processed file! */
408                          sprintf(errmsg, "(%s): already processed!", nm);
409                          error(WARNING, errmsg);
# Line 503 | Line 449 | printopts(void)                        /* print out option default values *
449          printf("d=%d\t\t\t\t# density of points\n", thisillum.sampdens);
450          printf("s=%d\t\t\t\t# samples per point\n", thisillum.nsamps);
451          printf("b=%f\t\t\t# minimum average brightness\n", thisillum.minbrt);
452 +        switch (thisillum.udir) {
453 +        case UDzneg:
454 +                fputs("u=-Z\t\t\t\t# up is negative Z\n", stdout);
455 +                break;
456 +        case UDyneg:
457 +                fputs("u=-Y\t\t\t\t# up is negative Y\n", stdout);
458 +                break;
459 +        case UDxneg:
460 +                fputs("u=-X\t\t\t\t# up is negative X\n", stdout);
461 +                break;
462 +        case UDxpos:
463 +                fputs("u=+X\t\t\t\t# up is positive X\n", stdout);
464 +                break;
465 +        case UDypos:
466 +                fputs("u=+Y\t\t\t\t# up is positive Y\n", stdout);
467 +                break;
468 +        case UDzpos:
469 +                fputs("u=+Z\t\t\t\t# up is positive Z\n", stdout);
470 +                break;
471 +        case UDunknown:
472 +                break;
473 +        }
474 +        printf("t=%f\t\t\t# object thickness\n", thisillum.thick);
475   }
476  
477  
# Line 576 | Line 545 | xobject(                               /* translate an object from fp */
545          checkhead();
546                                                  /* process object */
547          if (doit)
548 <                (*ofun[thisobj.otype].funp)(&thisobj, &thisillum, &rt0, nm);
548 >                switch (thisobj.otype) {
549 >                case OBJ_FACE:
550 >                        my_face(&thisobj, &thisillum, nm);
551 >                        break;
552 >                case OBJ_SPHERE:
553 >                        my_sphere(&thisobj, &thisillum, nm);
554 >                        break;
555 >                case OBJ_RING:
556 >                        my_ring(&thisobj, &thisillum, nm);
557 >                        break;
558 >                default:
559 >                        my_default(&thisobj, &thisillum, nm);
560 >                        break;
561 >                }
562          else
563                  printobj(thisillum.altmat, &thisobj);
564                                                  /* free arguments */
565          freefargs(&thisobj.oargs);
566          return;
567   readerr:
568 <        sprintf(errmsg, "(%s): error reading scene", nm);
568 >        sprintf(errmsg, "(%s): error reading input", nm);
569          error(USER, errmsg);
570   }

Diff Legend

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