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

Comparing ray/src/util/rad.c (file contents):
Revision 2.37 by greg, Sat Jan 1 09:16:50 1994 UTC vs.
Revision 2.50 by greg, Fri Jan 5 14:43:03 1996 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1994 Regents of the University of California */
1 > /* Copyright (c) 1995 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 10 | Line 10 | static char SCCSid[] = "$SunId$ LBL";
10  
11   #include "standard.h"
12   #include "paths.h"
13 + #include "vars.h"
14   #include <ctype.h>
15   #include <sys/types.h>
16  
16 –
17 – typedef struct {
18 –        char    *name;          /* variable name */
19 –        short   nick;           /* # characters required for nickname */
20 –        short   nass;           /* # assignments made */
21 –        char    *value;         /* assigned value(s) */
22 –        int     (*fixval)();    /* assignment checking function */
23 – } VARIABLE;
24 –
25 – int     onevalue(), catvalues(), boolvalue(),
26 –        qualvalue(), fltvalue(), intvalue();
27 –
17                                  /* variables */
18   #define OBJECT          0               /* object files */
19   #define SCENE           1               /* scene files */
# Line 38 | Line 27 | int    onevalue(), catvalues(), boolvalue(),
27   #define ZONE            9               /* simulation zone */
28   #define QUALITY         10              /* desired rendering quality */
29   #define OCTREE          11              /* octree file name */
30 < #define PICTURE         12              /* picture file name */
30 > #define PICTURE         12              /* picture file root name */
31   #define AMBFILE         13              /* ambient file name */
32   #define OPTFILE         14              /* rendering options file */
33   #define EXPOSURE        15              /* picture exposure setting */
# Line 49 | Line 38 | int    onevalue(), catvalues(), boolvalue(),
38   #define PENUMBRAS       20              /* shadow penumbras are desired */
39   #define VARIABILITY     21              /* level of light variability */
40   #define REPORT          22              /* report frequency and errfile */
41 + #define RAWFILE         23              /* raw picture file root name */
42 + #define ZFILE           24              /* distance file root name */
43                                  /* total number of variables */
44 < #define NVARS           23
44 > int NVARS = 25;
45  
46 < VARIABLE        vv[NVARS] = {           /* variable-value pairs */
46 > VARIABLE        vv[] = {                /* variable-value pairs */
47          {"objects",     3,      0,      NULL,   catvalues},
48          {"scene",       3,      0,      NULL,   catvalues},
49          {"materials",   3,      0,      NULL,   catvalues},
# Line 76 | Line 67 | VARIABLE       vv[NVARS] = {           /* variable-value pairs */
67          {"PENUMBRAS",   3,      0,      NULL,   boolvalue},
68          {"VARIABILITY", 3,      0,      NULL,   qualvalue},
69          {"REPORT",      3,      0,      NULL,   onevalue},
70 +        {"RAWFILE",     3,      0,      NULL,   onevalue},
71 +        {"ZFILE",       2,      0,      NULL,   onevalue},
72   };
73  
81 – VARIABLE        *matchvar();
82 – char    *nvalue();
83 –
84 – #define UPPER(c)        ((c)&~0x20)     /* ASCII trick */
85 –
86 – #define vnam(vc)        (vv[vc].name)
87 – #define vdef(vc)        (vv[vc].nass)
88 – #define vval(vc)        (vv[vc].value)
89 – #define vint(vc)        atoi(vval(vc))
90 – #define vlet(vc)        UPPER(vval(vc)[0])
91 – #define vscale          vlet
92 – #define vbool(vc)       (vlet(vc)=='T')
93 –
94 – #define HIGH            'H'
95 – #define MEDIUM          'M'
96 – #define LOW             'L'
97 –
74                                  /* overture calculation file */
75   #ifdef NIX
76 < char    overfile[] = "overture.raw";
76 > char    overfile[] = "overture.unf";
77   #else
78   char    overfile[] = "/dev/null";
79   #endif
# Line 117 | Line 93 | time_t oct1date;               /* date of post-mkillum octree (>= m
93   int     nowarn = 0;             /* no warnings */
94   int     explicate = 0;          /* explicate variables */
95   int     silent = 0;             /* do work silently */
96 + int     touchonly = 0;          /* touch files only */
97   int     noaction = 0;           /* don't do anything */
98   int     sayview = 0;            /* print view out */
99   char    *rvdevice = NULL;       /* rview output device */
# Line 135 | Line 112 | int    argc;
112   char    *argv[];
113   {
114          char    ropts[512];
115 +        char    popts[64];
116          int     i;
117  
118          progname = argv[0];
# Line 147 | Line 125 | char   *argv[];
125                  case 'n':
126                          noaction++;
127                          break;
128 +                case 't':
129 +                        touchonly++;
130 +                        break;
131                  case 'e':
132                          explicate++;
133                          break;
# Line 171 | Line 152 | char   *argv[];
152                                  /* assign Radiance root file name */
153          rootname(radname, rifname);
154                                  /* load variable values */
155 <        load(rifname);
155 >        loadvars(rifname);
156                                  /* get any additional assignments */
157          for (i++; i < argc; i++)
158                  setvariable(argv[i]);
# Line 183 | Line 164 | char   *argv[];
164          setdefaults();
165                                  /* print all values if requested */
166          if (explicate)
167 <                printvals();
167 >                printvars(stdout);
168                                  /* build octree (and run mkillum) */
169          oconv();
170                                  /* check date on ambient file */
171          checkambfile();
172                                  /* run simulation */
173 <        renderopts(ropts);
173 >        renderopts(ropts, popts);
174          xferopts(ropts);
175          if (rvdevice != NULL)
176 <                rview(ropts);
176 >                rview(ropts, popts);
177          else
178 <                rpict(ropts);
178 >                rpict(ropts, popts);
179          exit(0);
180   userr:
181          fprintf(stderr,
182 <        "Usage: %s [-s][-n][-e][-V][-v view][-o dev] rfile [VAR=value ..]\n",
182 > "Usage: %s [-s][-n][-t][-e][-V][-v view][-o dev] rfile [VAR=value ..]\n",
183                          progname);
184          exit(1);
185   }
# Line 219 | Line 200 | register char  *rn, *fn;
200   }
201  
202  
222 – #define NOCHAR  127             /* constant for character to delete */
223 –
224 –
225 – load(rfname)                    /* load Radiance simulation file */
226 – char    *rfname;
227 – {
228 –        FILE    *fp;
229 –        char    buf[512];
230 –        register char   *cp;
231 –
232 –        if (rfname == NULL)
233 –                fp = stdin;
234 –        else if ((fp = fopen(rfname, "r")) == NULL)
235 –                syserr(rfname);
236 –        while (fgetline(buf, sizeof(buf), fp) != NULL) {
237 –                for (cp = buf; *cp; cp++) {
238 –                        switch (*cp) {
239 –                        case '\\':
240 –                                *cp++ = NOCHAR;
241 –                                continue;
242 –                        case '#':
243 –                                *cp = '\0';
244 –                                break;
245 –                        default:
246 –                                continue;
247 –                        }
248 –                        break;
249 –                }
250 –                setvariable(buf);
251 –        }
252 –        fclose(fp);
253 – }
254 –
255 –
256 – setvariable(ass)                /* assign variable according to string */
257 – register char   *ass;
258 – {
259 –        char    varname[32];
260 –        int     n;
261 –        register char   *cp;
262 –        register VARIABLE       *vp;
263 –        register int    i;
264 –
265 –        while (isspace(*ass))           /* skip leading space */
266 –                ass++;
267 –        cp = varname;                   /* extract name */
268 –        while (cp < varname+sizeof(varname)-1
269 –                        && *ass && !isspace(*ass) && *ass != '=')
270 –                *cp++ = *ass++;
271 –        *cp = '\0';
272 –        if (!varname[0])
273 –                return;         /* no variable name! */
274 –                                        /* trim value */
275 –        while (isspace(*ass) || *ass == '=')
276 –                ass++;
277 –        for (n = strlen(ass); n > 0; n--)
278 –                if (!isspace(ass[n-1]))
279 –                        break;
280 –        if (!n && !nowarn) {
281 –                fprintf(stderr, "%s: warning - missing value for variable '%s'\n",
282 –                                progname, varname);
283 –                return;
284 –        }
285 –                                        /* match variable from list */
286 –        vp = matchvar(varname);
287 –        if (vp == NULL) {
288 –                fprintf(stderr, "%s: unknown variable '%s'\n",
289 –                                progname, varname);
290 –                exit(1);
291 –        }
292 –                                        /* assign new value */
293 –        if (i = vp->nass) {
294 –                cp = vp->value;
295 –                while (i--)
296 –                        while (*cp++)
297 –                                ;
298 –                i = cp - vp->value;
299 –                vp->value = realloc(vp->value, i+n+1);
300 –        } else
301 –                vp->value = malloc(n+1);
302 –        if (vp->value == NULL)
303 –                syserr(progname);
304 –        cp = vp->value+i;               /* copy value, squeezing spaces */
305 –        *cp = *ass;
306 –        for (i = 1; i <= n; i++) {
307 –                if (ass[i] == NOCHAR)
308 –                        continue;
309 –                if (isspace(*cp))
310 –                        while (isspace(ass[i]))
311 –                                i++;
312 –                *++cp = ass[i];
313 –        }
314 –        if (isspace(*cp))               /* remove trailing space */
315 –                *cp = '\0';
316 –        vp->nass++;
317 – }
318 –
319 –
320 – VARIABLE *
321 – matchvar(nam)                   /* match a variable by its name */
322 – char    *nam;
323 – {
324 –        int     n = strlen(nam);
325 –        register int    i;
326 –
327 –        for (i = 0; i < NVARS; i++)
328 –                if (n >= vv[i].nick && !strncmp(nam, vv[i].name, n))
329 –                        return(vv+i);
330 –        return(NULL);
331 – }
332 –
333 –
334 – char *
335 – nvalue(vp, n)                   /* return nth variable value */
336 – VARIABLE        *vp;
337 – register int    n;
338 – {
339 –        register char   *cp;
340 –
341 –        if (vp == NULL | n < 0 | n >= vp->nass)
342 –                return(NULL);
343 –        cp = vp->value;
344 –        while (n--)
345 –                while (*cp++)
346 –                        ;
347 –        return(cp);
348 – }
349 –
350 –
351 – checkvalues()                   /* check assignments */
352 – {
353 –        register int    i;
354 –
355 –        for (i = 0; i < NVARS; i++)
356 –                if (vv[i].fixval != NULL)
357 –                        (*vv[i].fixval)(vv+i);
358 – }
359 –
360 –
361 – onevalue(vp)                    /* only one assignment for this variable */
362 – register VARIABLE       *vp;
363 – {
364 –        if (vp->nass < 2)
365 –                return;
366 –        if (!nowarn)
367 –                fprintf(stderr,
368 –                "%s: warning - multiple assignment of variable '%s'\n",
369 –                        progname, vp->name);
370 –        do
371 –                vp->value += strlen(vp->value)+1;
372 –        while (--vp->nass > 1);
373 – }
374 –
375 –
376 – catvalues(vp)                   /* concatenate variable values */
377 – register VARIABLE       *vp;
378 – {
379 –        register char   *cp;
380 –
381 –        if (vp->nass < 2)
382 –                return;
383 –        for (cp = vp->value; vp->nass > 1; vp->nass--) {
384 –                while (*cp)
385 –                        cp++;
386 –                *cp++ = ' ';
387 –        }
388 – }
389 –
390 –
391 – int
392 – badmatch(tv, cv)                /* case insensitive truncated comparison */
393 – register char   *tv, *cv;
394 – {
395 –        if (!*tv) return(1);            /* null string cannot match */
396 –        do
397 –                if (UPPER(*tv) != *cv++)
398 –                        return(1);
399 –        while (*++tv);
400 –        return(0);              /* OK */
401 – }
402 –
403 –
404 – boolvalue(vp)                   /* check boolean for legal values */
405 – register VARIABLE       *vp;
406 – {
407 –        if (!vp->nass) return;
408 –        onevalue(vp);
409 –        switch (UPPER(vp->value[0])) {
410 –        case 'T':
411 –                if (badmatch(vp->value, "TRUE")) break;
412 –                return;
413 –        case 'F':
414 –                if (badmatch(vp->value, "FALSE")) break;
415 –                return;
416 –        }
417 –        fprintf(stderr, "%s: illegal value for boolean variable '%s'\n",
418 –                        progname, vp->name);
419 –        exit(1);
420 – }
421 –
422 –
423 – qualvalue(vp)                   /* check qualitative var. for legal values */
424 – register VARIABLE       *vp;
425 – {
426 –        if (!vp->nass) return;
427 –        onevalue(vp);
428 –        switch (UPPER(vp->value[0])) {
429 –        case 'L':
430 –                if (badmatch(vp->value, "LOW")) break;
431 –                return;
432 –        case 'M':
433 –                if (badmatch(vp->value, "MEDIUM")) break;
434 –                return;
435 –        case 'H':
436 –                if (badmatch(vp->value, "HIGH")) break;
437 –                return;
438 –        }
439 –        fprintf(stderr, "%s: illegal value for qualitative variable '%s'\n",
440 –                        progname, vp->name);
441 –        exit(1);
442 – }
443 –
444 –
445 – intvalue(vp)                    /* check integer variable for legal values */
446 – register VARIABLE       *vp;
447 – {
448 –        if (!vp->nass) return;
449 –        onevalue(vp);
450 –        if (isint(vp->value)) return;
451 –        fprintf(stderr, "%s: illegal value for integer variable '%s'\n",
452 –                        progname, vp->name);
453 –        exit(1);
454 – }
455 –
456 –
457 – fltvalue(vp)                    /* check float variable for legal values */
458 – register VARIABLE       *vp;
459 – {
460 –        if (!vp->nass) return;
461 –        onevalue(vp);
462 –        if (isflt(vp->value)) return;
463 –        fprintf(stderr, "%s: illegal value for real variable '%s'\n",
464 –                        progname, vp->name);
465 –        exit(1);
466 – }
467 –
468 –
203   time_t
204   checklast(fnames)                       /* check files and find most recent */
205   register char   *fnames;
# Line 554 | Line 288 | double org[3], *sizp;
288          extern FILE     *popen();
289          static double   oorg[3], osiz = 0.;
290          double  min[3], max[3];
291 <        char    buf[512];
291 >        char    buf[1024];
292          FILE    *fp;
293          register int    i;
294  
# Line 645 | Line 379 | setdefaults()                  /* set default values for unassigned v
379   }
380  
381  
648 – printvals()                     /* print variable values */
649 – {
650 –        int     i, j, clipline;
651 –        register char   *cp;
652 –        register int    k;
653 –
654 –        for (i = 0; i < NVARS; i++)             /* print each variable */
655 –            for (j = 0; j < vdef(i); j++) {     /* print each assignment */
656 –                fputs(vnam(i), stdout);
657 –                fputs("= ", stdout);
658 –                k = clipline = ( vv[i].fixval == catvalues ? 64 : 320 )
659 –                                - strlen(vnam(i)) ;
660 –                cp = nvalue(vv+i, j);
661 –                while (*cp) {
662 –                    putchar(*cp++);
663 –                    if (--k <= 0) {             /* line too long */
664 –                        while (*cp && !isspace(*cp))
665 –                            putchar(*cp++);     /* finish this word */
666 –                        if (*cp) {              /* start new line */
667 –                            putchar('\n');
668 –                            fputs(vnam(i), stdout);
669 –                            putchar('=');
670 –                            k = clipline;
671 –                        }
672 –                    }
673 –                }
674 –                putchar('\n');
675 –            }
676 –        fflush(stdout);
677 – }
678 –
679 –
382   oconv()                         /* run oconv and mkillum if necessary */
383   {
384          static char     illumtmp[] = "ilXXXXXX";
385 <        char    combuf[512], ocopts[64], mkopts[64];
385 >        char    combuf[1024], ocopts[64], mkopts[64];
386  
387          oconvopts(ocopts);              /* get options */
388          if (octreedate < scenedate) {   /* check date on original octree */
389 <                                                /* build command */
390 <                if (vdef(MATERIAL))
391 <                        sprintf(combuf, "oconv%s %s %s > %s", ocopts,
392 <                                vval(MATERIAL), vval(SCENE), vval(OCTREE));
393 <                else
394 <                        sprintf(combuf, "oconv%s %s > %s", ocopts,
395 <                                        vval(SCENE), vval(OCTREE));
396 <                
397 <                if (runcom(combuf)) {           /* run it */
398 <                        fprintf(stderr,
389 >                if (touchonly && octreedate)
390 >                        touch(vval(OCTREE));
391 >                else {                          /* build command */
392 >                        if (vdef(MATERIAL))
393 >                                sprintf(combuf, "oconv%s %s %s > %s", ocopts,
394 >                                                vval(MATERIAL), vval(SCENE),
395 >                                                vval(OCTREE));
396 >                        else
397 >                                sprintf(combuf, "oconv%s %s > %s", ocopts,
398 >                                                vval(SCENE), vval(OCTREE));
399 >                        
400 >                        if (runcom(combuf)) {           /* run it */
401 >                                fprintf(stderr,
402                                  "%s: error generating octree\n\t%s removed\n",
403 <                                        progname, vval(OCTREE));
404 <                        unlink(vval(OCTREE));
405 <                        exit(1);
403 >                                                progname, vval(OCTREE));
404 >                                unlink(vval(OCTREE));
405 >                                exit(1);
406 >                        }
407                  }
408                  octreedate = time((time_t *)NULL);
409 +                if (octreedate < scenedate)     /* in case clock is off */
410 +                        octreedate = scenedate;
411          }
412          if (oct1name == vval(OCTREE))           /* no mkillum? */
413                  oct1date = octreedate > matdate ? octreedate : matdate;
# Line 708 | Line 416 | oconv()                                /* run oconv and mkillum if necessary */
416                  return;
417                                                  /* make octree0 */
418          if (oct0date < scenedate | oct0date < illumdate) {
419 <                                                /* build command */
419 >                if (touchonly && oct0date)
420 >                        touch(oct0name);
421 >                else {                          /* build command */
422 >                        if (octreedate)
423 >                                sprintf(combuf, "oconv%s -i %s %s > %s", ocopts,
424 >                                        vval(OCTREE), vval(ILLUM), oct0name);
425 >                        else if (vdef(MATERIAL))
426 >                                sprintf(combuf, "oconv%s %s %s > %s", ocopts,
427 >                                        vval(MATERIAL), vval(ILLUM), oct0name);
428 >                        else
429 >                                sprintf(combuf, "oconv%s %s > %s", ocopts,
430 >                                        vval(ILLUM), oct0name);
431 >                        if (runcom(combuf)) {           /* run it */
432 >                                fprintf(stderr,
433 >                                "%s: error generating octree\n\t%s removed\n",
434 >                                                progname, oct0name);
435 >                                unlink(oct0name);
436 >                                exit(1);
437 >                        }
438 >                }
439 >                oct0date = time((time_t *)NULL);
440 >                if (oct0date < octreedate)      /* in case clock is off */
441 >                        oct0date = octreedate;
442 >                if (oct0date < illumdate)       /* ditto */
443 >                        oct0date = illumdate;
444 >                }
445 >        if (touchonly && oct1date)
446 >                touch(oct1name);
447 >        else {
448 >                mkillumopts(mkopts);            /* build mkillum command */
449 >                mktemp(illumtmp);
450 >                sprintf(combuf, "mkillum%s %s \"<\" %s > %s", mkopts,
451 >                                oct0name, vval(ILLUM), illumtmp);
452 >                if (runcom(combuf)) {                   /* run it */
453 >                        fprintf(stderr, "%s: error running mkillum\n",
454 >                                        progname);
455 >                        unlink(illumtmp);
456 >                        exit(1);
457 >                }
458 >                                                /* make octree1 (frozen) */
459                  if (octreedate)
460 <                        sprintf(combuf, "oconv%s -i %s %s > %s", ocopts,
461 <                                vval(OCTREE), vval(ILLUM), oct0name);
460 >                        sprintf(combuf, "oconv%s -f -i %s %s > %s", ocopts,
461 >                                vval(OCTREE), illumtmp, oct1name);
462                  else if (vdef(MATERIAL))
463 <                        sprintf(combuf, "oconv%s %s %s > %s", ocopts,
464 <                                vval(MATERIAL), vval(ILLUM), oct0name);
463 >                        sprintf(combuf, "oconv%s -f %s %s > %s", ocopts,
464 >                                vval(MATERIAL), illumtmp, oct1name);
465                  else
466 <                        sprintf(combuf, "oconv%s %s > %s", ocopts,
467 <                                vval(ILLUM), oct0name);
466 >                        sprintf(combuf, "oconv%s -f %s > %s", ocopts,
467 >                                illumtmp, oct1name);
468                  if (runcom(combuf)) {           /* run it */
469                          fprintf(stderr,
470                                  "%s: error generating octree\n\t%s removed\n",
471 <                                        progname, oct0name);
472 <                        unlink(oct0name);
471 >                                        progname, oct1name);
472 >                        unlink(oct1name);
473 >                        unlink(illumtmp);
474                          exit(1);
475                  }
476 <                oct0date = time((time_t *)NULL);
476 >                rmfile(illumtmp);
477          }
730 –        mkillumopts(mkopts);                    /* build mkillum command */
731 –        mktemp(illumtmp);
732 –        sprintf(combuf, "mkillum%s %s \"<\" %s > %s", mkopts,
733 –                        oct0name, vval(ILLUM), illumtmp);
734 –        if (runcom(combuf)) {                   /* run it */
735 –                fprintf(stderr, "%s: error running mkillum\n", progname);
736 –                unlink(illumtmp);
737 –                exit(1);
738 –        }
739 –                                                /* make octree1 (frozen) */
740 –        if (octreedate)
741 –                sprintf(combuf, "oconv%s -f -i %s %s > %s", ocopts,
742 –                        vval(OCTREE), illumtmp, oct1name);
743 –        else if (vdef(MATERIAL))
744 –                sprintf(combuf, "oconv%s -f %s %s > %s", ocopts,
745 –                        vval(MATERIAL), illumtmp, oct1name);
746 –        else
747 –                sprintf(combuf, "oconv%s -f %s > %s", ocopts,
748 –                        illumtmp, oct1name);
749 –        if (runcom(combuf)) {           /* run it */
750 –                fprintf(stderr,
751 –                        "%s: error generating octree\n\t%s removed\n",
752 –                                progname, oct1name);
753 –                unlink(oct1name);
754 –                exit(1);
755 –        }
478          oct1date = time((time_t *)NULL);
479 <        rmfile(illumtmp);
479 >        if (oct1date < oct0date)        /* in case clock is off */
480 >                oct1date = oct0date;
481   }
482  
483  
# Line 800 | Line 523 | checkambfile()                 /* check date on ambient file */
523          if (!(afdate = fdate(vval(AMBFILE))))
524                  return;
525          if (oct1date > afdate)
526 <                rmfile(vval(AMBFILE));
526 >                if (touchonly)
527 >                        touch(vval(AMBFILE));
528 >                else
529 >                        rmfile(vval(AMBFILE));
530   }
531  
532  
# Line 809 | Line 535 | ambval()                               /* compute ambient value */
535   {
536          if (vdef(EXPOSURE)) {
537                  if (vval(EXPOSURE)[0] == '+' || vval(EXPOSURE)[0] == '-')
538 <                        return(.5/pow(2.,atof(vval(EXPOSURE))));
539 <                return(.5/atof(vval(EXPOSURE)));
538 >                        return(.5/pow(2.,vflt(EXPOSURE)));
539 >                return(.5/vflt(EXPOSURE));
540          }
541          if (vlet(ZONE) == 'E')
542                  return(10.);
# Line 820 | Line 546 | ambval()                               /* compute ambient value */
546   }
547  
548  
549 < renderopts(op)                          /* set rendering options */
550 < char    *op;
549 > renderopts(op, po)                      /* set rendering options */
550 > char    *op, *po;
551   {
552          switch(vscale(QUALITY)) {
553          case LOW:
554 <                lowqopts(op);
554 >                lowqopts(op, po);
555                  break;
556          case MEDIUM:
557 <                medqopts(op);
557 >                medqopts(op, po);
558                  break;
559          case HIGH:
560 <                hiqopts(op);
560 >                hiqopts(op, po);
561                  break;
562          }
563   }
564  
565  
566 < lowqopts(op)                            /* low quality rendering options */
566 > lowqopts(op, po)                        /* low quality rendering options */
567   register char   *op;
568 + char    *po;
569   {
570          double  d, org[3], siz[3];
571  
572          *op = '\0';
573 +        *po = '\0';
574          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
575                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
576                  badvalue(ZONE);
577          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
578 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
579 +                badvalue(ZONE);
580          getoctcube(org, &d);
581          d *= 3./(siz[0]+siz[1]+siz[2]);
582          switch (vscale(DETAIL)) {
583          case LOW:
584 <                op = addarg(op, "-ps 16 -dp 64");
584 >                po = addarg(po, "-ps 16");
585 >                op = addarg(op, "-dp 64");
586                  sprintf(op, " -ar %d", (int)(4*d));
587                  op += strlen(op);
588                  break;
589          case MEDIUM:
590 <                op = addarg(op, "-ps 8 -dp 128");
590 >                po = addarg(po, "-ps 8");
591 >                op = addarg(op, "-dp 128");
592                  sprintf(op, " -ar %d", (int)(8*d));
593                  op += strlen(op);
594                  break;
595          case HIGH:
596 <                op = addarg(op, "-ps 4 -dp 256");
596 >                po = addarg(po, "-ps 4");
597 >                op = addarg(op, "-dp 256");
598                  sprintf(op, " -ar %d", (int)(16*d));
599                  op += strlen(op);
600                  break;
601          }
602 <        op = addarg(op, "-pt .16");
602 >        po = addarg(po, "-pt .16");
603          if (vbool(PENUMBRAS))
604                  op = addarg(op, "-ds .4");
605          else
# Line 898 | Line 631 | register char  *op;
631   }
632  
633  
634 < medqopts(op)                            /* medium quality rendering options */
634 > medqopts(op, po)                        /* medium quality rendering options */
635   register char   *op;
636 + char    *po;
637   {
638 <        double  d, org[3], siz[3];
638 >        double  d, org[3], siz[3], asz;
639  
640          *op = '\0';
641 +        *po = '\0';
642          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
643                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
644                  badvalue(ZONE);
645          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
646 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
647 +                badvalue(ZONE);
648          getoctcube(org, &d);
649 <        d *= 3./(siz[0]+siz[1]+siz[2]);
649 >        asz = (siz[0]+siz[1]+siz[2])/3.;
650 >        d /= asz;
651          switch (vscale(DETAIL)) {
652          case LOW:
653 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
653 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
654                  op = addarg(op, "-dp 256");
655                  sprintf(op, " -ar %d", (int)(8*d));
656                  op += strlen(op);
657 +                sprintf(op, " -ms %.2g", asz/20.);
658 +                op += strlen(op);
659                  break;
660          case MEDIUM:
661 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
661 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
662                  op = addarg(op, "-dp 512");
663                  sprintf(op, " -ar %d", (int)(16*d));
664                  op += strlen(op);
665 +                sprintf(op, " -ms %.2g", asz/40.);
666 +                op += strlen(op);
667                  break;
668          case HIGH:
669 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
669 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
670                  op = addarg(op, "-dp 1024");
671                  sprintf(op, " -ar %d", (int)(32*d));
672                  op += strlen(op);
673 +                sprintf(op, " -ms %.2g", asz/80.);
674 +                op += strlen(op);
675                  break;
676          }
677 <        op = addarg(op, "-pt .08");
677 >        po = addarg(po, "-pt .08");
678          if (vbool(PENUMBRAS))
679                  op = addarg(op, "-ds .2 -dj .5");
680          else
# Line 965 | Line 709 | register char  *op;
709   }
710  
711  
712 < hiqopts(op)                             /* high quality rendering options */
712 > hiqopts(op, po)                         /* high quality rendering options */
713   register char   *op;
714 + char    *po;
715   {
716 <        double  d, org[3], siz[3];
716 >        double  d, org[3], siz[3], asz;
717  
718          *op = '\0';
719 +        *po = '\0';
720          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
721                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
722                  badvalue(ZONE);
723          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
724 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
725 +                badvalue(ZONE);
726          getoctcube(org, &d);
727 <        d *= 3./(siz[0]+siz[1]+siz[2]);
727 >        asz = (siz[0]+siz[1]+siz[2])/3.;
728 >        d /= asz;
729          switch (vscale(DETAIL)) {
730          case LOW:
731 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
731 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
732                  op = addarg(op, "-dp 1024");
733                  sprintf(op, " -ar %d", (int)(16*d));
734                  op += strlen(op);
735 +                sprintf(op, " -ms %.2g", asz/40.);
736 +                op += strlen(op);
737                  break;
738          case MEDIUM:
739 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
739 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
740                  op = addarg(op, "-dp 2048");
741                  sprintf(op, " -ar %d", (int)(32*d));
742                  op += strlen(op);
743 +                sprintf(op, " -ms %.2g", asz/80.);
744 +                op += strlen(op);
745                  break;
746          case HIGH:
747 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
747 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
748                  op = addarg(op, "-dp 4096");
749                  sprintf(op, " -ar %d", (int)(64*d));
750                  op += strlen(op);
751 +                sprintf(op, " -ms %.2g", asz/160.);
752 +                op += strlen(op);
753                  break;
754          }
755 <        op = addarg(op, "-pt .04");
755 >        po = addarg(po, "-pt .04");
756          if (vbool(PENUMBRAS))
757 <                op = addarg(op, "-ds .1 -dj .7");
757 >                op = addarg(op, "-ds .1 -dj .65");
758          else
759                  op = addarg(op, "-ds .2");
760          op = addarg(op, "-dt .05 -dc .75 -dr 3 -sj 1 -st .01");
# Line 1041 | Line 796 | char   *ro;
796                  return;
797          if (vdef(OPTFILE)) {
798                  for (cp = ro; cp[1]; cp++)
799 <                        if (isspace(cp[1]) && cp[2] == '-' && isalpha(cp[3]))
799 >                        if (isspace(cp[1]) && (cp[2] == '@' ||
800 >                                        (cp[2] == '-' && isalpha(cp[3]))))
801                                  *cp = '\n';
802                          else
803                                  *cp = cp[1];
# Line 1136 | Line 892 | register char  *vs;
892                  zpos = -1; vs++;
893          }
894          viewtype = 'v';
895 <        if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h')
895 >        if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h' | *vs == 'c')
896                  viewtype = *vs++;
897          cp = viewopts;
898          if ((!*vs || isspace(*vs)) && (xpos|ypos|zpos)) {       /* got one! */
# Line 1147 | Line 903 | register char  *vs;
903                          badvalue(ZONE);
904                  for (i = 0; i < 3; i++) {
905                          dim[i] -= cent[i];
906 +                        if (dim[i] <= FTINY)
907 +                                badvalue(ZONE);
908                          cent[i] += .5*dim[i];
909                  }
910                  mult = vlet(ZONE)=='E' ? 2. : .45 ;
# Line 1189 | Line 947 | register char  *vs;
947                  case 'h':
948                          cp = addarg(cp, "-vh 180 -vv 180");
949                          break;
950 +                case 'c':
951 +                        cp = addarg(cp, "-vh 180 -vv 90");
952 +                        break;
953                  }
954          } else {
955                  while (!isspace(*vs))           /* else skip id */
# Line 1310 | Line 1071 | register char  *vopts;
1071   }
1072  
1073  
1074 < rview(opts)                             /* run rview with first view */
1075 < char    *opts;
1074 > rview(opts, po)                         /* run rview with first view */
1075 > char    *opts, *po;
1076   {
1077          char    *vw;
1078          char    combuf[512];
1079                                          /* build command */
1080 <        if ((vw = getview(0, NULL)) == NULL)
1080 >        if (touchonly || (vw = getview(0, NULL)) == NULL)
1081                  return;
1082          if (sayview)
1083                  printview(vw);
1084 <        sprintf(combuf, "rview %s%s -R %s ", vw, opts, rifname);
1084 >        sprintf(combuf, "rview %s%s%s -R %s ", vw, po, opts, rifname);
1085          if (rvdevice != NULL)
1086                  sprintf(combuf+strlen(combuf), "-o %s ", rvdevice);
1087          if (vdef(EXPOSURE))
# Line 1333 | Line 1094 | char   *opts;
1094   }
1095  
1096  
1097 < rpict(opts)                             /* run rpict and pfilt for each view */
1098 < char    *opts;
1097 > rpict(opts, po)                         /* run rpict and pfilt for each view */
1098 > char    *opts, *po;
1099   {
1100          char    combuf[1024];
1101 <        char    rawfile[MAXPATH], picfile[MAXPATH], rep[MAXPATH+16], res[32];
1101 >        char    rawfile[MAXPATH], picfile[MAXPATH];
1102 >        char    zopt[MAXPATH+4], rep[MAXPATH+16], res[32];
1103          char    pfopts[128];
1104          char    vs[32], *vw;
1105          int     vn, mult;
1106 +        time_t  rfdt, pfdt;
1107                                          /* get pfilt options */
1108          pfiltopts(pfopts);
1109                                          /* get resolution, reporting */
# Line 1389 | Line 1152 | char   *opts;
1152                  if (!vs[0])
1153                          sprintf(vs, "%d", vn);
1154                  sprintf(picfile, "%s_%s.pic", vval(PICTURE), vs);
1155 +                if (vdef(ZFILE))
1156 +                        sprintf(zopt, " -z %s_%s.zbf", vval(ZFILE), vs);
1157 +                else
1158 +                        zopt[0] = '\0';
1159                                                  /* check date on picture */
1160 <                if (fdate(picfile) >= oct1date)
1160 >                pfdt = fdate(picfile);
1161 >                if (pfdt >= oct1date)
1162                          continue;
1163 +                                                /* get raw file name */
1164 +                sprintf(rawfile, "%s_%s.unf",
1165 +                        vdef(RAWFILE) ? vval(RAWFILE) : vval(PICTURE), vs);
1166 +                rfdt = fdate(rawfile);
1167 +                if (touchonly) {                /* update times only */
1168 +                        if (rfdt) {
1169 +                                if (rfdt < oct1date)
1170 +                                        touch(rawfile);
1171 +                        } else if (pfdt && pfdt < oct1date)
1172 +                                touch(picfile);
1173 +                        continue;
1174 +                }
1175                                                  /* build rpict command */
1176 <                sprintf(rawfile, "%s_%s.raw", vval(PICTURE), vs);
1177 <                if (fdate(rawfile) >= oct1date)         /* recover */
1178 <                        sprintf(combuf, "rpict%s%s -ro %s %s",
1399 <                                        rep, opts, rawfile, oct1name);
1176 >                if (rfdt >= oct1date)           /* recover */
1177 >                        sprintf(combuf, "rpict%s%s%s%s -ro %s %s",
1178 >                                rep, po, opts, zopt, rawfile, oct1name);
1179                  else {
1180                          if (overture) {         /* run overture calculation */
1181                                  sprintf(combuf,
# Line 1413 | Line 1192 | char   *opts;
1192                                  rmfile(overfile);
1193   #endif
1194                          }
1195 <                        sprintf(combuf, "rpict%s %s %s%s %s > %s",
1196 <                                        rep, vw, res, opts,
1197 <                                oct1name, rawfile);
1195 >                        sprintf(combuf, "rpict%s %s %s%s%s%s %s > %s",
1196 >                                        rep, vw, res, po, opts, zopt,
1197 >                                        oct1name, rawfile);
1198                  }
1199                  if (runcom(combuf)) {           /* run rpict */
1200                          fprintf(stderr, "%s: error rendering view %s\n",
1201                                          progname, vs);
1202                          exit(1);
1203                  }
1204 +                if (!vdef(RAWFILE) || strcmp(vval(RAWFILE),vval(PICTURE))) {
1205                                                  /* build pfilt command */
1206 <                if (mult > 1)
1207 <                        sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
1206 >                        if (mult > 1)
1207 >                                sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
1208                                          pfopts, mult, mult, rawfile, picfile);
1209 <                else
1210 <                        sprintf(combuf, "pfilt%s %s > %s", pfopts,
1211 <                                        rawfile, picfile);
1212 <                if (runcom(combuf)) {           /* run pfilt */
1213 <                        fprintf(stderr,
1214 <                        "%s: error filtering view %s\n\t%s removed\n",
1215 <                                        progname, vs, picfile);
1216 <                        unlink(picfile);
1217 <                        exit(1);
1209 >                        else
1210 >                                sprintf(combuf, "pfilt%s %s > %s", pfopts,
1211 >                                                rawfile, picfile);
1212 >                        if (runcom(combuf)) {           /* run pfilt */
1213 >                                fprintf(stderr,
1214 >                                "%s: error filtering view %s\n\t%s removed\n",
1215 >                                                progname, vs, picfile);
1216 >                                unlink(picfile);
1217 >                                exit(1);
1218 >                        }
1219                  }
1220 <                                                /* remove raw file */
1221 <                rmfile(rawfile);
1220 >                                                /* remove/rename raw file */
1221 >                if (vdef(RAWFILE)) {
1222 >                        sprintf(combuf, "%s_%s.pic", vval(RAWFILE), vs);
1223 >                        mvfile(rawfile, combuf);
1224 >                } else
1225 >                        rmfile(rawfile);
1226          }
1227   }
1228  
1229  
1230 + touch(fn)                       /* update a file */
1231 + char    *fn;
1232 + {
1233 +        if (!silent)
1234 +                printf("\ttouch %s\n", fn);
1235 +        if (noaction)
1236 +                return(0);
1237 + #ifdef notused
1238 +        if (access(fn, F_OK) == -1)             /* create it */
1239 +                if (close(open(fn, O_WRONLY|O_CREAT, 0666)) == -1)
1240 +                        return(-1);
1241 + #endif
1242 +        return(setfdate(fn, time((time_t *)NULL)));
1243 + }
1244 +
1245 +
1246   runcom(cs)                      /* run command */
1247   char    *cs;
1248   {
# Line 1466 | Line 1267 | char   *fn;
1267          if (noaction)
1268                  return(0);
1269          return(unlink(fn));
1270 + }
1271 +
1272 +
1273 + mvfile(fold, fnew)              /* move a file */
1274 + char    *fold, *fnew;
1275 + {
1276 +        if (!silent)
1277 + #ifdef MSDOS
1278 +                printf("\trename %s %s\n", fold, fnew);
1279 + #else
1280 +                printf("\tmv %s %s\n", fold, fnew);
1281 + #endif
1282 +        if (noaction)
1283 +                return(0);
1284 +        return(rename(fold, fnew));
1285   }
1286  
1287  

Diff Legend

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