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

Comparing ray/src/util/rad.c (file contents):
Revision 2.110 by greg, Wed Dec 4 20:07:07 2013 UTC vs.
Revision 2.128 by greg, Thu Jan 21 21:36:28 2021 UTC

# Line 9 | Line 9 | static const char      RCSid[] = "$Id$";
9  
10   #include <ctype.h>
11   #include <time.h>
12 + #include <signal.h>
13  
14   #include "platform.h"
15   #include "rtprocess.h"
# Line 16 | Line 17 | static const char      RCSid[] = "$Id$";
17   #include "paths.h"
18   #include "vars.h"
19  
20 < #ifdef _WIN32
20 > #if defined(_WIN32) || defined(_WIN64)
21    #define DELCMD "del"
22    #define RENAMECMD "rename"
23   #else
# Line 24 | Line 25 | static const char      RCSid[] = "$Id$";
25    #define RENAMECMD "mv"
26    #include <sys/types.h>
27    #include <sys/wait.h>
28 +  #include <signal.h>
29   #endif
30  
31                                  /* variables (alphabetical by name) */
# Line 35 | Line 37 | static const char      RCSid[] = "$Id$";
37   #define INDIRECT        5               /* indirection in lighting */
38   #define MATERIAL        6               /* material files */
39   #define MKILLUM         7               /* mkillum options */
40 < #define OBJECT          8               /* object files */
41 < #define OCONV           9               /* oconv options */
42 < #define OCTREE          10              /* octree file name */
43 < #define OPTFILE         11              /* rendering options file */
44 < #define PENUMBRAS       12              /* shadow penumbras are desired */
45 < #define PFILT           13              /* pfilt options */
46 < #define PICTURE         14              /* picture file root name */
47 < #define QUALITY         15              /* desired rendering quality */
48 < #define RAWFILE         16              /* raw picture file root name */
49 < #define RENDER          17              /* rendering options */
50 < #define REPORT          18              /* report frequency and errfile */
51 < #define RESOLUTION      19              /* maximum picture resolution */
52 < #define RPICT           20              /* rpict parameters */
53 < #define RVU             21              /* rvu parameters */
54 < #define SCENE           22              /* scene files */
55 < #define UP              23              /* view up (X, Y or Z) */
56 < #define VARIABILITY     24              /* level of light variability */
57 < #define VIEWS           25              /* view(s) for picture(s) */
58 < #define ZFILE           26              /* distance file root name */
59 < #define ZONE            27              /* simulation zone */
40 > #define MKPMAP          8               /* mkpmap options */
41 > #define OBJECT          9               /* object files */
42 > #define OCONV           10              /* oconv options */
43 > #define OCTREE          11              /* octree file name */
44 > #define OPTFILE         12              /* rendering options file */
45 > #define PCMAP           13              /* caustic photon map */
46 > #define PENUMBRAS       14              /* shadow penumbras are desired */
47 > #define PFILT           15              /* pfilt options */
48 > #define PGMAP           16              /* global photon map */
49 > #define PICTURE         17              /* picture file root name */
50 > #define QUALITY         18              /* desired rendering quality */
51 > #define RAWFILE         19              /* raw picture file root name */
52 > #define RENDER          20              /* rendering options */
53 > #define REPORT          21              /* report frequency and errfile */
54 > #define RESOLUTION      22              /* maximum picture resolution */
55 > #define RPICT           23              /* rpict parameters */
56 > #define RVU             24              /* rvu parameters */
57 > #define SCENE           25              /* scene files */
58 > #define UP              26              /* view up (X, Y or Z) */
59 > #define VARIABILITY     27              /* level of light variability */
60 > #define VIEWS           28              /* view(s) for picture(s) */
61 > #define ZFILE           29              /* distance file root name */
62 > #define ZONE            30              /* simulation zone */
63                                  /* total number of variables */
64 < int NVARS = 28;
64 > int NVARS = 31;
65  
66   VARIABLE        vv[] = {                /* variable-value pairs */
67          {"AMBFILE",     3,      0,      NULL,   onevalue},
# Line 67 | Line 72 | VARIABLE       vv[] = {                /* variable-value pairs */
72          {"INDIRECT",    3,      0,      NULL,   intvalue},
73          {"materials",   3,      0,      NULL,   catvalues},
74          {"mkillum",     3,      0,      NULL,   catvalues},
75 +        {"mkpmap",      3,      0,      NULL,   catvalues},
76          {"objects",     3,      0,      NULL,   catvalues},
77          {"oconv",       3,      0,      NULL,   catvalues},
78          {"OCTREE",      3,      0,      NULL,   onevalue},
79          {"OPTFILE",     3,      0,      NULL,   onevalue},
80 +        {"PCMAP",       2,      0,      NULL,   onevalue},
81          {"PENUMBRAS",   3,      0,      NULL,   boolvalue},
82          {"pfilt",       2,      0,      NULL,   catvalues},
83 +        {"PGMAP",       2,      0,      NULL,   onevalue},
84          {"PICTURE",     3,      0,      NULL,   onevalue},
85          {"QUALITY",     3,      0,      NULL,   qualvalue},
86          {"RAWFILE",     3,      0,      NULL,   onevalue},
# Line 107 | Line 115 | time_t oct0date;               /* date of pre-mkillum octree */
115   char    *oct1name;              /* name of post-mkillum octree */
116   time_t  oct1date;               /* date of post-mkillum octree (>= matdate) */
117  
118 + char    *pgmapname;             /* name of global photon map */
119 + time_t  pgmapdate;              /* date of global photon map (>= oct1date) */
120 + char    *pcmapname;             /* name of caustic photon map */
121 + time_t  pcmapdate;              /* date of caustic photon map (>= oct1date) */
122 +
123   int     nowarn = 0;             /* no warnings */
124   int     explicate = 0;          /* explicate variables */
125   int     silent = 0;             /* do work silently */
# Line 121 | Line 134 | char   *viewselect = NULL;     /* specific view only */
134                                  /* command paths */
135   char    c_oconv[256] = "oconv";
136   char    c_mkillum[256] = "mkillum";
137 + char    c_mkpmap[256] = "mkpmap";
138   char    c_rvu[256] = "rvu";
139   char    c_rpict[256] = DEF_RPICT_PATH;
140   char    c_rpiece[] = "rpiece";
# Line 144 | Line 158 | static void checkfiles(void);
158   static void getoctcube(double   org[3], double  *sizp);
159   static void setdefaults(void);
160   static void oconv(void);
161 + static void mkpmap(void);
162   static char * addarg(char       *op, char       *arg);
163   static void oconvopts(char      *oo);
164   static void mkillumopts(char    *mo);
165 + static void mkpmapopts(char     *mo);
166   static void checkambfile(void);
167   static double ambval(void);
168   static void renderopts(char     *op, char       *po);
# Line 243 | Line 259 | main(
259                  printvars(stdout);
260                                  /* build octree (and run mkillum) */
261          oconv();
262 +                                /* run mkpmap if indicated */
263 +        mkpmap();
264                                  /* check date on ambient file */
265          checkambfile();
266                                  /* run simulation */
# Line 335 | Line 353 | newfname(              /* create modified file name */
353   static void
354   checkfiles(void)                        /* check for existence and modified times */
355   {
356 +        char    fntemp[256];
357          time_t  objdate;
358  
359          if (!vdef(OCTREE)) {
# Line 364 | Line 383 | checkfiles(void)                       /* check for existence and modified
383                                  vnam(OCTREE), vnam(SCENE), vnam(ILLUM));
384                  quit(1);
385          }
386 +        if (vdef(PGMAP)) {
387 +                if (!*sskip2(vval(PGMAP),1)) {
388 +                        fprintf(stderr, "%s: '%s' missing # photons argument\n",
389 +                                        progname, vnam(PGMAP));
390 +                        quit(1);
391 +                }
392 +                atos(fntemp, sizeof(fntemp), vval(PGMAP));
393 +                pgmapname = savqstr(fntemp);
394 +                pgmapdate = fdate(pgmapname);
395 +        }
396 +        if (vdef(PCMAP)) {
397 +                if (!*sskip2(vval(PCMAP),1)) {
398 +                        fprintf(stderr, "%s: '%s' missing # photons argument\n",
399 +                                        progname, vnam(PCMAP));
400 +                        quit(1);
401 +                }
402 +                atos(fntemp, sizeof(fntemp), vval(PCMAP));
403 +                pcmapname = savqstr(fntemp);
404 +                pcmapdate = fdate(pcmapname);
405 +        }
406          matdate = checklast(vval(MATERIAL));
407   }      
408  
# Line 517 | Line 556 | oconv(void)                            /* run oconv and mkillum if necessary *
556                  return;
557                                                  /* make octree0 */
558          if ((oct0date < scenedate) | (oct0date < illumdate)) {
559 <                if (touchonly && oct0date)
560 <                        touch(oct0name);
561 <                else {                          /* build command */
559 >                if (touchonly && (oct0date || oct1date)) {
560 >                        if (oct0date)
561 >                                touch(oct0name);
562 >                } else {                        /* build command */
563                          if (octreedate)
564                                  sprintf(combuf, "%s%s -i %s %s > %s", c_oconv,
565                                          ocopts, vval(OCTREE),
# Line 544 | Line 584 | oconv(void)                            /* run oconv and mkillum if necessary *
584                          oct0date = octreedate;
585                  if (oct0date < illumdate)       /* ditto */
586                          oct0date = illumdate;
587 <                }
587 >        }
588          if (touchonly && oct1date)
589                  touch(oct1name);
590          else {
# Line 585 | Line 625 | oconv(void)                            /* run oconv and mkillum if necessary *
625   }
626  
627  
628 + static void
629 + mkpmap(void)                    /* run mkpmap if indicated */
630 + {
631 +        char    combuf[2048], *cp;
632 +        time_t  tnow;
633 +                                /* nothing to do? */
634 +        if ((pgmapname == NULL) | (pgmapdate >= oct1date) &&
635 +                        (pcmapname == NULL) | (pcmapdate >= oct1date))
636 +                return;
637 +                                /* just update existing file dates? */
638 +        if (touchonly && (pgmapname == NULL) | (pgmapdate > 0) &&
639 +                        (pcmapname == NULL) | (pcmapdate > 0)) {
640 +                if (pgmapname != NULL)
641 +                        touch(pgmapname);
642 +                if (pcmapname != NULL)
643 +                        touch(pcmapname);
644 +        } else {                /* else need to (re)run pkpmap */
645 +                strcpy(combuf, c_mkpmap);
646 +                for (cp = combuf; *cp; cp++)
647 +                        ;
648 +                mkpmapopts(cp);
649 +                                /* force file overwrite */
650 +                cp = addarg(cp, "-fo+");
651 +                if (vdef(REPORT)) {
652 +                        char    errfile[256];
653 +                        int     n;
654 +                        double  minutes;
655 +                        n = sscanf(vval(REPORT), "%lf %s", &minutes, errfile);
656 +                        if (n == 2)
657 +                                sprintf(cp, " -t %d -e %s", (int)(minutes*60), errfile);
658 +                        else if (n == 1)
659 +                                sprintf(cp, " -t %d", (int)(minutes*60));
660 +                        else
661 +                                badvalue(REPORT);
662 +                }
663 +                if (pgmapname != NULL && pgmapdate < oct1date) {
664 +                        cp = addarg(cp, "-apg");
665 +                        addarg(cp, vval(PGMAP));
666 +                        cp = sskip(sskip(cp));  /* remove any bandwidth */
667 +                        *cp = '\0';
668 +                }
669 +                if (pcmapname != NULL && pcmapdate < oct1date) {
670 +                        cp = addarg(cp, "-apc");
671 +                        addarg(cp, vval(PCMAP));
672 +                        cp = sskip(sskip(cp));  /* remove any bandwidth */
673 +                        *cp = '\0';
674 +                }
675 +                cp = addarg(cp, oct1name);
676 +                if (runcom(combuf)) {
677 +                        fprintf(stderr, "%s: error running %s\n",
678 +                                        progname, c_mkpmap);
679 +                        if (pgmapname != NULL && pgmapdate < oct1date)
680 +                                unlink(pgmapname);
681 +                        if (pcmapname != NULL && pcmapdate < oct1date)
682 +                                unlink(pcmapname);
683 +                        quit(1);
684 +                }
685 +        }
686 +        tnow = time((time_t *)NULL);
687 +        if (pgmapname != NULL)
688 +                pgmapdate = tnow;
689 +        if (pcmapname != NULL)
690 +                pcmapdate = tnow;
691 +        oct1date = tnow;        /* trigger ambient file removal if needed */
692 + }
693 +
694 +
695   static char *
696   addarg(                         /* append argument and advance pointer */
697   char    *op,
# Line 608 | Line 715 | oconvopts(                             /* get oconv options */
715          /* BEWARE:  This may be called via setdefaults(), so no assumptions */
716  
717          *oo = '\0';
718 <        if (vdef(OCONV))
719 <                if (vval(OCONV)[0] != '-') {
720 <                        atos(c_oconv, sizeof(c_oconv), vval(OCONV));
721 <                        oo = addarg(oo, sskip2(vval(OCONV), 1));
722 <                } else
723 <                        oo = addarg(oo, vval(OCONV));
718 >        if (!vdef(OCONV))
719 >                return;
720 >        if (vval(OCONV)[0] != '-') {
721 >                atos(c_oconv, sizeof(c_oconv), vval(OCONV));
722 >                oo = addarg(oo, sskip2(vval(OCONV), 1));
723 >        } else
724 >                oo = addarg(oo, vval(OCONV));
725   }
726  
727  
# Line 628 | Line 736 | mkillumopts(                           /* get mkillum options */
736                  sprintf(mo, " -n %d", nprocs);
737          else
738                  *mo = '\0';
739 <        if (vdef(MKILLUM))
740 <                if (vval(MKILLUM)[0] != '-') {
741 <                        atos(c_mkillum, sizeof(c_mkillum), vval(MKILLUM));
742 <                        mo = addarg(mo, sskip2(vval(MKILLUM), 1));
743 <                } else
744 <                        mo = addarg(mo, vval(MKILLUM));
739 >        if (!vdef(MKILLUM))
740 >                return;
741 >        if (vval(MKILLUM)[0] != '-') {
742 >                atos(c_mkillum, sizeof(c_mkillum), vval(MKILLUM));
743 >                mo = addarg(mo, sskip2(vval(MKILLUM), 1));
744 >        } else
745 >                mo = addarg(mo, vval(MKILLUM));
746   }
747  
748  
749   static void
750 + mkpmapopts(                             /* get mkpmap options */
751 +        char    *mo
752 + )
753 + {
754 +        /* BEWARE:  This may be called via setdefaults(), so no assumptions */
755 +
756 +        if (nprocs > 1)
757 +                sprintf(mo, " -n %d", nprocs);
758 +        else
759 +                *mo = '\0';
760 +        if (!vdef(MKPMAP))
761 +                return;
762 +        if (vval(MKPMAP)[0] != '-') {
763 +                atos(c_mkpmap, sizeof(c_mkpmap), vval(MKPMAP));
764 +                mo = addarg(mo, sskip2(vval(MKPMAP), 1));
765 +        } else
766 +                mo = addarg(mo, vval(MKPMAP));
767 + }
768 +
769 +
770 + static void
771   checkambfile(void)                      /* check date on ambient file */
772   {
773          time_t  afdate;
# Line 660 | Line 790 | ambval(void)                           /* compute ambient value */
790   {
791          if (vdef(EXPOSURE)) {
792                  if (vval(EXPOSURE)[0] == '+' || vval(EXPOSURE)[0] == '-')
793 <                        return(.5/pow(2.,vflt(EXPOSURE)));
794 <                return(.5/vflt(EXPOSURE));
793 >                        return(.18/pow(2.,vflt(EXPOSURE)));
794 >                return(.18/vflt(EXPOSURE));
795          }
796          if (vlet(ZONE) == 'E')
797                  return(10.);
# Line 678 | Line 808 | renderopts(                    /* set rendering options */
808          char    *po
809   )
810   {
811 +        char    pmapf[256], *bw;
812 +
813 +        if (vdef(PGMAP)) {
814 +                *op = '\0';
815 +                bw = sskip2(vval(PGMAP), 2);
816 +                atos(pmapf, sizeof(pmapf), vval(PGMAP));
817 +                op = addarg(addarg(op, "-ap"), pmapf);
818 +                if (atoi(bw) > 0) op = addarg(op, bw);
819 +        }
820          switch(vscale(QUALITY)) {
821          case LOW:
822                  lowqopts(op, po);
# Line 689 | Line 828 | renderopts(                    /* set rendering options */
828                  hiqopts(op, po);
829                  break;
830          }
831 <        if (vdef(RENDER))
831 >        if (vdef(PCMAP)) {
832 >                bw = sskip2(vval(PCMAP), 2);
833 >                atos(pmapf, sizeof(pmapf), vval(PCMAP));
834 >                op = addarg(addarg(op, "-ap"), pmapf);
835 >                if (atoi(bw) > 0) op = addarg(op, bw);
836 >        }
837 >        if (vdef(RENDER)) {
838                  op = addarg(op, vval(RENDER));
839 +                bw = strstr(vval(RENDER), "-aa ");
840 +                if (bw != NULL && atof(bw+4) <= FTINY)
841 +                        overture = 0;
842 +        }
843          if (rvdevice != NULL) {
844 <                if (vdef(RVU))
844 >                if (vdef(RVU)) {
845                          if (vval(RVU)[0] != '-') {
846                                  atos(c_rvu, sizeof(c_rvu), vval(RVU));
847                                  po = addarg(po, sskip2(vval(RVU), 1));
848                          } else
849                                  po = addarg(po, vval(RVU));
850 +                }
851          } else {
852 <                if (vdef(RPICT))
852 >                if (vdef(RPICT)) {
853                          if (vval(RPICT)[0] != '-') {
854                                  atos(c_rpict, sizeof(c_rpict), vval(RPICT));
855                                  po = addarg(po, sskip2(vval(RPICT), 1));
856                          } else
857                                  po = addarg(po, vval(RPICT));
858 +                }
859          }
860   }
861  
# Line 931 | Line 1082 | hiqopts(                               /* high quality rendering options */
1082   }
1083  
1084  
1085 < #ifdef _WIN32
1085 > #if defined(_WIN32) || defined(_WIN64)
1086   static void
1087   setenv(                 /* set an environment variable */
1088          char    *vname,
# Line 978 | Line 1129 | xferopts(                              /* transfer options if indicated */
1129                          syserr(vval(OPTFILE));
1130                  sprintf(ro, " @%s", vval(OPTFILE));
1131          }
1132 < #ifdef _WIN32
1132 > #if defined(_WIN32) || defined(_WIN64)
1133          else if (n > 50) {
1134                  setenv("ROPT", ro+1);
1135                  strcpy(ro, " $ROPT");
# Line 1005 | Line 1156 | pfiltopts(                             /* get pfilt options */
1156                  po = addarg(po, "-m .25");
1157                  break;
1158          }
1159 <        if (vdef(PFILT))
1159 >        if (vdef(PFILT)) {
1160                  if (vval(PFILT)[0] != '-') {
1161                          atos(c_pfilt, sizeof(c_pfilt), vval(PFILT));
1162                          po = addarg(po, sskip2(vval(PFILT), 1));
1163                  } else
1164                          po = addarg(po, vval(PFILT));
1165 +        }
1166   }
1167  
1168  
# Line 1155 | Line 1307 | specview(                              /* get proper view spec from vs */
1307          if (cp == viewopts)             /* append any additional options */
1308                  vs++;           /* skip prefixed space if unneeded */
1309          strcpy(cp, vs);
1310 < #ifdef _WIN32
1310 > #if defined(_WIN32) || defined(_WIN64)
1311          if (strlen(viewopts) > 40) {
1312                  setenv("VIEW", viewopts);
1313                  return("$VIEW");
# Line 1201 | Line 1353 | getview(                               /* get view n, or NULL if none */
1353          }
1354   numview:
1355          mv = nvalue(VIEWS, n);          /* use view n */
1356 <        if ((vn != NULL) & (mv != NULL))
1356 >        if ((vn != NULL) & (mv != NULL)) {
1357                  if (*mv != '-') {
1358                          char    *mv2 = mv;
1359                          while (*mv2 && !isspace(*mv2))
# Line 1209 | Line 1361 | numview:
1361                          *vn = '\0';
1362                  } else
1363                          sprintf(vn, "%d", n+1);
1364 <
1364 >        }
1365          return(specview(mv));
1366   }
1367  
# Line 1223 | Line 1375 | myprintview(                   /* print out selected view */
1375          VIEW    vwr;
1376          char    buf[128];
1377          char    *cp;
1378 < #ifdef _WIN32
1378 > #if defined(_WIN32) || defined(_WIN64)
1379   /* XXX Should we allow something like this for all platforms? */
1380   /* XXX Or is it still required at all? */
1381   again:
1382   #endif
1383          if (vopts == NULL)
1384                  return(-1);
1385 < #ifdef _WIN32
1385 > #if defined(_WIN32) || defined(_WIN64)
1386          if (vopts[0] == '$') {
1387                  vopts = getenv(vopts+1);
1388                  goto again;
# Line 1377 | Line 1529 | rpict(                         /* run rpict and pfilt for each view */
1529                                  fprintf(stderr, "%s: cannot create\n", pfile);
1530                                  quit(1);
1531                          }
1532 +                        pfile[-5] = '\0';
1533                          pfile = NULL;
1534                  }
1535          }

Diff Legend

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