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.36 by greg, Thu Dec 30 08:20:45 1993 UTC vs.
Revision 2.50 by greg, Fri Jan 5 14:43:03 1996 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1993 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 114 | Line 90 | time_t oct0date;               /* date of pre-mkillum octree */
90   char    *oct1name;              /* name of post-mkillum octree */
91   time_t  oct1date;               /* date of post-mkillum octree (>= matdate) */
92  
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 134 | 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 146 | 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 158 | Line 140 | char   *argv[];
140                  case 'v':
141                          viewselect = argv[++i];
142                          break;
143 +                case 'w':
144 +                        nowarn++;
145 +                        break;
146                  default:
147                          goto userr;
148                  }
# Line 167 | 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 179 | 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 215 | Line 200 | register char  *rn, *fn;
200   }
201  
202  
218 #define NOCHAR  127             /* constant for character to delete */
219
220
221 load(rfname)                    /* load Radiance simulation file */
222 char    *rfname;
223 {
224        FILE    *fp;
225        char    buf[512];
226        register char   *cp;
227
228        if (rfname == NULL)
229                fp = stdin;
230        else if ((fp = fopen(rfname, "r")) == NULL)
231                syserr(rfname);
232        while (fgetline(buf, sizeof(buf), fp) != NULL) {
233                for (cp = buf; *cp; cp++) {
234                        switch (*cp) {
235                        case '\\':
236                                *cp++ = NOCHAR;
237                                continue;
238                        case '#':
239                                *cp = '\0';
240                                break;
241                        default:
242                                continue;
243                        }
244                        break;
245                }
246                setvariable(buf);
247        }
248        fclose(fp);
249 }
250
251
252 setvariable(ass)                /* assign variable according to string */
253 register char   *ass;
254 {
255        char    varname[32];
256        int     n;
257        register char   *cp;
258        register VARIABLE       *vp;
259        register int    i;
260
261        while (isspace(*ass))           /* skip leading space */
262                ass++;
263        cp = varname;                   /* extract name */
264        while (cp < varname+sizeof(varname)-1
265                        && *ass && !isspace(*ass) && *ass != '=')
266                *cp++ = *ass++;
267        *cp = '\0';
268        if (!varname[0])
269                return;         /* no variable name! */
270                                        /* trim value */
271        while (isspace(*ass) || *ass == '=')
272                ass++;
273        for (n = strlen(ass); n > 0; n--)
274                if (!isspace(ass[n-1]))
275                        break;
276        if (!n) {
277                fprintf(stderr, "%s: warning - missing value for variable '%s'\n",
278                                progname, varname);
279                return;
280        }
281                                        /* match variable from list */
282        vp = matchvar(varname);
283        if (vp == NULL) {
284                fprintf(stderr, "%s: unknown variable '%s'\n",
285                                progname, varname);
286                exit(1);
287        }
288                                        /* assign new value */
289        if (i = vp->nass) {
290                cp = vp->value;
291                while (i--)
292                        while (*cp++)
293                                ;
294                i = cp - vp->value;
295                vp->value = realloc(vp->value, i+n+1);
296        } else
297                vp->value = malloc(n+1);
298        if (vp->value == NULL)
299                syserr(progname);
300        cp = vp->value+i;               /* copy value, squeezing spaces */
301        *cp = *ass;
302        for (i = 1; i <= n; i++) {
303                if (ass[i] == NOCHAR)
304                        continue;
305                if (isspace(*cp))
306                        while (isspace(ass[i]))
307                                i++;
308                *++cp = ass[i];
309        }
310        if (isspace(*cp))               /* remove trailing space */
311                *cp = '\0';
312        vp->nass++;
313 }
314
315
316 VARIABLE *
317 matchvar(nam)                   /* match a variable by its name */
318 char    *nam;
319 {
320        int     n = strlen(nam);
321        register int    i;
322
323        for (i = 0; i < NVARS; i++)
324                if (n >= vv[i].nick && !strncmp(nam, vv[i].name, n))
325                        return(vv+i);
326        return(NULL);
327 }
328
329
330 char *
331 nvalue(vp, n)                   /* return nth variable value */
332 VARIABLE        *vp;
333 register int    n;
334 {
335        register char   *cp;
336
337        if (vp == NULL | n < 0 | n >= vp->nass)
338                return(NULL);
339        cp = vp->value;
340        while (n--)
341                while (*cp++)
342                        ;
343        return(cp);
344 }
345
346
347 checkvalues()                   /* check assignments */
348 {
349        register int    i;
350
351        for (i = 0; i < NVARS; i++)
352                if (vv[i].fixval != NULL)
353                        (*vv[i].fixval)(vv+i);
354 }
355
356
357 onevalue(vp)                    /* only one assignment for this variable */
358 register VARIABLE       *vp;
359 {
360        if (vp->nass < 2)
361                return;
362        fprintf(stderr, "%s: warning - multiple assignment of variable '%s'\n",
363                        progname, vp->name);
364        do
365                vp->value += strlen(vp->value)+1;
366        while (--vp->nass > 1);
367 }
368
369
370 catvalues(vp)                   /* concatenate variable values */
371 register VARIABLE       *vp;
372 {
373        register char   *cp;
374
375        if (vp->nass < 2)
376                return;
377        for (cp = vp->value; vp->nass > 1; vp->nass--) {
378                while (*cp)
379                        cp++;
380                *cp++ = ' ';
381        }
382 }
383
384
385 int
386 badmatch(tv, cv)                /* case insensitive truncated comparison */
387 register char   *tv, *cv;
388 {
389        if (!*tv) return(1);            /* null string cannot match */
390        do
391                if (UPPER(*tv) != *cv++)
392                        return(1);
393        while (*++tv);
394        return(0);              /* OK */
395 }
396
397
398 boolvalue(vp)                   /* check boolean for legal values */
399 register VARIABLE       *vp;
400 {
401        if (!vp->nass) return;
402        onevalue(vp);
403        switch (UPPER(vp->value[0])) {
404        case 'T':
405                if (badmatch(vp->value, "TRUE")) break;
406                return;
407        case 'F':
408                if (badmatch(vp->value, "FALSE")) break;
409                return;
410        }
411        fprintf(stderr, "%s: illegal value for boolean variable '%s'\n",
412                        progname, vp->name);
413        exit(1);
414 }
415
416
417 qualvalue(vp)                   /* check qualitative var. for legal values */
418 register VARIABLE       *vp;
419 {
420        if (!vp->nass) return;
421        onevalue(vp);
422        switch (UPPER(vp->value[0])) {
423        case 'L':
424                if (badmatch(vp->value, "LOW")) break;
425                return;
426        case 'M':
427                if (badmatch(vp->value, "MEDIUM")) break;
428                return;
429        case 'H':
430                if (badmatch(vp->value, "HIGH")) break;
431                return;
432        }
433        fprintf(stderr, "%s: illegal value for qualitative variable '%s'\n",
434                        progname, vp->name);
435        exit(1);
436 }
437
438
439 intvalue(vp)                    /* check integer variable for legal values */
440 register VARIABLE       *vp;
441 {
442        if (!vp->nass) return;
443        onevalue(vp);
444        if (isint(vp->value)) return;
445        fprintf(stderr, "%s: illegal value for integer variable '%s'\n",
446                        progname, vp->name);
447        exit(1);
448 }
449
450
451 fltvalue(vp)                    /* check float variable for legal values */
452 register VARIABLE       *vp;
453 {
454        if (!vp->nass) return;
455        onevalue(vp);
456        if (isflt(vp->value)) return;
457        fprintf(stderr, "%s: illegal value for real variable '%s'\n",
458                        progname, vp->name);
459        exit(1);
460 }
461
462
203   time_t
204   checklast(fnames)                       /* check files and find most recent */
205   register char   *fnames;
# Line 548 | 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 639 | Line 379 | setdefaults()                  /* set default values for unassigned v
379   }
380  
381  
642 printvals()                     /* print variable values */
643 {
644        int     i, j, clipline;
645        register char   *cp;
646        register int    k;
647
648        for (i = 0; i < NVARS; i++)             /* print each variable */
649            for (j = 0; j < vdef(i); j++) {     /* print each assignment */
650                fputs(vnam(i), stdout);
651                fputs("= ", stdout);
652                k = clipline = ( vv[i].fixval == catvalues ? 64 : 320 )
653                                - strlen(vnam(i)) ;
654                cp = nvalue(vv+i, j);
655                while (*cp) {
656                    putchar(*cp++);
657                    if (--k <= 0) {             /* line too long */
658                        while (*cp && !isspace(*cp))
659                            putchar(*cp++);     /* finish this word */
660                        if (*cp) {              /* start new line */
661                            putchar('\n');
662                            fputs(vnam(i), stdout);
663                            putchar('=');
664                            k = clipline;
665                        }
666                    }
667                }
668                putchar('\n');
669            }
670        fflush(stdout);
671 }
672
673
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 702 | 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          }
724        mkillumopts(mkopts);                    /* build mkillum command */
725        mktemp(illumtmp);
726        sprintf(combuf, "mkillum%s %s \"<\" %s > %s", mkopts,
727                        oct0name, vval(ILLUM), illumtmp);
728        if (runcom(combuf)) {                   /* run it */
729                fprintf(stderr, "%s: error running mkillum\n", progname);
730                unlink(illumtmp);
731                exit(1);
732        }
733                                                /* make octree1 (frozen) */
734        if (octreedate)
735                sprintf(combuf, "oconv%s -f -i %s %s > %s", ocopts,
736                        vval(OCTREE), illumtmp, oct1name);
737        else if (vdef(MATERIAL))
738                sprintf(combuf, "oconv%s -f %s %s > %s", ocopts,
739                        vval(MATERIAL), illumtmp, oct1name);
740        else
741                sprintf(combuf, "oconv%s -f %s > %s", ocopts,
742                        illumtmp, oct1name);
743        if (runcom(combuf)) {           /* run it */
744                fprintf(stderr,
745                        "%s: error generating octree\n\t%s removed\n",
746                                progname, oct1name);
747                unlink(oct1name);
748                exit(1);
749        }
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 794 | 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 803 | 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 814 | 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 892 | 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 959 | 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 1035 | 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 1130 | 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 1141 | 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 1183 | 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 1304 | 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 1327 | 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 1383 | 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",
1393 <                                        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 1407 | 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 1460 | 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