ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rad.c
Revision: 2.69
Committed: Mon Jul 14 20:02:30 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.68: +12 -12 lines
Log Message: 
Moved some more platform dependencies to common header files.
Included a few necessary system headers.

File Contents

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