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.27 by greg, Sat Sep 25 13:21:35 1993 UTC vs.
Revision 2.47 by greg, Fri Dec 8 19:55:50 1995 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 11 | Line 11 | static char SCCSid[] = "$SunId$ LBL";
11   #include "standard.h"
12   #include "paths.h"
13   #include <ctype.h>
14 + #include <sys/types.h>
15  
16  
17   typedef struct {
# Line 21 | Line 22 | typedef struct {
22          int     (*fixval)();    /* assignment checking function */
23   } VARIABLE;
24  
25 < int     onevalue(), catvalues();
25 > int     onevalue(), catvalues(), boolvalue(),
26 >        qualvalue(), fltvalue(), intvalue();
27  
28                                  /* variables */
29   #define OBJECT          0               /* object files */
# Line 36 | Line 38 | int    onevalue(), catvalues();
38   #define ZONE            9               /* simulation zone */
39   #define QUALITY         10              /* desired rendering quality */
40   #define OCTREE          11              /* octree file name */
41 < #define PICTURE         12              /* picture file name */
41 > #define PICTURE         12              /* picture file root name */
42   #define AMBFILE         13              /* ambient file name */
43   #define OPTFILE         14              /* rendering options file */
44   #define EXPOSURE        15              /* picture exposure setting */
# Line 47 | Line 49 | int    onevalue(), catvalues();
49   #define PENUMBRAS       20              /* shadow penumbras are desired */
50   #define VARIABILITY     21              /* level of light variability */
51   #define REPORT          22              /* report frequency and errfile */
52 + #define RAWFILE         23              /* raw picture file root name */
53 + #define ZFILE           24              /* distance file root name */
54                                  /* total number of variables */
55 < #define NVARS           23
55 > #define NVARS           25
56  
57   VARIABLE        vv[NVARS] = {           /* variable-value pairs */
58          {"objects",     3,      0,      NULL,   catvalues},
# Line 61 | Line 65 | VARIABLE       vv[NVARS] = {           /* variable-value pairs */
65          {"pfilt",       2,      0,      NULL,   catvalues},
66          {"view",        2,      0,      NULL,   NULL},
67          {"ZONE",        2,      0,      NULL,   onevalue},
68 <        {"QUALITY",     3,      0,      NULL,   onevalue},
68 >        {"QUALITY",     3,      0,      NULL,   qualvalue},
69          {"OCTREE",      3,      0,      NULL,   onevalue},
70          {"PICTURE",     3,      0,      NULL,   onevalue},
71          {"AMBFILE",     3,      0,      NULL,   onevalue},
72          {"OPTFILE",     3,      0,      NULL,   onevalue},
73 <        {"EXPOSURE",    3,      0,      NULL,   onevalue},
73 >        {"EXPOSURE",    3,      0,      NULL,   fltvalue},
74          {"RESOLUTION",  3,      0,      NULL,   onevalue},
75          {"UP",          2,      0,      NULL,   onevalue},
76 <        {"INDIRECT",    3,      0,      NULL,   onevalue},
77 <        {"DETAIL",      3,      0,      NULL,   onevalue},
78 <        {"PENUMBRAS",   3,      0,      NULL,   onevalue},
79 <        {"VARIABILITY", 3,      0,      NULL,   onevalue},
76 >        {"INDIRECT",    3,      0,      NULL,   intvalue},
77 >        {"DETAIL",      3,      0,      NULL,   qualvalue},
78 >        {"PENUMBRAS",   3,      0,      NULL,   boolvalue},
79 >        {"VARIABILITY", 3,      0,      NULL,   qualvalue},
80          {"REPORT",      3,      0,      NULL,   onevalue},
81 +        {"RAWFILE",     3,      0,      NULL,   onevalue},
82 +        {"ZFILE",       2,      0,      NULL,   onevalue},
83   };
84  
85   VARIABLE        *matchvar();
86   char    *nvalue();
81 – int     vscale();
87  
88   #define UPPER(c)        ((c)&~0x20)     /* ASCII trick */
89  
# Line 87 | Line 92 | int    vscale();
92   #define vval(vc)        (vv[vc].value)
93   #define vint(vc)        atoi(vval(vc))
94   #define vlet(vc)        UPPER(vval(vc)[0])
95 + #define vscale          vlet
96   #define vbool(vc)       (vlet(vc)=='T')
97  
98 < #define HIGH            2
99 < #define MEDIUM          1
100 < #define LOW             0
98 > #define HIGH            'H'
99 > #define MEDIUM          'M'
100 > #define LOW             'L'
101  
96 – int     lowqopts(), medqopts(), hiqopts();
97 – int     (*setqopts[3])() = {lowqopts, medqopts, hiqopts};
98 –
99 – #define renderopts      (*setqopts[vscale(QUALITY)])
100 –
102                                  /* overture calculation file */
103   #ifdef NIX
104 < char    overfile[] = "overture.raw";
104 > char    overfile[] = "overture.unf";
105   #else
106   char    overfile[] = "/dev/null";
107   #endif
108  
109 < extern unsigned long    fdate(), time();
109 > extern time_t   fdate(), time();
110  
111 < unsigned long   scenedate;      /* date of latest scene or object file */
112 < unsigned long   octreedate;     /* date of octree */
113 < unsigned long   matdate;        /* date of latest material file */
114 < unsigned long   illumdate;      /* date of last illum file */
111 > time_t  scenedate;              /* date of latest scene or object file */
112 > time_t  octreedate;             /* date of octree */
113 > time_t  matdate;                /* date of latest material file */
114 > time_t  illumdate;              /* date of last illum file */
115  
116   char    *oct0name;              /* name of pre-mkillum octree */
117 < unsigned long   oct0date;       /* date of pre-mkillum octree */
117 > time_t  oct0date;               /* date of pre-mkillum octree */
118   char    *oct1name;              /* name of post-mkillum octree */
119 < unsigned long   oct1date;       /* date of post-mkillum octree (>= matdate) */
119 > time_t  oct1date;               /* date of post-mkillum octree (>= matdate) */
120  
121 + int     nowarn = 0;             /* no warnings */
122   int     explicate = 0;          /* explicate variables */
123   int     silent = 0;             /* do work silently */
124 + int     touchonly = 0;          /* touch files only */
125   int     noaction = 0;           /* don't do anything */
126   int     sayview = 0;            /* print view out */
127   char    *rvdevice = NULL;       /* rview output device */
# Line 137 | Line 140 | int    argc;
140   char    *argv[];
141   {
142          char    ropts[512];
143 +        char    popts[64];
144          int     i;
145  
146          progname = argv[0];
# Line 149 | Line 153 | char   *argv[];
153                  case 'n':
154                          noaction++;
155                          break;
156 +                case 't':
157 +                        touchonly++;
158 +                        break;
159                  case 'e':
160                          explicate++;
161                          break;
# Line 161 | Line 168 | char   *argv[];
168                  case 'v':
169                          viewselect = argv[++i];
170                          break;
171 +                case 'w':
172 +                        nowarn++;
173 +                        break;
174                  default:
175                          goto userr;
176                  }
# Line 188 | Line 198 | char   *argv[];
198                                  /* check date on ambient file */
199          checkambfile();
200                                  /* run simulation */
201 <        renderopts(ropts);
201 >        renderopts(ropts, popts);
202          xferopts(ropts);
203          if (rvdevice != NULL)
204 <                rview(ropts);
204 >                rview(ropts, popts);
205          else
206 <                rpict(ropts);
206 >                rpict(ropts, popts);
207          exit(0);
208   userr:
209          fprintf(stderr,
210 <        "Usage: %s [-s][-n][-e][-V][-v view][-o dev] rfile [VAR=value ..]\n",
210 > "Usage: %s [-s][-n][-t][-e][-V][-v view][-o dev] rfile [VAR=value ..]\n",
211                          progname);
212          exit(1);
213   }
# Line 256 | Line 266 | setvariable(ass)               /* assign variable according to stri
266   register char   *ass;
267   {
268          char    varname[32];
259 –        char    varval[512];
269          int     n;
270          register char   *cp;
271          register VARIABLE       *vp;
# Line 277 | Line 286 | register char  *ass;
286          for (n = strlen(ass); n > 0; n--)
287                  if (!isspace(ass[n-1]))
288                          break;
289 <        if (!n) {
289 >        if (!n && !nowarn) {
290                  fprintf(stderr, "%s: warning - missing value for variable '%s'\n",
291                                  progname, varname);
292                  return;
# Line 348 | Line 357 | register int   n;
357   }
358  
359  
351 – int
352 – vscale(vc)                      /* return scale for variable vc */
353 – int     vc;
354 – {
355 –        switch(vlet(vc)) {
356 –        case 'H':
357 –                return(HIGH);
358 –        case 'M':
359 –                return(MEDIUM);
360 –        case 'L':
361 –                return(LOW);
362 –        }
363 –        badvalue(vc);
364 – }
365 –
366 –
360   checkvalues()                   /* check assignments */
361   {
362          register int    i;
# Line 379 | Line 372 | register VARIABLE      *vp;
372   {
373          if (vp->nass < 2)
374                  return;
375 <        fprintf(stderr, "%s: warning - multiple assignment of variable '%s'\n",
375 >        if (!nowarn)
376 >                fprintf(stderr,
377 >                "%s: warning - multiple assignment of variable '%s'\n",
378                          progname, vp->name);
379          do
380                  vp->value += strlen(vp->value)+1;
# Line 402 | Line 397 | register VARIABLE      *vp;
397   }
398  
399  
400 < unsigned long
400 > int
401 > badmatch(tv, cv)                /* case insensitive truncated comparison */
402 > register char   *tv, *cv;
403 > {
404 >        if (!*tv) return(1);            /* null string cannot match */
405 >        do
406 >                if (UPPER(*tv) != *cv++)
407 >                        return(1);
408 >        while (*++tv);
409 >        return(0);              /* OK */
410 > }
411 >
412 >
413 > boolvalue(vp)                   /* check boolean for legal values */
414 > register VARIABLE       *vp;
415 > {
416 >        if (!vp->nass) return;
417 >        onevalue(vp);
418 >        switch (UPPER(vp->value[0])) {
419 >        case 'T':
420 >                if (badmatch(vp->value, "TRUE")) break;
421 >                return;
422 >        case 'F':
423 >                if (badmatch(vp->value, "FALSE")) break;
424 >                return;
425 >        }
426 >        fprintf(stderr, "%s: illegal value for boolean variable '%s'\n",
427 >                        progname, vp->name);
428 >        exit(1);
429 > }
430 >
431 >
432 > qualvalue(vp)                   /* check qualitative var. for legal values */
433 > register VARIABLE       *vp;
434 > {
435 >        if (!vp->nass) return;
436 >        onevalue(vp);
437 >        switch (UPPER(vp->value[0])) {
438 >        case 'L':
439 >                if (badmatch(vp->value, "LOW")) break;
440 >                return;
441 >        case 'M':
442 >                if (badmatch(vp->value, "MEDIUM")) break;
443 >                return;
444 >        case 'H':
445 >                if (badmatch(vp->value, "HIGH")) break;
446 >                return;
447 >        }
448 >        fprintf(stderr, "%s: illegal value for qualitative variable '%s'\n",
449 >                        progname, vp->name);
450 >        exit(1);
451 > }
452 >
453 >
454 > intvalue(vp)                    /* check integer variable for legal values */
455 > register VARIABLE       *vp;
456 > {
457 >        if (!vp->nass) return;
458 >        onevalue(vp);
459 >        if (isint(vp->value)) return;
460 >        fprintf(stderr, "%s: illegal value for integer variable '%s'\n",
461 >                        progname, vp->name);
462 >        exit(1);
463 > }
464 >
465 >
466 > fltvalue(vp)                    /* check float variable for legal values */
467 > register VARIABLE       *vp;
468 > {
469 >        if (!vp->nass) return;
470 >        onevalue(vp);
471 >        if (isflt(vp->value)) return;
472 >        fprintf(stderr, "%s: illegal value for real variable '%s'\n",
473 >                        progname, vp->name);
474 >        exit(1);
475 > }
476 >
477 >
478 > time_t
479   checklast(fnames)                       /* check files and find most recent */
480   register char   *fnames;
481   {
482          char    thisfile[MAXPATH];
483 <        unsigned long   thisdate, lastdate = 0;
483 >        time_t  thisdate, lastdate = 0;
484          register char   *cp;
485  
486          if (fnames == NULL)
# Line 455 | Line 528 | int    pred;
528  
529   checkfiles()                    /* check for existence and modified times */
530   {
531 <        unsigned long   objdate;
531 >        time_t  objdate;
532  
533          if (!vdef(OCTREE)) {
534                  if ((vval(OCTREE) = bmalloc(strlen(radname)+5)) == NULL)
# Line 490 | Line 563 | double org[3], *sizp;
563          extern FILE     *popen();
564          static double   oorg[3], osiz = 0.;
565          double  min[3], max[3];
566 <        char    buf[512];
566 >        char    buf[1024];
567          FILE    *fp;
568          register int    i;
569  
# Line 583 | Line 656 | setdefaults()                  /* set default values for unassigned v
656  
657   printvals()                     /* print variable values */
658   {
659 <        register int    i, j;
659 >        int     i, j, clipline;
660 >        register char   *cp;
661 >        register int    k;
662  
663 <        for (i = 0; i < NVARS; i++)
664 <                for (j = 0; j < vdef(i); j++)
665 <                        printf("%s= %s\n", vnam(i), nvalue(vv+i, j));
663 >        for (i = 0; i < NVARS; i++)             /* print each variable */
664 >            for (j = 0; j < vdef(i); j++) {     /* print each assignment */
665 >                fputs(vnam(i), stdout);
666 >                fputs("= ", stdout);
667 >                k = clipline = ( vv[i].fixval == catvalues ? 64 : 320 )
668 >                                - strlen(vnam(i)) ;
669 >                cp = nvalue(vv+i, j);
670 >                while (*cp) {
671 >                    putchar(*cp++);
672 >                    if (--k <= 0) {             /* line too long */
673 >                        while (*cp && !isspace(*cp))
674 >                            putchar(*cp++);     /* finish this word */
675 >                        if (*cp) {              /* start new line */
676 >                            putchar('\n');
677 >                            fputs(vnam(i), stdout);
678 >                            putchar('=');
679 >                            k = clipline;
680 >                        }
681 >                    }
682 >                }
683 >                putchar('\n');
684 >            }
685          fflush(stdout);
686   }
687  
# Line 595 | Line 689 | printvals()                    /* print variable values */
689   oconv()                         /* run oconv and mkillum if necessary */
690   {
691          static char     illumtmp[] = "ilXXXXXX";
692 <        char    combuf[512], ocopts[64], mkopts[64];
692 >        char    combuf[1024], ocopts[64], mkopts[64];
693  
694          oconvopts(ocopts);              /* get options */
695          if (octreedate < scenedate) {   /* check date on original octree */
696 <                                                /* build command */
697 <                if (vdef(MATERIAL))
698 <                        sprintf(combuf, "oconv%s %s %s > %s", ocopts,
699 <                                vval(MATERIAL), vval(SCENE), vval(OCTREE));
700 <                else
701 <                        sprintf(combuf, "oconv%s %s > %s", ocopts,
702 <                                        vval(SCENE), vval(OCTREE));
703 <                
704 <                if (runcom(combuf)) {           /* run it */
705 <                        fprintf(stderr,
696 >                if (touchonly && octreedate)
697 >                        touch(vval(OCTREE));
698 >                else {                          /* build command */
699 >                        if (vdef(MATERIAL))
700 >                                sprintf(combuf, "oconv%s %s %s > %s", ocopts,
701 >                                                vval(MATERIAL), vval(SCENE),
702 >                                                vval(OCTREE));
703 >                        else
704 >                                sprintf(combuf, "oconv%s %s > %s", ocopts,
705 >                                                vval(SCENE), vval(OCTREE));
706 >                        
707 >                        if (runcom(combuf)) {           /* run it */
708 >                                fprintf(stderr,
709                                  "%s: error generating octree\n\t%s removed\n",
710 <                                        progname, vval(OCTREE));
711 <                        unlink(vval(OCTREE));
712 <                        exit(1);
710 >                                                progname, vval(OCTREE));
711 >                                unlink(vval(OCTREE));
712 >                                exit(1);
713 >                        }
714                  }
715 <                octreedate = time(0);
715 >                octreedate = time((time_t *)NULL);
716 >                if (octreedate < scenedate)     /* in case clock is off */
717 >                        octreedate = scenedate;
718          }
719          if (oct1name == vval(OCTREE))           /* no mkillum? */
720                  oct1date = octreedate > matdate ? octreedate : matdate;
# Line 623 | Line 723 | oconv()                                /* run oconv and mkillum if necessary */
723                  return;
724                                                  /* make octree0 */
725          if (oct0date < scenedate | oct0date < illumdate) {
726 <                                                /* build command */
726 >                if (touchonly && oct0date)
727 >                        touch(oct0name);
728 >                else {                          /* build command */
729 >                        if (octreedate)
730 >                                sprintf(combuf, "oconv%s -i %s %s > %s", ocopts,
731 >                                        vval(OCTREE), vval(ILLUM), oct0name);
732 >                        else if (vdef(MATERIAL))
733 >                                sprintf(combuf, "oconv%s %s %s > %s", ocopts,
734 >                                        vval(MATERIAL), vval(ILLUM), oct0name);
735 >                        else
736 >                                sprintf(combuf, "oconv%s %s > %s", ocopts,
737 >                                        vval(ILLUM), oct0name);
738 >                        if (runcom(combuf)) {           /* run it */
739 >                                fprintf(stderr,
740 >                                "%s: error generating octree\n\t%s removed\n",
741 >                                                progname, oct0name);
742 >                                unlink(oct0name);
743 >                                exit(1);
744 >                        }
745 >                }
746 >                oct0date = time((time_t *)NULL);
747 >                if (oct0date < octreedate)      /* in case clock is off */
748 >                        oct0date = octreedate;
749 >                if (oct0date < illumdate)       /* ditto */
750 >                        oct0date = illumdate;
751 >                }
752 >        if (touchonly && oct1date)
753 >                touch(oct1name);
754 >        else {
755 >                mkillumopts(mkopts);            /* build mkillum command */
756 >                mktemp(illumtmp);
757 >                sprintf(combuf, "mkillum%s %s \"<\" %s > %s", mkopts,
758 >                                oct0name, vval(ILLUM), illumtmp);
759 >                if (runcom(combuf)) {                   /* run it */
760 >                        fprintf(stderr, "%s: error running mkillum\n",
761 >                                        progname);
762 >                        unlink(illumtmp);
763 >                        exit(1);
764 >                }
765 >                                                /* make octree1 (frozen) */
766                  if (octreedate)
767 <                        sprintf(combuf, "oconv%s -i %s %s > %s", ocopts,
768 <                                vval(OCTREE), vval(ILLUM), oct0name);
767 >                        sprintf(combuf, "oconv%s -f -i %s %s > %s", ocopts,
768 >                                vval(OCTREE), illumtmp, oct1name);
769                  else if (vdef(MATERIAL))
770 <                        sprintf(combuf, "oconv%s %s %s > %s", ocopts,
771 <                                vval(MATERIAL), vval(ILLUM), oct0name);
770 >                        sprintf(combuf, "oconv%s -f %s %s > %s", ocopts,
771 >                                vval(MATERIAL), illumtmp, oct1name);
772                  else
773 <                        sprintf(combuf, "oconv%s %s > %s", ocopts,
774 <                                vval(ILLUM), oct0name);
773 >                        sprintf(combuf, "oconv%s -f %s > %s", ocopts,
774 >                                illumtmp, oct1name);
775                  if (runcom(combuf)) {           /* run it */
776                          fprintf(stderr,
777                                  "%s: error generating octree\n\t%s removed\n",
778 <                                        progname, oct0name);
779 <                        unlink(oct0name);
778 >                                        progname, oct1name);
779 >                        unlink(oct1name);
780 >                        unlink(illumtmp);
781                          exit(1);
782                  }
783 <                oct0date = time(0);
783 >                rmfile(illumtmp);
784          }
785 <        mkillumopts(mkopts);                    /* build mkillum command */
786 <        mktemp(illumtmp);
787 <        sprintf(combuf, "mkillum%s %s \"<\" %s > %s", mkopts,
648 <                        oct0name, vval(ILLUM), illumtmp);
649 <        if (runcom(combuf)) {                   /* run it */
650 <                fprintf(stderr, "%s: error running mkillum\n", progname);
651 <                unlink(illumtmp);
652 <                exit(1);
653 <        }
654 <                                                /* make octree1 (frozen) */
655 <        if (octreedate)
656 <                sprintf(combuf, "oconv%s -f -i %s %s > %s", ocopts,
657 <                        vval(OCTREE), illumtmp, oct1name);
658 <        else if (vdef(MATERIAL))
659 <                sprintf(combuf, "oconv%s -f %s %s > %s", ocopts,
660 <                        vval(MATERIAL), illumtmp, oct1name);
661 <        else
662 <                sprintf(combuf, "oconv%s -f %s > %s", ocopts,
663 <                        illumtmp, oct1name);
664 <        if (runcom(combuf)) {           /* run it */
665 <                fprintf(stderr,
666 <                        "%s: error generating octree\n\t%s removed\n",
667 <                                progname, oct1name);
668 <                unlink(oct1name);
669 <                exit(1);
670 <        }
671 <        oct1date = time(0);
672 <        rmfile(illumtmp);
785 >        oct1date = time((time_t *)NULL);
786 >        if (oct1date < oct0date)        /* in case clock is off */
787 >                oct1date = oct0date;
788   }
789  
790  
# Line 708 | Line 823 | register char  *mo;
823  
824   checkambfile()                  /* check date on ambient file */
825   {
826 <        unsigned long   afdate;
826 >        time_t  afdate;
827  
828          if (!vdef(AMBFILE))
829                  return;
830          if (!(afdate = fdate(vval(AMBFILE))))
831                  return;
832          if (oct1date > afdate)
833 <                rmfile(vval(AMBFILE));
833 >                if (touchonly)
834 >                        touch(vval(AMBFILE));
835 >                else
836 >                        rmfile(vval(AMBFILE));
837   }
838  
839  
# Line 723 | Line 841 | double
841   ambval()                                /* compute ambient value */
842   {
843          if (vdef(EXPOSURE)) {
726 –                if (!isflt(vval(EXPOSURE)))
727 –                        badvalue(EXPOSURE);
844                  if (vval(EXPOSURE)[0] == '+' || vval(EXPOSURE)[0] == '-')
845                          return(.5/pow(2.,atof(vval(EXPOSURE))));
846                  return(.5/atof(vval(EXPOSURE)));
# Line 737 | Line 853 | ambval()                               /* compute ambient value */
853   }
854  
855  
856 < lowqopts(op)                            /* low quality rendering options */
856 > renderopts(op, po)                      /* set rendering options */
857 > char    *op, *po;
858 > {
859 >        switch(vscale(QUALITY)) {
860 >        case LOW:
861 >                lowqopts(op, po);
862 >                break;
863 >        case MEDIUM:
864 >                medqopts(op, po);
865 >                break;
866 >        case HIGH:
867 >                hiqopts(op, po);
868 >                break;
869 >        }
870 > }
871 >
872 >
873 > lowqopts(op, po)                        /* low quality rendering options */
874   register char   *op;
875 + char    *po;
876   {
877          double  d, org[3], siz[3];
878  
879          *op = '\0';
880 +        *po = '\0';
881          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
882                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
883                  badvalue(ZONE);
884          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
885 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
886 +                badvalue(ZONE);
887          getoctcube(org, &d);
888          d *= 3./(siz[0]+siz[1]+siz[2]);
889          switch (vscale(DETAIL)) {
890          case LOW:
891 <                op = addarg(op, "-ps 16 -dp 64");
891 >                po = addarg(po, "-ps 16");
892 >                op = addarg(op, "-dp 64");
893                  sprintf(op, " -ar %d", (int)(4*d));
894                  op += strlen(op);
895                  break;
896          case MEDIUM:
897 <                op = addarg(op, "-ps 8 -dp 128");
897 >                po = addarg(po, "-ps 8");
898 >                op = addarg(op, "-dp 128");
899                  sprintf(op, " -ar %d", (int)(8*d));
900                  op += strlen(op);
901                  break;
902          case HIGH:
903 <                op = addarg(op, "-ps 4 -dp 256");
903 >                po = addarg(po, "-ps 4");
904 >                op = addarg(op, "-dp 256");
905                  sprintf(op, " -ar %d", (int)(16*d));
906                  op += strlen(op);
907                  break;
908          }
909 <        op = addarg(op, "-pt .16");
909 >        po = addarg(po, "-pt .16");
910          if (vbool(PENUMBRAS))
911                  op = addarg(op, "-ds .4");
912          else
# Line 798 | Line 938 | register char  *op;
938   }
939  
940  
941 < medqopts(op)                            /* medium quality rendering options */
941 > medqopts(op, po)                        /* medium quality rendering options */
942   register char   *op;
943 + char    *po;
944   {
945 <        double  d, org[3], siz[3];
945 >        double  d, org[3], siz[3], asz;
946  
947          *op = '\0';
948 +        *po = '\0';
949          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
950                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
951                  badvalue(ZONE);
952          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
953 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
954 +                badvalue(ZONE);
955          getoctcube(org, &d);
956 <        d *= 3./(siz[0]+siz[1]+siz[2]);
956 >        asz = (siz[0]+siz[1]+siz[2])/3.;
957 >        d /= asz;
958          switch (vscale(DETAIL)) {
959          case LOW:
960 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
960 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
961                  op = addarg(op, "-dp 256");
962                  sprintf(op, " -ar %d", (int)(8*d));
963                  op += strlen(op);
964 +                sprintf(op, " -ms %.2g", asz/20.);
965 +                op += strlen(op);
966                  break;
967          case MEDIUM:
968 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
968 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
969                  op = addarg(op, "-dp 512");
970                  sprintf(op, " -ar %d", (int)(16*d));
971                  op += strlen(op);
972 +                sprintf(op, " -ms %.2g", asz/40.);
973 +                op += strlen(op);
974                  break;
975          case HIGH:
976 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
976 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
977                  op = addarg(op, "-dp 1024");
978                  sprintf(op, " -ar %d", (int)(32*d));
979                  op += strlen(op);
980 +                sprintf(op, " -ms %.2g", asz/80.);
981 +                op += strlen(op);
982                  break;
983          }
984 <        op = addarg(op, "-pt .08");
984 >        po = addarg(po, "-pt .08");
985          if (vbool(PENUMBRAS))
986 <                op = addarg(op, "-ds .2 -dj .35");
986 >                op = addarg(op, "-ds .2 -dj .5");
987          else
988                  op = addarg(op, "-ds .3");
989          op = addarg(op, "-dt .1 -dc .5 -dr 1 -sj .7 -st .1");
# Line 865 | Line 1016 | register char  *op;
1016   }
1017  
1018  
1019 < hiqopts(op)                             /* high quality rendering options */
1019 > hiqopts(op, po)                         /* high quality rendering options */
1020   register char   *op;
1021 + char    *po;
1022   {
1023 <        double  d, org[3], siz[3];
1023 >        double  d, org[3], siz[3], asz;
1024  
1025          *op = '\0';
1026 +        *po = '\0';
1027          if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
1028                          &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
1029                  badvalue(ZONE);
1030          siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
1031 +        if (siz[0] <= FTINY | siz[1] <= FTINY | siz[2] <= FTINY)
1032 +                badvalue(ZONE);
1033          getoctcube(org, &d);
1034 <        d *= 3./(siz[0]+siz[1]+siz[2]);
1034 >        asz = (siz[0]+siz[1]+siz[2])/3.;
1035 >        d /= asz;
1036          switch (vscale(DETAIL)) {
1037          case LOW:
1038 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
1038 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
1039                  op = addarg(op, "-dp 1024");
1040                  sprintf(op, " -ar %d", (int)(16*d));
1041                  op += strlen(op);
1042 +                sprintf(op, " -ms %.2g", asz/40.);
1043 +                op += strlen(op);
1044                  break;
1045          case MEDIUM:
1046 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
1046 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
1047                  op = addarg(op, "-dp 2048");
1048                  sprintf(op, " -ar %d", (int)(32*d));
1049                  op += strlen(op);
1050 +                sprintf(op, " -ms %.2g", asz/80.);
1051 +                op += strlen(op);
1052                  break;
1053          case HIGH:
1054 <                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
1054 >                po = addarg(po, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
1055                  op = addarg(op, "-dp 4096");
1056                  sprintf(op, " -ar %d", (int)(64*d));
1057                  op += strlen(op);
1058 +                sprintf(op, " -ms %.2g", asz/160.);
1059 +                op += strlen(op);
1060                  break;
1061          }
1062 <        op = addarg(op, "-pt .04");
1062 >        po = addarg(po, "-pt .04");
1063          if (vbool(PENUMBRAS))
1064 <                op = addarg(op, "-ds .1 -dj .7");
1064 >                op = addarg(op, "-ds .1 -dj .65");
1065          else
1066                  op = addarg(op, "-ds .2");
1067          op = addarg(op, "-dt .05 -dc .75 -dr 3 -sj 1 -st .01");
# Line 1036 | Line 1198 | register char  *vs;
1198                  zpos = -1; vs++;
1199          }
1200          viewtype = 'v';
1201 <        if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h')
1201 >        if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h' | *vs == 'c')
1202                  viewtype = *vs++;
1203          cp = viewopts;
1204          if ((!*vs || isspace(*vs)) && (xpos|ypos|zpos)) {       /* got one! */
# Line 1047 | Line 1209 | register char  *vs;
1209                          badvalue(ZONE);
1210                  for (i = 0; i < 3; i++) {
1211                          dim[i] -= cent[i];
1212 +                        if (dim[i] <= FTINY)
1213 +                                badvalue(ZONE);
1214                          cent[i] += .5*dim[i];
1215                  }
1216                  mult = vlet(ZONE)=='E' ? 2. : .45 ;
# Line 1089 | Line 1253 | register char  *vs;
1253                  case 'h':
1254                          cp = addarg(cp, "-vh 180 -vv 180");
1255                          break;
1256 +                case 'c':
1257 +                        cp = addarg(cp, "-vh 180 -vv 90");
1258 +                        break;
1259                  }
1260          } else {
1261                  while (!isspace(*vs))           /* else skip id */
# Line 1210 | Line 1377 | register char  *vopts;
1377   }
1378  
1379  
1380 < rview(opts)                             /* run rview with first view */
1381 < char    *opts;
1380 > rview(opts, po)                         /* run rview with first view */
1381 > char    *opts, *po;
1382   {
1383          char    *vw;
1384          char    combuf[512];
1385                                          /* build command */
1386 <        if ((vw = getview(0, NULL)) == NULL)
1386 >        if (touchonly || (vw = getview(0, NULL)) == NULL)
1387                  return;
1388          if (sayview)
1389                  printview(vw);
1390 <        sprintf(combuf, "rview %s%s -R %s ", vw, opts, rifname);
1390 >        sprintf(combuf, "rview %s%s%s -R %s ", vw, po, opts, rifname);
1391          if (rvdevice != NULL)
1392                  sprintf(combuf+strlen(combuf), "-o %s ", rvdevice);
1393 +        if (vdef(EXPOSURE))
1394 +                sprintf(combuf+strlen(combuf), "-pe %s ", vval(EXPOSURE));
1395          strcat(combuf, oct1name);
1396          if (runcom(combuf)) {           /* run it */
1397                  fprintf(stderr, "%s: error running rview\n", progname);
# Line 1231 | Line 1400 | char   *opts;
1400   }
1401  
1402  
1403 < rpict(opts)                             /* run rpict and pfilt for each view */
1404 < char    *opts;
1403 > rpict(opts, po)                         /* run rpict and pfilt for each view */
1404 > char    *opts, *po;
1405   {
1406          char    combuf[1024];
1407 <        char    rawfile[MAXPATH], picfile[MAXPATH], rep[MAXPATH+16], res[32];
1407 >        char    rawfile[MAXPATH], picfile[MAXPATH];
1408 >        char    zopt[MAXPATH+4], rep[MAXPATH+16], res[32];
1409          char    pfopts[128];
1410          char    vs[32], *vw;
1411          int     vn, mult;
1412 +        time_t  rfdt, pfdt;
1413                                          /* get pfilt options */
1414          pfiltopts(pfopts);
1415                                          /* get resolution, reporting */
1416 <        mult = vscale(QUALITY)+1;
1416 >        switch (vscale(QUALITY)) {
1417 >        case LOW:
1418 >                mult = 1;
1419 >                break;
1420 >        case MEDIUM:
1421 >                mult = 2;
1422 >                break;
1423 >        case HIGH:
1424 >                mult = 3;
1425 >                break;
1426 >        }
1427          {
1428                  int     xres, yres;
1429                  double  aspect;
# Line 1277 | Line 1458 | char   *opts;
1458                  if (!vs[0])
1459                          sprintf(vs, "%d", vn);
1460                  sprintf(picfile, "%s_%s.pic", vval(PICTURE), vs);
1461 +                if (vdef(ZFILE))
1462 +                        sprintf(zopt, " -z %s_%s.zbf", vval(ZFILE), vs);
1463 +                else
1464 +                        zopt[0] = '\0';
1465                                                  /* check date on picture */
1466 <                if (fdate(picfile) >= oct1date)
1466 >                pfdt = fdate(picfile);
1467 >                if (pfdt >= oct1date)
1468                          continue;
1469 +                                                /* get raw file name */
1470 +                sprintf(rawfile, "%s_%s.unf",
1471 +                        vdef(RAWFILE) ? vval(RAWFILE) : vval(PICTURE), vs);
1472 +                rfdt = fdate(rawfile);
1473 +                if (touchonly) {                /* update times only */
1474 +                        if (rfdt) {
1475 +                                if (rfdt < oct1date)
1476 +                                        touch(rawfile);
1477 +                        } else if (pfdt && pfdt < oct1date)
1478 +                                touch(picfile);
1479 +                        continue;
1480 +                }
1481                                                  /* build rpict command */
1482 <                sprintf(rawfile, "%s_%s.raw", vval(PICTURE), vs);
1483 <                if (fdate(rawfile) >= oct1date)         /* recover */
1484 <                        sprintf(combuf, "rpict%s%s -ro %s %s",
1287 <                                        rep, opts, rawfile, oct1name);
1482 >                if (rfdt >= oct1date)           /* recover */
1483 >                        sprintf(combuf, "rpict%s%s%s%s -ro %s %s",
1484 >                                rep, po, opts, zopt, rawfile, oct1name);
1485                  else {
1486                          if (overture) {         /* run overture calculation */
1487                                  sprintf(combuf,
# Line 1301 | Line 1498 | char   *opts;
1498                                  rmfile(overfile);
1499   #endif
1500                          }
1501 <                        sprintf(combuf, "rpict%s %s %s%s %s > %s",
1502 <                                        rep, vw, res, opts,
1501 >                        sprintf(combuf, "rpict%s %s %s%s%s%s %s > %s",
1502 >                                        rep, vw, res, po, opts, zopt,
1503                                          oct1name, rawfile);
1504                  }
1505                  if (runcom(combuf)) {           /* run rpict */
# Line 1310 | Line 1507 | char   *opts;
1507                                          progname, vs);
1508                          exit(1);
1509                  }
1510 +                if (!vdef(RAWFILE) || strcmp(vval(RAWFILE),vval(PICTURE))) {
1511                                                  /* build pfilt command */
1512 <                if (mult > 1)
1513 <                        sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
1512 >                        if (mult > 1)
1513 >                                sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
1514                                          pfopts, mult, mult, rawfile, picfile);
1515 <                else
1516 <                        sprintf(combuf, "pfilt%s %s > %s", pfopts,
1517 <                                        rawfile, picfile);
1518 <                if (runcom(combuf)) {           /* run pfilt */
1519 <                        fprintf(stderr,
1520 <                        "%s: error filtering view %s\n\t%s removed\n",
1521 <                                        progname, vs, picfile);
1522 <                        unlink(picfile);
1523 <                        exit(1);
1515 >                        else
1516 >                                sprintf(combuf, "pfilt%s %s > %s", pfopts,
1517 >                                                rawfile, picfile);
1518 >                        if (runcom(combuf)) {           /* run pfilt */
1519 >                                fprintf(stderr,
1520 >                                "%s: error filtering view %s\n\t%s removed\n",
1521 >                                                progname, vs, picfile);
1522 >                                unlink(picfile);
1523 >                                exit(1);
1524 >                        }
1525                  }
1526 <                                                /* remove raw file */
1527 <                rmfile(rawfile);
1526 >                                                /* remove/rename raw file */
1527 >                if (vdef(RAWFILE)) {
1528 >                        sprintf(combuf, "%s_%s.pic", vval(RAWFILE), vs);
1529 >                        mvfile(rawfile, combuf);
1530 >                } else
1531 >                        rmfile(rawfile);
1532          }
1533   }
1534  
1535  
1536 + touch(fn)                       /* update a file */
1537 + char    *fn;
1538 + {
1539 +        if (!silent)
1540 +                printf("\ttouch %s\n", fn);
1541 +        if (noaction)
1542 +                return(0);
1543 + #ifdef notused
1544 +        if (access(fn, F_OK) == -1)             /* create it */
1545 +                if (close(open(fn, O_WRONLY|O_CREAT, 0666)) == -1)
1546 +                        return(-1);
1547 + #endif
1548 +        return(setfdate(fn, time((time_t *)NULL)));
1549 + }
1550 +
1551 +
1552   runcom(cs)                      /* run command */
1553   char    *cs;
1554   {
# Line 1354 | Line 1573 | char   *fn;
1573          if (noaction)
1574                  return(0);
1575          return(unlink(fn));
1576 + }
1577 +
1578 +
1579 + mvfile(fold, fnew)              /* move a file */
1580 + char    *fold, *fnew;
1581 + {
1582 +        if (!silent)
1583 + #ifdef MSDOS
1584 +                printf("\trename %s %s\n", fold, fnew);
1585 + #else
1586 +                printf("\tmv %s %s\n", fold, fnew);
1587 + #endif
1588 +        if (noaction)
1589 +                return(0);
1590 +        return(rename(fold, fnew));
1591   }
1592  
1593  

Diff Legend

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