ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rad.c
Revision: 2.64
Committed: Thu Jun 26 00:58:11 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.63: +5 -2 lines
Log Message:
Abstracted process and path handling for Windows.
Renamed FLOAT to RREAL because of conflict on Windows.
Added conditional compiles for some signal handlers.

File Contents

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