ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rad.c
Revision: 2.22
Committed: Wed Aug 25 13:48:37 1993 UTC (30 years, 7 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.21: +5 -3 lines
Log Message:
added -R option to rview

File Contents

# User Rev Content
1 greg 2.1 /* Copyright (c) 1993 Regents of the University of California */
2    
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * Executive program for oconv, rpict and pfilt
9     */
10    
11     #include "standard.h"
12     #include "paths.h"
13     #include <ctype.h>
14    
15    
16     typedef struct {
17     char *name; /* variable name */
18     short nick; /* # characters required for nickname */
19     short nass; /* # assignments made */
20     char *value; /* assigned value(s) */
21     int (*fixval)(); /* assignment checking function */
22     } VARIABLE;
23    
24     int onevalue(), catvalues();
25    
26     /* variables */
27     #define OBJECT 0 /* object files */
28     #define SCENE 1 /* scene files */
29     #define MATERIAL 2 /* material files */
30     #define RENDER 3 /* rendering options */
31     #define OCONV 4 /* oconv options */
32     #define PFILT 5 /* pfilt options */
33     #define VIEW 6 /* view(s) for picture(s) */
34     #define ZONE 7 /* simulation zone */
35     #define QUALITY 8 /* desired rendering quality */
36     #define OCTREE 9 /* octree file name */
37     #define PICTURE 10 /* picture file name */
38     #define AMBFILE 11 /* ambient file name */
39     #define OPTFILE 12 /* rendering options file */
40     #define EXPOSURE 13 /* picture exposure setting */
41     #define RESOLUTION 14 /* maximum picture resolution */
42     #define UP 15 /* view up (X, Y or Z) */
43     #define INDIRECT 16 /* indirection in lighting */
44     #define DETAIL 17 /* level of scene detail */
45     #define PENUMBRAS 18 /* shadow penumbras are desired */
46     #define VARIABILITY 19 /* level of light variability */
47     #define REPORT 20 /* report frequency and errfile */
48     /* total number of variables */
49     #define NVARS 21
50    
51     VARIABLE vv[NVARS] = { /* variable-value pairs */
52     {"objects", 3, 0, NULL, catvalues},
53     {"scene", 3, 0, NULL, catvalues},
54     {"materials", 3, 0, NULL, catvalues},
55     {"render", 3, 0, NULL, catvalues},
56     {"oconv", 3, 0, NULL, catvalues},
57     {"pfilt", 2, 0, NULL, catvalues},
58     {"view", 2, 0, NULL, NULL},
59     {"ZONE", 2, 0, NULL, onevalue},
60     {"QUALITY", 3, 0, NULL, onevalue},
61     {"OCTREE", 3, 0, NULL, onevalue},
62     {"PICTURE", 3, 0, NULL, onevalue},
63     {"AMBFILE", 3, 0, NULL, onevalue},
64     {"OPTFILE", 3, 0, NULL, onevalue},
65     {"EXPOSURE", 3, 0, NULL, onevalue},
66     {"RESOLUTION", 3, 0, NULL, onevalue},
67     {"UP", 2, 0, NULL, onevalue},
68     {"INDIRECT", 3, 0, NULL, onevalue},
69     {"DETAIL", 3, 0, NULL, onevalue},
70     {"PENUMBRAS", 3, 0, NULL, onevalue},
71     {"VARIABILITY", 3, 0, NULL, onevalue},
72     {"REPORT", 3, 0, NULL, onevalue},
73     };
74    
75     VARIABLE *matchvar();
76     char *nvalue();
77     int vscale();
78    
79 greg 2.3 #define UPPER(c) ((c)&~0x20) /* ASCII trick */
80    
81 greg 2.1 #define vnam(vc) (vv[vc].name)
82     #define vdef(vc) (vv[vc].nass)
83     #define vval(vc) (vv[vc].value)
84     #define vint(vc) atoi(vval(vc))
85 greg 2.3 #define vlet(vc) UPPER(vval(vc)[0])
86 greg 2.1 #define vbool(vc) (vlet(vc)=='T')
87    
88     #define HIGH 2
89     #define MEDIUM 1
90     #define LOW 0
91    
92     int lowqopts(), medqopts(), hiqopts();
93     int (*setqopts[3])() = {lowqopts, medqopts, hiqopts};
94    
95     #define renderopts (*setqopts[vscale(QUALITY)])
96    
97 greg 2.6 /* overture calculation file */
98     #ifdef NIX
99     char overfile[] = "overture.raw";
100     #else
101     char overfile[] = "/dev/null";
102     #endif
103    
104 greg 2.2 extern long fdate(), time();
105 greg 2.1
106     long scenedate; /* date of latest scene or object file */
107     long octreedate; /* date of octree */
108 greg 2.11 long matdate; /* date of latest material file */
109 greg 2.1
110     int explicate = 0; /* explicate variables */
111     int silent = 0; /* do work silently */
112     int noaction = 0; /* don't do anything */
113 greg 2.21 int sayview = 0; /* print view out */
114 greg 2.1 char *rvdevice = NULL; /* rview output device */
115     char *viewselect = NULL; /* specific view only */
116    
117     int overture = 0; /* overture calculation needed */
118    
119     char *progname; /* global argv[0] */
120 greg 2.22 char *rifname; /* global rad input file name */
121 greg 2.1
122     char radname[MAXPATH]; /* root Radiance file name */
123    
124    
125     main(argc, argv)
126     int argc;
127     char *argv[];
128     {
129     char ropts[512];
130     int i;
131    
132     progname = argv[0];
133     /* get options */
134     for (i = 1; i < argc && argv[i][0] == '-'; i++)
135     switch (argv[i][1]) {
136     case 's':
137     silent++;
138     break;
139     case 'n':
140     noaction++;
141     break;
142     case 'e':
143     explicate++;
144     break;
145     case 'o':
146     rvdevice = argv[++i];
147     break;
148 greg 2.20 case 'V':
149 greg 2.21 sayview++;
150     break;
151 greg 2.1 case 'v':
152     viewselect = argv[++i];
153     break;
154     default:
155     goto userr;
156     }
157     if (i >= argc)
158     goto userr;
159 greg 2.22 rifname = argv[i];
160 greg 2.1 /* assign Radiance root file name */
161 greg 2.22 rootname(radname, rifname);
162 greg 2.1 /* load variable values */
163 greg 2.22 load(rifname);
164 greg 2.1 /* get any additional assignments */
165     for (i++; i < argc; i++)
166     setvariable(argv[i]);
167     /* check assignments */
168     checkvalues();
169     /* check files and dates */
170     checkfiles();
171     /* set default values as necessary */
172     setdefaults();
173     /* print all values if requested */
174     if (explicate)
175     printvals();
176 greg 2.7 /* build octree */
177     oconv();
178     /* check date on ambient file */
179     checkambfile();
180 greg 2.1 /* run simulation */
181     renderopts(ropts);
182     xferopts(ropts);
183     if (rvdevice != NULL)
184     rview(ropts);
185     else
186     rpict(ropts);
187     exit(0);
188     userr:
189     fprintf(stderr,
190 greg 2.3 "Usage: %s [-s][-n][-e][-v view][-o dev] rfile [VAR=value ..]\n",
191 greg 2.1 progname);
192     exit(1);
193     }
194    
195    
196     rootname(rn, fn) /* remove tail from end of fn */
197     register char *rn, *fn;
198     {
199     char *tp, *dp;
200    
201     for (tp = NULL, dp = rn; *rn = *fn++; rn++)
202     if (ISDIRSEP(*rn))
203     dp = rn;
204     else if (*rn == '.')
205     tp = rn;
206     if (tp != NULL && tp > dp)
207     *tp = '\0';
208     }
209    
210    
211     load(rfname) /* load Radiance simulation file */
212     char *rfname;
213     {
214     FILE *fp;
215 greg 2.4 char buf[512];
216 greg 2.1 register char *cp;
217    
218     if (rfname == NULL)
219     fp = stdin;
220 greg 2.3 else if ((fp = fopen(rfname, "r")) == NULL)
221     syserr(rfname);
222 greg 2.1 while (fgetline(buf, sizeof(buf), fp) != NULL) {
223     for (cp = buf; *cp; cp++) {
224     switch (*cp) {
225     case '\\':
226     case '\n':
227     *cp = ' ';
228     continue;
229     case '#':
230     *cp = '\0';
231     break;
232 greg 2.4 default:
233     continue;
234 greg 2.1 }
235     break;
236     }
237     setvariable(buf);
238     }
239     fclose(fp);
240     }
241    
242    
243     setvariable(ass) /* assign variable according to string */
244     register char *ass;
245     {
246     char varname[32];
247     register char *cp;
248     register VARIABLE *vp;
249     register int i;
250     int n;
251    
252     while (isspace(*ass)) /* skip leading space */
253     ass++;
254     cp = varname; /* extract name */
255     while (cp < varname+sizeof(varname)-1
256     && *ass && !isspace(*ass) && *ass != '=')
257     *cp++ = *ass++;
258     *cp = '\0';
259     if (!varname[0])
260     return; /* no variable name! */
261     /* trim value */
262     while (isspace(*ass) || *ass == '=')
263     ass++;
264     cp = ass + strlen(ass);
265     do
266     *cp-- = '\0';
267     while (cp >= ass && isspace(*cp));
268     n = cp - ass + 1;
269     if (!n) {
270     fprintf(stderr, "%s: warning - missing value for variable '%s'\n",
271     progname, varname);
272     return;
273     }
274     /* match variable from list */
275     vp = matchvar(varname);
276     if (vp == NULL) {
277     fprintf(stderr, "%s: unknown variable '%s'\n",
278     progname, varname);
279     exit(1);
280     }
281     /* assign new value */
282 greg 2.4 if (i = vp->nass) {
283     cp = vp->value;
284     while (i--)
285     while (*cp++)
286     ;
287     i = cp - vp->value;
288     vp->value = realloc(vp->value, i+n+1);
289     } else
290     vp->value = malloc(n+1);
291 greg 2.3 if (vp->value == NULL)
292     syserr(progname);
293 greg 2.1 strcpy(vp->value+i, ass);
294     vp->nass++;
295     }
296    
297    
298     VARIABLE *
299     matchvar(nam) /* match a variable by its name */
300     char *nam;
301     {
302     int n = strlen(nam);
303     register int i;
304    
305     for (i = 0; i < NVARS; i++)
306     if (n >= vv[i].nick && !strncmp(nam, vv[i].name, n))
307     return(vv+i);
308     return(NULL);
309     }
310    
311    
312     char *
313     nvalue(vp, n) /* return nth variable value */
314     VARIABLE *vp;
315     register int n;
316     {
317     register char *cp;
318    
319 greg 2.19 if (vp == NULL | n < 0 | n >= vp->nass)
320 greg 2.1 return(NULL);
321     cp = vp->value;
322     while (n--)
323     while (*cp++)
324     ;
325     return(cp);
326     }
327    
328    
329     int
330     vscale(vc) /* return scale for variable vc */
331     int vc;
332     {
333     switch(vlet(vc)) {
334     case 'H':
335     return(HIGH);
336     case 'M':
337     return(MEDIUM);
338     case 'L':
339     return(LOW);
340     }
341 greg 2.3 badvalue(vc);
342 greg 2.1 }
343    
344    
345     checkvalues() /* check assignments */
346     {
347     register int i;
348    
349     for (i = 0; i < NVARS; i++)
350     if (vv[i].fixval != NULL)
351     (*vv[i].fixval)(vv+i);
352     }
353    
354    
355     onevalue(vp) /* only one assignment for this variable */
356     register VARIABLE *vp;
357     {
358     if (vp->nass < 2)
359     return;
360     fprintf(stderr, "%s: warning - multiple assignment of variable '%s'\n",
361     progname, vp->name);
362 greg 2.2 do
363 greg 2.1 vp->value += strlen(vp->value)+1;
364 greg 2.2 while (--vp->nass > 1);
365 greg 2.1 }
366    
367    
368     catvalues(vp) /* concatenate variable values */
369     register VARIABLE *vp;
370     {
371     register char *cp;
372    
373     if (vp->nass < 2)
374     return;
375     for (cp = vp->value; vp->nass > 1; vp->nass--) {
376     while (*cp)
377     cp++;
378     *cp++ = ' ';
379     }
380     }
381    
382    
383     long
384     checklast(fnames) /* check files and find most recent */
385     register char *fnames;
386     {
387     char thisfile[MAXPATH];
388     long thisdate, lastdate = -1;
389     register char *cp;
390    
391     while (*fnames) {
392     while (isspace(*fnames)) fnames++;
393     cp = thisfile;
394 greg 2.4 while (*fnames && !isspace(*fnames))
395     *cp++ = *fnames++;
396 greg 2.1 *cp = '\0';
397 greg 2.3 if ((thisdate = fdate(thisfile)) < 0)
398     syserr(thisfile);
399 greg 2.1 if (thisdate > lastdate)
400     lastdate = thisdate;
401     }
402     return(lastdate);
403     }
404    
405    
406     checkfiles() /* check for existence and modified times */
407     {
408 greg 2.2 char *cp;
409 greg 2.1 long objdate;
410    
411 greg 2.2 if (!vdef(OCTREE)) {
412 greg 2.3 if ((cp = bmalloc(strlen(radname)+5)) == NULL)
413     syserr(progname);
414 greg 2.2 sprintf(cp, "%s.oct", radname);
415     vval(OCTREE) = cp;
416     vdef(OCTREE)++;
417     }
418     octreedate = fdate(vval(OCTREE));
419 greg 2.1 scenedate = -1;
420     if (vdef(SCENE)) {
421     scenedate = checklast(vval(SCENE));
422     if (vdef(OBJECT)) {
423     objdate = checklast(vval(OBJECT));
424     if (objdate > scenedate)
425     scenedate = objdate;
426     }
427     }
428     if (octreedate < 0 & scenedate < 0) {
429 greg 2.2 fprintf(stderr, "%s: need '%s' or '%s'\n", progname,
430     vnam(OCTREE), vnam(SCENE));
431 greg 2.1 exit(1);
432     }
433 greg 2.11 matdate = -1;
434     if (vdef(MATERIAL))
435     matdate = checklast(vval(MATERIAL));
436 greg 2.1 }
437    
438    
439 greg 2.2 getoctcube(org, sizp) /* get octree bounding cube */
440     double org[3], *sizp;
441 greg 2.1 {
442 greg 2.2 extern FILE *popen();
443     static double oorg[3], osiz = 0.;
444     char buf[MAXPATH+16];
445 greg 2.1 FILE *fp;
446    
447 greg 2.2 if (osiz <= FTINY) {
448     oconv(); /* does nothing if done already */
449     sprintf(buf, "getinfo -d < %s", vval(OCTREE));
450 greg 2.3 if ((fp = popen(buf, "r")) == NULL)
451     syserr("getinfo");
452 greg 2.2 if (fscanf(fp, "%lf %lf %lf %lf", &oorg[0], &oorg[1],
453     &oorg[2], &osiz) != 4) {
454     fprintf(stderr,
455     "%s: error reading bounding cube from getinfo\n",
456     progname);
457     exit(1);
458     }
459     pclose(fp);
460 greg 2.1 }
461 greg 2.2 org[0] = oorg[0]; org[1] = oorg[1]; org[2] = oorg[2]; *sizp = osiz;
462     }
463    
464    
465     setdefaults() /* set default values for unassigned var's */
466     {
467     double org[3], size;
468     char buf[128];
469    
470 greg 2.1 if (!vdef(ZONE)) {
471 greg 2.2 getoctcube(org, &size);
472     sprintf(buf, "E %g %g %g %g %g %g", org[0], org[0]+size,
473     org[1], org[1]+size, org[2], org[2]+size);
474 greg 2.1 vval(ZONE) = savqstr(buf);
475     vdef(ZONE)++;
476     }
477     if (!vdef(INDIRECT)) {
478     vval(INDIRECT) = "0";
479     vdef(INDIRECT)++;
480     }
481     if (!vdef(QUALITY)) {
482     vval(QUALITY) = "L";
483     vdef(QUALITY)++;
484     }
485     if (!vdef(RESOLUTION)) {
486     vval(RESOLUTION) = "512";
487     vdef(RESOLUTION)++;
488     }
489     if (!vdef(PICTURE)) {
490     vval(PICTURE) = radname;
491     vdef(PICTURE)++;
492     }
493     if (!vdef(VIEW)) {
494     vval(VIEW) = "X";
495     vdef(VIEW)++;
496     }
497     if (!vdef(DETAIL)) {
498     vval(DETAIL) = "M";
499     vdef(DETAIL)++;
500     }
501     if (!vdef(PENUMBRAS)) {
502     vval(PENUMBRAS) = "F";
503     vdef(PENUMBRAS)++;
504     }
505     if (!vdef(VARIABILITY)) {
506     vval(VARIABILITY) = "L";
507     vdef(VARIABILITY)++;
508     }
509     }
510    
511    
512     printvals() /* print variable values */
513     {
514     register int i, j;
515    
516     for (i = 0; i < NVARS; i++)
517 greg 2.4 for (j = 0; j < vdef(i); j++)
518     printf("%s= %s\n", vnam(i), nvalue(vv+i, j));
519 greg 2.1 fflush(stdout);
520     }
521    
522    
523     oconv() /* run oconv if necessary */
524     {
525     char combuf[512], ocopts[64];
526    
527 greg 2.2 if (octreedate >= scenedate) /* check dates */
528 greg 2.1 return;
529     /* build command */
530     oconvopts(ocopts);
531 greg 2.4 if (vdef(MATERIAL))
532     sprintf(combuf, "oconv%s %s %s > %s", ocopts,
533     vval(MATERIAL), vval(SCENE), vval(OCTREE));
534     else
535     sprintf(combuf, "oconv%s %s > %s", ocopts,
536     vval(SCENE), vval(OCTREE));
537    
538     if (runcom(combuf)) { /* run it */
539 greg 2.1 fprintf(stderr, "%s: error generating octree\n\t%s removed\n",
540     progname, vval(OCTREE));
541     unlink(vval(OCTREE));
542     exit(1);
543     }
544 greg 2.2 octreedate = time(0);
545 greg 2.1 }
546    
547    
548     char *
549     addarg(op, arg) /* add argument and advance pointer */
550     register char *op, *arg;
551     {
552     *op = ' ';
553     while (*++op = *arg++)
554     ;
555     return(op);
556     }
557    
558    
559     oconvopts(oo) /* get oconv options */
560 greg 2.2 register char *oo;
561 greg 2.1 {
562 greg 2.2 /* BEWARE: This may be called via setdefaults(), so no assumptions */
563    
564 greg 2.1 *oo = '\0';
565     if (vdef(OCONV))
566     addarg(oo, vval(OCONV));
567     }
568    
569    
570 greg 2.7 checkambfile() /* check date on ambient file */
571     {
572     long afdate;
573    
574 greg 2.11 if (!vdef(AMBFILE))
575     return;
576     if ((afdate = fdate(vval(AMBFILE))) < 0)
577     return;
578     if (octreedate > afdate | matdate > afdate)
579     rmfile(vval(AMBFILE));
580 greg 2.7 }
581    
582    
583 greg 2.2 double
584     ambval() /* compute ambient value */
585 greg 2.1 {
586 greg 2.3 if (vdef(EXPOSURE)) {
587 greg 2.20 if (!isflt(vval(EXPOSURE)))
588     badvalue(EXPOSURE);
589 greg 2.2 if (vval(EXPOSURE)[0] == '+' || vval(EXPOSURE)[0] == '-')
590     return(.5/pow(2.,atof(vval(EXPOSURE))));
591 greg 2.20 return(.5/atof(vval(EXPOSURE)));
592 greg 2.3 }
593 greg 2.2 if (vlet(ZONE) == 'E')
594     return(10.);
595 greg 2.3 if (vlet(ZONE) == 'I')
596 greg 2.2 return(.01);
597 greg 2.3 badvalue(ZONE);
598 greg 2.2 }
599 greg 2.1
600 greg 2.2
601     lowqopts(op) /* low quality rendering options */
602     register char *op;
603     {
604     double d, org[3], siz[3];
605    
606 greg 2.1 *op = '\0';
607 greg 2.2 if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
608 greg 2.3 &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
609     badvalue(ZONE);
610 greg 2.2 siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
611     getoctcube(org, &d);
612     d *= 3./(siz[0]+siz[1]+siz[2]);
613     switch (vscale(DETAIL)) {
614     case LOW:
615 greg 2.4 op = addarg(op, "-ps 16 -dp 16");
616 greg 2.2 sprintf(op, " -ar %d", (int)(4*d));
617     op += strlen(op);
618     break;
619     case MEDIUM:
620 greg 2.4 op = addarg(op, "-ps 8 -dp 32");
621 greg 2.2 sprintf(op, " -ar %d", (int)(8*d));
622     op += strlen(op);
623     break;
624     case HIGH:
625 greg 2.4 op = addarg(op, "-ps 4 -dp 64");
626 greg 2.2 sprintf(op, " -ar %d", (int)(16*d));
627     op += strlen(op);
628     break;
629     }
630     op = addarg(op, "-pt .16");
631     if (vbool(PENUMBRAS))
632     op = addarg(op, "-ds .4");
633 greg 2.4 else
634     op = addarg(op, "-ds 0");
635 greg 2.17 op = addarg(op, "-dt .2 -dc .25 -dr 0 -sj 0 -st .5");
636 greg 2.2 if (vdef(AMBFILE)) {
637     sprintf(op, " -af %s", vval(AMBFILE));
638     op += strlen(op);
639     } else
640     overture = 0;
641     switch (vscale(VARIABILITY)) {
642     case LOW:
643 greg 2.16 op = addarg(op, "-aa .4 -ad 64");
644 greg 2.2 break;
645     case MEDIUM:
646 greg 2.16 op = addarg(op, "-aa .3 -ad 128");
647 greg 2.2 break;
648     case HIGH:
649 greg 2.16 op = addarg(op, "-aa .25 -ad 256");
650 greg 2.2 break;
651     }
652     op = addarg(op, "-as 0");
653     d = ambval();
654     sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
655     op += strlen(op);
656 greg 2.4 op = addarg(op, "-lr 3 -lw .02");
657 greg 2.1 if (vdef(RENDER))
658     op = addarg(op, vval(RENDER));
659     }
660    
661    
662 greg 2.2 medqopts(op) /* medium quality rendering options */
663     register char *op;
664 greg 2.1 {
665 greg 2.2 double d, org[3], siz[3];
666 greg 2.1
667     *op = '\0';
668 greg 2.2 if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
669 greg 2.3 &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
670     badvalue(ZONE);
671 greg 2.2 siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
672     getoctcube(org, &d);
673     d *= 3./(siz[0]+siz[1]+siz[2]);
674     switch (vscale(DETAIL)) {
675     case LOW:
676     op = addarg(op, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
677     op = addarg(op, "-dp 64");
678     sprintf(op, " -ar %d", (int)(8*d));
679     op += strlen(op);
680     break;
681     case MEDIUM:
682     op = addarg(op, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
683     op = addarg(op, "-dp 128");
684     sprintf(op, " -ar %d", (int)(16*d));
685     op += strlen(op);
686     break;
687     case HIGH:
688     op = addarg(op, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
689     op = addarg(op, "-dp 256");
690     sprintf(op, " -ar %d", (int)(32*d));
691     op += strlen(op);
692     break;
693     }
694     op = addarg(op, "-pt .08");
695 greg 2.4 if (vbool(PENUMBRAS))
696     op = addarg(op, "-ds .2 -dj .35");
697     else
698 greg 2.2 op = addarg(op, "-ds .3");
699 greg 2.17 op = addarg(op, "-dt .1 -dc .5 -dr 1 -sj .7 -st .1");
700 greg 2.5 if (overture = vint(INDIRECT)) {
701     sprintf(op, " -ab %d", overture);
702     op += strlen(op);
703     }
704 greg 2.2 if (vdef(AMBFILE)) {
705     sprintf(op, " -af %s", vval(AMBFILE));
706     op += strlen(op);
707     } else
708     overture = 0;
709     switch (vscale(VARIABILITY)) {
710     case LOW:
711 greg 2.16 op = addarg(op, "-aa .25 -ad 196 -as 0");
712 greg 2.2 break;
713     case MEDIUM:
714 greg 2.16 op = addarg(op, "-aa .2 -ad 400 -as 64");
715 greg 2.2 break;
716     case HIGH:
717 greg 2.16 op = addarg(op, "-aa .15 -ad 768 -as 196");
718 greg 2.2 break;
719     }
720     d = ambval();
721     sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
722     op += strlen(op);
723 greg 2.4 op = addarg(op, "-lr 6 -lw .002");
724 greg 2.1 if (vdef(RENDER))
725     op = addarg(op, vval(RENDER));
726     }
727    
728    
729 greg 2.2 hiqopts(op) /* high quality rendering options */
730     register char *op;
731 greg 2.1 {
732 greg 2.2 double d, org[3], siz[3];
733 greg 2.1
734     *op = '\0';
735 greg 2.2 if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
736 greg 2.3 &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
737     badvalue(ZONE);
738 greg 2.2 siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
739     getoctcube(org, &d);
740     d *= 3./(siz[0]+siz[1]+siz[2]);
741     switch (vscale(DETAIL)) {
742     case LOW:
743     op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
744     op = addarg(op, "-dp 256");
745     sprintf(op, " -ar %d", (int)(16*d));
746     op += strlen(op);
747     break;
748     case MEDIUM:
749     op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
750     op = addarg(op, "-dp 512");
751     sprintf(op, " -ar %d", (int)(32*d));
752     op += strlen(op);
753     break;
754     case HIGH:
755     op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
756     op = addarg(op, "-dp 1024");
757     sprintf(op, " -ar %d", (int)(64*d));
758     op += strlen(op);
759     break;
760     }
761     op = addarg(op, "-pt .04");
762 greg 2.4 if (vbool(PENUMBRAS))
763     op = addarg(op, "-ds .1 -dj .7");
764     else
765 greg 2.2 op = addarg(op, "-ds .2");
766 greg 2.17 op = addarg(op, "-dt .05 -dc .75 -dr 3 -sj 1 -st .01");
767 greg 2.2 sprintf(op, " -ab %d", overture=vint(INDIRECT)+1);
768     op += strlen(op);
769     if (vdef(AMBFILE)) {
770     sprintf(op, " -af %s", vval(AMBFILE));
771     op += strlen(op);
772     } else
773     overture = 0;
774     switch (vscale(VARIABILITY)) {
775     case LOW:
776 greg 2.16 op = addarg(op, "-aa .15 -ad 256 -as 0");
777 greg 2.2 break;
778     case MEDIUM:
779 greg 2.16 op = addarg(op, "-aa .125 -ad 512 -as 256");
780 greg 2.2 break;
781     case HIGH:
782 greg 2.16 op = addarg(op, "-aa .08 -ad 1024 -as 512");
783 greg 2.2 break;
784     }
785     d = ambval();
786     sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
787     op += strlen(op);
788 greg 2.4 op = addarg(op, "-lr 12 -lw .0005");
789 greg 2.1 if (vdef(RENDER))
790     op = addarg(op, vval(RENDER));
791     }
792    
793    
794     xferopts(ro) /* transfer options if indicated */
795     char *ro;
796     {
797     int fd, n;
798 greg 2.10 register char *cp;
799 greg 2.1
800     n = strlen(ro);
801     if (n < 2)
802     return;
803     if (vdef(OPTFILE)) {
804 greg 2.10 for (cp = ro; cp[1]; cp++)
805     if (isspace(cp[1]) && cp[2] == '-' && isalpha(cp[3]))
806     *cp = '\n';
807     else
808     *cp = cp[1];
809     *cp = '\n';
810     fd = open(vval(OPTFILE), O_WRONLY|O_CREAT|O_TRUNC, 0666);
811     if (fd < 0 || write(fd, ro, n) != n || close(fd) < 0)
812 greg 2.3 syserr(vval(OPTFILE));
813 greg 2.15 sprintf(ro, " @%s", vval(OPTFILE));
814 greg 2.1 }
815     #ifdef MSDOS
816     else if (n > 50) {
817 greg 2.8 setenv("ROPT", ro+1);
818 greg 2.1 strcpy(ro, " $ROPT");
819     }
820     #endif
821     }
822    
823    
824     pfiltopts(po) /* get pfilt options */
825     register char *po;
826     {
827     *po = '\0';
828     if (vdef(EXPOSURE)) {
829     po = addarg(po, "-1 -e");
830     po = addarg(po, vval(EXPOSURE));
831     }
832 greg 2.12 switch (vscale(QUALITY)) {
833     case MEDIUM:
834     po = addarg(po, "-r 1");
835     break;
836     case HIGH:
837 greg 2.14 po = addarg(po, "-m .25");
838 greg 2.12 break;
839     }
840 greg 2.1 if (vdef(PFILT))
841     po = addarg(po, vval(PFILT));
842     }
843    
844    
845     matchword(s1, s2) /* match white-delimited words */
846     register char *s1, *s2;
847     {
848     while (isspace(*s1)) s1++;
849     while (isspace(*s2)) s2++;
850     while (*s1 && !isspace(*s1))
851     if (*s1++ != *s2++)
852     return(0);
853     return(!*s2 || isspace(*s2));
854     }
855    
856    
857     char *
858     specview(vs) /* get proper view spec from vs */
859     register char *vs;
860     {
861 greg 2.3 static char vup[7][12] = {"-vu 0 0 -1","-vu 0 -1 0","-vu -1 0 0",
862     "-vu 0 0 1", "-vu 1 0 0","-vu 0 1 0","-vu 0 0 1"};
863 greg 2.1 static char viewopts[128];
864     register char *cp;
865 greg 2.3 int xpos, ypos, zpos, viewtype, upax;
866     register int i;
867 greg 2.1 double cent[3], dim[3], mult, d;
868    
869 greg 2.9 if (vs == NULL || *vs == '-')
870 greg 2.1 return(vs);
871 greg 2.3 upax = 0; /* get the up vector */
872     if (vdef(UP)) {
873     if (vval(UP)[0] == '-' || vval(UP)[0] == '+')
874     upax = 1-'X'+UPPER(vval(UP)[1]);
875     else
876     upax = 1-'X'+vlet(UP);
877     if (upax < 1 | upax > 3)
878     badvalue(UP);
879     if (vval(UP)[0] == '-')
880     upax = -upax;
881     }
882 greg 2.1 /* check standard view names */
883     xpos = ypos = zpos = 0;
884     if (*vs == 'X') {
885     xpos = 1; vs++;
886     } else if (*vs == 'x') {
887     xpos = -1; vs++;
888     }
889     if (*vs == 'Y') {
890     ypos = 1; vs++;
891     } else if (*vs == 'y') {
892     ypos = -1; vs++;
893     }
894     if (*vs == 'Z') {
895     zpos = 1; vs++;
896     } else if (*vs == 'z') {
897     zpos = -1; vs++;
898     }
899 greg 2.3 viewtype = 'v';
900 greg 2.1 if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h')
901     viewtype = *vs++;
902     cp = viewopts;
903 greg 2.3 if ((!*vs || isspace(*vs)) && (xpos|ypos|zpos)) { /* got one! */
904 greg 2.1 *cp++ = '-'; *cp++ = 'v'; *cp++ = 't'; *cp++ = viewtype;
905     if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf",
906     &cent[0], &dim[0], &cent[1], &dim[1],
907 greg 2.3 &cent[2], &dim[2]) != 6)
908     badvalue(ZONE);
909     for (i = 0; i < 3; i++) {
910     dim[i] -= cent[i];
911     cent[i] += .5*dim[i];
912 greg 2.1 }
913 greg 2.3 mult = vlet(ZONE)=='E' ? 2. : .45 ;
914 greg 2.1 sprintf(cp, " -vp %.2g %.2g %.2g -vd %.2g %.2g %.2g",
915     cent[0]+xpos*mult*dim[0],
916     cent[1]+ypos*mult*dim[1],
917     cent[2]+zpos*mult*dim[2],
918     -xpos*dim[0], -ypos*dim[1], -zpos*dim[2]);
919     cp += strlen(cp);
920 greg 2.3 /* redirect up axis if necessary */
921     switch (upax) {
922     case 3: /* plus or minus Z axis */
923     case -3:
924     case 0:
925     if (!(xpos|ypos))
926     upax = 2;
927 greg 2.1 break;
928 greg 2.3 case 2: /* plus or minus Y axis */
929     case -2:
930     if (!(xpos|zpos))
931     upax = 1;
932     break;
933     case 1: /* plus or minus X axis */
934     case -1:
935     if (!(ypos|zpos))
936     upax = 3;
937     break;
938 greg 2.1 }
939 greg 2.3 cp = addarg(cp, vup[upax+3]);
940 greg 2.1 switch (viewtype) {
941     case 'v':
942     cp = addarg(cp, "-vh 45 -vv 45");
943     break;
944     case 'l':
945     d = sqrt(dim[0]*dim[0]+dim[1]*dim[1]+dim[2]*dim[2]);
946     sprintf(cp, " -vh %.2g -vv %.2g", d, d);
947     cp += strlen(cp);
948     break;
949     case 'a':
950     case 'h':
951     cp = addarg(cp, "-vh 180 -vv 180");
952     break;
953     }
954 greg 2.3 } else {
955 greg 2.4 while (!isspace(*vs)) /* else skip id */
956     if (!*vs++)
957     return(NULL);
958 greg 2.3 if (upax) { /* specify up vector */
959     strcpy(cp, vup[upax+3]);
960     cp += strlen(cp);
961     }
962     }
963 greg 2.8 strcpy(cp, vs); /* append any additional options */
964     #ifdef MSDOS
965     if (strlen(viewopts) > 40) {
966     setenv("VIEW", viewopts);
967     return("$VIEW");
968     }
969     #endif
970 greg 2.1 return(viewopts);
971     }
972    
973    
974     char *
975     getview(n, vn) /* get view n, or NULL if none */
976     int n;
977 greg 2.8 char *vn; /* returned view name */
978 greg 2.1 {
979 greg 2.19 register char *mv;
980 greg 2.1
981 greg 2.8 if (viewselect != NULL) { /* command-line selected */
982 greg 2.1 if (n) /* only do one */
983     return(NULL);
984     if (viewselect[0] == '-') { /* already specified */
985     if (vn != NULL) *vn = '\0';
986     return(viewselect);
987     }
988     if (vn != NULL) {
989     for (mv = viewselect; *mv && !isspace(*mv);
990     *vn++ = *mv++)
991     ;
992     *vn = '\0';
993     }
994 greg 2.19 /* view number? */
995 greg 2.20 if (isint(viewselect))
996 greg 2.19 return(specview(nvalue(vv+VIEW, atoi(viewselect)-1)));
997 greg 2.1 /* check list */
998     while ((mv = nvalue(vv+VIEW, n++)) != NULL)
999     if (matchword(viewselect, mv))
1000     return(specview(mv));
1001     return(specview(viewselect)); /* standard view? */
1002     }
1003 greg 2.8 mv = nvalue(vv+VIEW, n); /* use view n */
1004     if (vn != NULL & mv != NULL) {
1005 greg 2.19 register char *mv2 = mv;
1006     if (*mv2 != '-')
1007     while (*mv2 && !isspace(*mv2))
1008     *vn++ = *mv2++;
1009 greg 2.1 *vn = '\0';
1010     }
1011 greg 2.8 return(specview(mv));
1012 greg 2.20 }
1013    
1014    
1015 greg 2.21 printview(vopts) /* print out selected view */
1016     register char *vopts;
1017 greg 2.20 {
1018     extern char *atos();
1019     char buf[256];
1020     FILE *fp;
1021 greg 2.21 register char *cp;
1022 greg 2.20
1023     if (vopts == NULL)
1024     return(-1);
1025     fputs("VIEW=", stdout);
1026     do {
1027     if (matchword(vopts, "-vf")) { /* expand view file */
1028     vopts = sskip(vopts);
1029     if (!*atos(buf, sizeof(buf), vopts))
1030     return(-1);
1031     if ((fp = fopen(buf, "r")) == NULL)
1032     return(-1);
1033     for (buf[sizeof(buf)-2] = '\n';
1034     fgets(buf, sizeof(buf), fp) != NULL &&
1035     buf[0] != '\n';
1036     buf[sizeof(buf)-2] = '\n') {
1037     if (buf[sizeof(buf)-2] != '\n') {
1038     ungetc(buf[sizeof(buf)-2], fp);
1039     buf[sizeof(buf)-2] = '\0';
1040     }
1041     if (matchword(buf, "VIEW=") ||
1042     matchword(buf, "rview")) {
1043     for (cp = sskip(buf); *cp && *cp != '\n'; cp++)
1044     putchar(*cp);
1045     }
1046     }
1047     fclose(fp);
1048     vopts = sskip(vopts);
1049     } else {
1050     while (isspace(*vopts))
1051     vopts++;
1052     putchar(' ');
1053     while (*vopts && !isspace(*vopts))
1054     putchar(*vopts++);
1055     }
1056     } while (*vopts++);
1057     putchar('\n');
1058     return(0);
1059 greg 2.1 }
1060    
1061    
1062     rview(opts) /* run rview with first view */
1063     char *opts;
1064     {
1065 greg 2.21 char *vw;
1066 greg 2.1 char combuf[512];
1067     /* build command */
1068 greg 2.21 if ((vw = getview(0, NULL)) == NULL)
1069     return;
1070     if (sayview)
1071     printview(vw);
1072 greg 2.22 sprintf(combuf, "rview %s%s -R %s ", vw, opts, rifname);
1073 greg 2.1 if (rvdevice != NULL)
1074     sprintf(combuf+strlen(combuf), "-o %s ", rvdevice);
1075     strcat(combuf, vval(OCTREE));
1076 greg 2.4 if (runcom(combuf)) { /* run it */
1077 greg 2.1 fprintf(stderr, "%s: error running rview\n", progname);
1078     exit(1);
1079     }
1080     }
1081    
1082    
1083     rpict(opts) /* run rpict and pfilt for each view */
1084     char *opts;
1085     {
1086 greg 2.2 char combuf[1024];
1087 greg 2.3 char rawfile[MAXPATH], picfile[MAXPATH], rep[MAXPATH+16], res[32];
1088     char pfopts[128];
1089 greg 2.1 char vs[32], *vw;
1090     int vn, mult;
1091 greg 2.13 long lastdate;
1092 greg 2.1 /* get pfilt options */
1093     pfiltopts(pfopts);
1094     /* get resolution, reporting */
1095     mult = vscale(QUALITY)+1;
1096     {
1097     int xres, yres;
1098     double aspect;
1099     int n;
1100     n = sscanf(vval(RESOLUTION), "%d %d %lf", &xres, &yres, &aspect);
1101     if (n == 3)
1102     sprintf(res, "-x %d -y %d -pa %.3f",
1103     mult*xres, mult*yres, aspect);
1104     else if (n) {
1105     if (n == 1) yres = xres;
1106     sprintf(res, "-x %d -y %d", mult*xres, mult*yres);
1107 greg 2.3 } else
1108     badvalue(RESOLUTION);
1109 greg 2.1 }
1110     rep[0] = '\0';
1111     if (vdef(REPORT)) {
1112     double minutes;
1113     int n;
1114     n = sscanf(vval(REPORT), "%lf %s", &minutes, rawfile);
1115     if (n == 2)
1116     sprintf(rep, " -t %d -e %s", (int)(minutes*60), rawfile);
1117     else if (n == 1)
1118     sprintf(rep, " -t %d", (int)(minutes*60));
1119 greg 2.3 else
1120     badvalue(REPORT);
1121 greg 2.1 }
1122 greg 2.13 /* get update time */
1123     lastdate = octreedate > matdate ? octreedate : matdate;
1124 greg 2.1 /* do each view */
1125     vn = 0;
1126     while ((vw = getview(vn++, vs)) != NULL) {
1127 greg 2.21 if (sayview)
1128     printview(vw);
1129 greg 2.1 if (!vs[0])
1130     sprintf(vs, "%d", vn);
1131     sprintf(picfile, "%s_%s.pic", vval(PICTURE), vs);
1132     /* check date on picture */
1133 greg 2.13 if (fdate(picfile) >= lastdate)
1134 greg 2.1 continue;
1135     /* build rpict command */
1136     sprintf(rawfile, "%s_%s.raw", vval(PICTURE), vs);
1137 greg 2.11 if (fdate(rawfile) >= octreedate) /* recover */
1138 greg 2.1 sprintf(combuf, "rpict%s%s -ro %s %s",
1139     rep, opts, rawfile, vval(OCTREE));
1140     else {
1141     if (overture) { /* run overture calculation */
1142     sprintf(combuf,
1143     "rpict%s %s%s -x 64 -y 64 -ps 1 %s > %s",
1144     rep, vw, opts,
1145 greg 2.6 vval(OCTREE), overfile);
1146 greg 2.4 if (runcom(combuf)) {
1147 greg 2.1 fprintf(stderr,
1148 greg 2.6 "%s: error in overture for view %s\n",
1149     progname, vs);
1150 greg 2.1 exit(1);
1151     }
1152 greg 2.6 #ifdef NIX
1153     rmfile(overfile);
1154     #endif
1155 greg 2.1 }
1156     sprintf(combuf, "rpict%s %s %s%s %s > %s",
1157     rep, vw, res, opts,
1158     vval(OCTREE), rawfile);
1159     }
1160 greg 2.4 if (runcom(combuf)) { /* run rpict */
1161 greg 2.1 fprintf(stderr, "%s: error rendering view %s\n",
1162     progname, vs);
1163     exit(1);
1164     }
1165     /* build pfilt command */
1166     if (mult > 1)
1167     sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
1168     pfopts, mult, mult, rawfile, picfile);
1169     else
1170     sprintf(combuf, "pfilt%s %s > %s", pfopts,
1171     rawfile, picfile);
1172 greg 2.4 if (runcom(combuf)) { /* run pfilt */
1173 greg 2.1 fprintf(stderr,
1174     "%s: error filtering view %s\n\t%s removed\n",
1175     progname, vs, picfile);
1176     unlink(picfile);
1177     exit(1);
1178     }
1179     /* remove raw file */
1180 greg 2.4 rmfile(rawfile);
1181     }
1182     }
1183    
1184    
1185     runcom(cs) /* run command */
1186     char *cs;
1187     {
1188     if (!silent) /* echo it */
1189     printf("\t%s\n", cs);
1190     if (noaction)
1191     return(0);
1192     fflush(stdout); /* flush output and pass to shell */
1193     return(system(cs));
1194     }
1195    
1196    
1197     rmfile(fn) /* remove a file */
1198     char *fn;
1199     {
1200     if (!silent)
1201 greg 2.1 #ifdef MSDOS
1202 greg 2.4 printf("\tdel %s\n", fn);
1203 greg 2.1 #else
1204 greg 2.4 printf("\trm -f %s\n", fn);
1205 greg 2.1 #endif
1206 greg 2.4 if (noaction)
1207     return(0);
1208     return(unlink(fn));
1209 greg 2.3 }
1210 greg 2.8
1211    
1212     #ifdef MSDOS
1213     setenv(vname, value) /* set an environment variable */
1214     char *vname, *value;
1215     {
1216     register char *evp;
1217    
1218     evp = bmalloc(strlen(vname)+strlen(value)+2);
1219     if (evp == NULL)
1220     syserr(progname);
1221     sprintf(evp, "%s=%s", vname, value);
1222     if (putenv(evp) != 0) {
1223     fprintf(stderr, "%s: out of environment space\n", progname);
1224     exit(1);
1225     }
1226     if (!silent)
1227     printf("set %s\n", evp);
1228     }
1229     #endif
1230 greg 2.3
1231    
1232     badvalue(vc) /* report bad variable value and exit */
1233     int vc;
1234     {
1235     fprintf(stderr, "%s: bad value for variable '%s'\n",
1236     progname, vnam(vc));
1237     exit(1);
1238     }
1239    
1240    
1241     syserr(s) /* report a system error and exit */
1242     char *s;
1243     {
1244     perror(s);
1245     exit(1);
1246 greg 2.1 }