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.111 by greg, Sun Oct 26 17:35:53 2014 UTC vs.
Revision 2.114 by greg, Wed May 27 13:29:57 2015 UTC

# Line 36 | Line 36 | static const char      RCSid[] = "$Id$";
36   #define INDIRECT        5               /* indirection in lighting */
37   #define MATERIAL        6               /* material files */
38   #define MKILLUM         7               /* mkillum options */
39 < #define OBJECT          8               /* object files */
40 < #define OCONV           9               /* oconv options */
41 < #define OCTREE          10              /* octree file name */
42 < #define OPTFILE         11              /* rendering options file */
43 < #define PENUMBRAS       12              /* shadow penumbras are desired */
44 < #define PFILT           13              /* pfilt options */
45 < #define PICTURE         14              /* picture file root name */
46 < #define QUALITY         15              /* desired rendering quality */
47 < #define RAWFILE         16              /* raw picture file root name */
48 < #define RENDER          17              /* rendering options */
49 < #define REPORT          18              /* report frequency and errfile */
50 < #define RESOLUTION      19              /* maximum picture resolution */
51 < #define RPICT           20              /* rpict parameters */
52 < #define RVU             21              /* rvu parameters */
53 < #define SCENE           22              /* scene files */
54 < #define UP              23              /* view up (X, Y or Z) */
55 < #define VARIABILITY     24              /* level of light variability */
56 < #define VIEWS           25              /* view(s) for picture(s) */
57 < #define ZFILE           26              /* distance file root name */
58 < #define ZONE            27              /* simulation zone */
39 > #define MKPMAP          8               /* mkpmap options */
40 > #define OBJECT          9               /* object files */
41 > #define OCONV           10              /* oconv options */
42 > #define OCTREE          11              /* octree file name */
43 > #define OPTFILE         12              /* rendering options file */
44 > #define PCMAP           13              /* caustic photon map */
45 > #define PENUMBRAS       14              /* shadow penumbras are desired */
46 > #define PFILT           15              /* pfilt options */
47 > #define PGMAP           16              /* global photon map */
48 > #define PICTURE         17              /* picture file root name */
49 > #define QUALITY         18              /* desired rendering quality */
50 > #define RAWFILE         19              /* raw picture file root name */
51 > #define RENDER          20              /* rendering options */
52 > #define REPORT          21              /* report frequency and errfile */
53 > #define RESOLUTION      22              /* maximum picture resolution */
54 > #define RPICT           23              /* rpict parameters */
55 > #define RVU             24              /* rvu parameters */
56 > #define SCENE           25              /* scene files */
57 > #define UP              26              /* view up (X, Y or Z) */
58 > #define VARIABILITY     27              /* level of light variability */
59 > #define VIEWS           28              /* view(s) for picture(s) */
60 > #define ZFILE           29              /* distance file root name */
61 > #define ZONE            30              /* simulation zone */
62                                  /* total number of variables */
63 < int NVARS = 28;
63 > int NVARS = 31;
64  
65   VARIABLE        vv[] = {                /* variable-value pairs */
66          {"AMBFILE",     3,      0,      NULL,   onevalue},
# Line 68 | Line 71 | VARIABLE       vv[] = {                /* variable-value pairs */
71          {"INDIRECT",    3,      0,      NULL,   intvalue},
72          {"materials",   3,      0,      NULL,   catvalues},
73          {"mkillum",     3,      0,      NULL,   catvalues},
74 +        {"mkpmap",      3,      0,      NULL,   catvalues},
75          {"objects",     3,      0,      NULL,   catvalues},
76          {"oconv",       3,      0,      NULL,   catvalues},
77          {"OCTREE",      3,      0,      NULL,   onevalue},
78          {"OPTFILE",     3,      0,      NULL,   onevalue},
79 +        {"PCMAP",       2,      0,      NULL,   onevalue},
80          {"PENUMBRAS",   3,      0,      NULL,   boolvalue},
81          {"pfilt",       2,      0,      NULL,   catvalues},
82 +        {"PGMAP",       2,      0,      NULL,   onevalue},
83          {"PICTURE",     3,      0,      NULL,   onevalue},
84          {"QUALITY",     3,      0,      NULL,   qualvalue},
85          {"RAWFILE",     3,      0,      NULL,   onevalue},
# Line 108 | Line 114 | time_t oct0date;               /* date of pre-mkillum octree */
114   char    *oct1name;              /* name of post-mkillum octree */
115   time_t  oct1date;               /* date of post-mkillum octree (>= matdate) */
116  
117 + char    *pgmapname;             /* name of global photon map */
118 + time_t  pgmapdate;              /* date of global photon map (>= oct1date) */
119 + char    *pcmapname;             /* name of caustic photon map */
120 + time_t  pcmapdate;              /* date of caustic photon map (>= oct1date) */
121 +
122   int     nowarn = 0;             /* no warnings */
123   int     explicate = 0;          /* explicate variables */
124   int     silent = 0;             /* do work silently */
# Line 122 | Line 133 | char   *viewselect = NULL;     /* specific view only */
133                                  /* command paths */
134   char    c_oconv[256] = "oconv";
135   char    c_mkillum[256] = "mkillum";
136 + char    c_mkpmap[256] = "mkpmap";
137   char    c_rvu[256] = "rvu";
138   char    c_rpict[256] = DEF_RPICT_PATH;
139   char    c_rpiece[] = "rpiece";
# Line 145 | Line 157 | static void checkfiles(void);
157   static void getoctcube(double   org[3], double  *sizp);
158   static void setdefaults(void);
159   static void oconv(void);
160 + static void mkpmap(void);
161   static char * addarg(char       *op, char       *arg);
162   static void oconvopts(char      *oo);
163   static void mkillumopts(char    *mo);
164 + static void mkpmapopts(char     *mo);
165   static void checkambfile(void);
166   static double ambval(void);
167   static void renderopts(char     *op, char       *po);
# Line 244 | Line 258 | main(
258                  printvars(stdout);
259                                  /* build octree (and run mkillum) */
260          oconv();
261 +                                /* run mkpmap if indicated */
262 +        mkpmap();
263                                  /* check date on ambient file */
264          checkambfile();
265                                  /* run simulation */
# Line 336 | Line 352 | newfname(              /* create modified file name */
352   static void
353   checkfiles(void)                        /* check for existence and modified times */
354   {
355 +        char    fntemp[256];
356          time_t  objdate;
357  
358          if (!vdef(OCTREE)) {
# Line 365 | Line 382 | checkfiles(void)                       /* check for existence and modified
382                                  vnam(OCTREE), vnam(SCENE), vnam(ILLUM));
383                  quit(1);
384          }
385 +        if (vdef(PGMAP)) {
386 +                if (!*sskip2(vval(PGMAP),1)) {
387 +                        fprintf(stderr, "%s: '%s' missing # photons argument\n",
388 +                                        progname, vnam(PGMAP));
389 +                        quit(1);
390 +                }
391 +                atos(fntemp, sizeof(fntemp), vval(PGMAP));
392 +                pgmapname = savqstr(fntemp);
393 +                pgmapdate = fdate(pgmapname);
394 +        }
395 +        if (vdef(PCMAP)) {
396 +                if (!*sskip2(vval(PCMAP),1)) {
397 +                        fprintf(stderr, "%s: '%s' missing # photons argument\n",
398 +                                        progname, vnam(PCMAP));
399 +                        quit(1);
400 +                }
401 +                atos(fntemp, sizeof(fntemp), vval(PCMAP));
402 +                pcmapname = savqstr(fntemp);
403 +                pcmapdate = fdate(pcmapname);
404 +        }
405          matdate = checklast(vval(MATERIAL));
406   }      
407  
# Line 586 | Line 623 | oconv(void)                            /* run oconv and mkillum if necessary *
623   }
624  
625  
626 + static void
627 + mkpmap(void)                    /* run mkpmap if indicated */
628 + {
629 +        char    combuf[2048], *cp;
630 +        time_t  tnow;
631 +                                /* nothing to do? */
632 +        if ((pgmapname == NULL) | (pgmapdate >= oct1date) &&
633 +                        (pcmapname == NULL) | (pcmapdate >= oct1date))
634 +                return;
635 +                                /* just update existing file dates? */
636 +        if (touchonly && (pgmapname == NULL) | (pgmapdate > 0) &&
637 +                        (pcmapname == NULL) | (pcmapdate > 0)) {
638 +                if (pgmapname != NULL)
639 +                        touch(pgmapname);
640 +                if (pcmapname != NULL)
641 +                        touch(pcmapname);
642 +        } else {                /* else need to (re)run pkpmap */
643 +                strcpy(combuf, c_mkpmap);
644 +                for (cp = combuf; *cp; cp++)
645 +                        ;
646 +                mkpmapopts(cp);
647 +                if (vdef(REPORT)) {
648 +                        char    errfile[256];
649 +                        int     n;
650 +                        double  minutes;
651 +                        n = sscanf(vval(REPORT), "%lf %s", &minutes, errfile);
652 +                        if (n == 2)
653 +                                sprintf(cp, " -t %d -e %s", (int)(minutes*60), errfile);
654 +                        else if (n == 1)
655 +                                sprintf(cp, " -t %d", (int)(minutes*60));
656 +                        else
657 +                                badvalue(REPORT);
658 +                }
659 +                if (pgmapname != NULL && pgmapdate < oct1date) {
660 +                        cp = addarg(cp, "-apg");
661 +                        addarg(cp, vval(PGMAP));
662 +                        cp = sskip(sskip(cp));  /* remove any bandwidth */
663 +                        *cp = '\0';
664 +                }
665 +                if (pcmapname != NULL && pcmapdate < oct1date) {
666 +                        cp = addarg(cp, "-apc");
667 +                        addarg(cp, vval(PCMAP));
668 +                        cp = sskip(sskip(cp));  /* remove any bandwidth */
669 +                        *cp = '\0';
670 +                }
671 +                cp = addarg(cp, oct1name);
672 +                if (runcom(combuf)) {
673 +                        fprintf(stderr, "%s: error running %s\n",
674 +                                        progname, c_mkpmap);
675 +                        if (pgmapname != NULL)
676 +                                unlink(pgmapname);
677 +                        if (pcmapname != NULL)
678 +                                unlink(pcmapname);
679 +                        quit(1);
680 +                }
681 +        }
682 +        tnow = time((time_t *)NULL);
683 +        if (pgmapname != NULL)
684 +                pgmapdate = tnow;
685 +        if (pcmapname != NULL)
686 +                pcmapdate = tnow;
687 +        oct1date = tnow;        /* trigger ambient file removal if needed */
688 + }
689 +
690 +
691   static char *
692   addarg(                         /* append argument and advance pointer */
693   char    *op,
# Line 609 | Line 711 | oconvopts(                             /* get oconv options */
711          /* BEWARE:  This may be called via setdefaults(), so no assumptions */
712  
713          *oo = '\0';
714 <        if (vdef(OCONV))
715 <                if (vval(OCONV)[0] != '-') {
716 <                        atos(c_oconv, sizeof(c_oconv), vval(OCONV));
717 <                        oo = addarg(oo, sskip2(vval(OCONV), 1));
718 <                } else
719 <                        oo = addarg(oo, vval(OCONV));
714 >        if (!vdef(OCONV))
715 >                return;
716 >        if (vval(OCONV)[0] != '-') {
717 >                atos(c_oconv, sizeof(c_oconv), vval(OCONV));
718 >                oo = addarg(oo, sskip2(vval(OCONV), 1));
719 >        } else
720 >                oo = addarg(oo, vval(OCONV));
721   }
722  
723  
# Line 629 | Line 732 | mkillumopts(                           /* get mkillum options */
732                  sprintf(mo, " -n %d", nprocs);
733          else
734                  *mo = '\0';
735 <        if (vdef(MKILLUM))
736 <                if (vval(MKILLUM)[0] != '-') {
737 <                        atos(c_mkillum, sizeof(c_mkillum), vval(MKILLUM));
738 <                        mo = addarg(mo, sskip2(vval(MKILLUM), 1));
739 <                } else
740 <                        mo = addarg(mo, vval(MKILLUM));
735 >        if (!vdef(MKILLUM))
736 >                return;
737 >        if (vval(MKILLUM)[0] != '-') {
738 >                atos(c_mkillum, sizeof(c_mkillum), vval(MKILLUM));
739 >                mo = addarg(mo, sskip2(vval(MKILLUM), 1));
740 >        } else
741 >                mo = addarg(mo, vval(MKILLUM));
742   }
743  
744  
745   static void
746 + mkpmapopts(                             /* get mkpmap options */
747 +        char    *mo
748 + )
749 + {
750 +        /* BEWARE:  This may be called via setdefaults(), so no assumptions */
751 +
752 +        *mo = '\0';
753 +        if (!vdef(MKPMAP))
754 +                return;
755 +        if (vval(MKPMAP)[0] != '-') {
756 +                atos(c_mkpmap, sizeof(c_mkpmap), vval(MKPMAP));
757 +                mo = addarg(mo, sskip2(vval(MKPMAP), 1));
758 +        } else
759 +                mo = addarg(mo, vval(MKPMAP));
760 + }
761 +
762 +
763 + static void
764   checkambfile(void)                      /* check date on ambient file */
765   {
766          time_t  afdate;
# Line 679 | Line 801 | renderopts(                    /* set rendering options */
801          char    *po
802   )
803   {
804 +        char    pmapf[256], *bw;
805 +
806          switch(vscale(QUALITY)) {
807          case LOW:
808                  lowqopts(op, po);
# Line 689 | Line 813 | renderopts(                    /* set rendering options */
813          case HIGH:
814                  hiqopts(op, po);
815                  break;
816 +        }
817 +        if (vdef(PGMAP)) {
818 +                bw = sskip2(vval(PGMAP), 2);
819 +                atos(pmapf, sizeof(pmapf), vval(PGMAP));
820 +                op = addarg(addarg(op, "-ap"), pmapf);
821 +                if (atoi(bw) > 0) op = addarg(op, bw);
822 +        }
823 +        if (vdef(PCMAP)) {
824 +                bw = sskip2(vval(PCMAP), 2);
825 +                atos(pmapf, sizeof(pmapf), vval(PCMAP));
826 +                op = addarg(addarg(op, "-ap"), pmapf);
827 +                if (atoi(bw) > 0) op = addarg(op, bw);
828          }
829          if (vdef(RENDER))
830                  op = addarg(op, vval(RENDER));

Diff Legend

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