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.91 by greg, Thu Dec 18 23:26:57 2008 UTC vs.
Revision 2.108 by greg, Mon Aug 12 21:01:00 2013 UTC

# Line 116 | Line 116 | int    sayview = 0;            /* print view out */
116   char    *rvdevice = NULL;       /* rvu output device */
117   char    *viewselect = NULL;     /* specific view only */
118  
119 + #define DEF_RPICT_PATH  "rpict"         /* default rpict path */
120 +
121                                  /* command paths */
122   char    c_oconv[256] = "oconv";
123   char    c_mkillum[256] = "mkillum";
124   char    c_rvu[256] = "rvu";
125 < char    c_rpict[256] = "rpict";
125 > char    c_rpict[256] = DEF_RPICT_PATH;
126 > char    c_rpiece[] = "rpiece";
127   char    c_pfilt[256] = "pfilt";
128  
129   int     overture = 0;           /* overture calculation needed */
# Line 162 | Line 165 | static int touch(char  *fn);
165   static int runcom(char  *cs);
166   static int rmfile(char  *fn);
167   static int mvfile(char  *fold, char     *fnew);
168 < static int next_process(void);
168 > static int next_process(int     reserve);
169   static void wait_process(int    all);
170   static void finish_process(void);
171   static void badvalue(int        vc);
# Line 218 | Line 221 | main(
221          if (i >= argc)
222                  goto userr;
223          rifname = argv[i];
221 –                                /* check command-line options */
222 –        if ((nprocs > 1) & (viewselect != NULL) & (rvdevice == NULL))
223 –                nprocs = 1;
224                                  /* assign Radiance root file name */
225          rootname(radname, rifname);
226                                  /* load variable values */
# Line 558 | Line 558 | oconv(void)                            /* run oconv and mkillum if necessary *
558                          unlink(illumtmp);
559                          quit(1);
560                  }
561 +                rmfile(oct0name);
562                                                  /* make octree1 (frozen) */
563                  if (octreedate)
564                          sprintf(combuf, "%s%s -f -i %s %s > %s", c_oconv,
# Line 751 | Line 752 | lowqopts(                      /* low quality rendering options */
752                  op = addarg(op, "-ds .4");
753          else
754                  op = addarg(op, "-ds 0");
755 <        op = addarg(op, "-dt .2 -dc .25 -dr 0 -sj 0 -st .5");
755 >        op = addarg(op, "-dt .2 -dc .25 -dr 0 -ss 0 -st .5");
756          if (vdef(AMBFILE)) {
757                  sprintf(op, " -af %s", vval(AMBFILE));
758                  op += strlen(op);
# Line 772 | Line 773 | lowqopts(                      /* low quality rendering options */
773          d = ambval();
774          sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
775          op += strlen(op);
776 <        op = addarg(op, "-lr 6 -lw .01");
776 >        op = addarg(op, "-lr 6 -lw .003");
777   }
778  
779  
# Line 826 | Line 827 | medqopts(                      /* medium quality rendering options */
827                  op = addarg(op, "-ds .2 -dj .9");
828          else
829                  op = addarg(op, "-ds .3");
830 <        op = addarg(op, "-dt .1 -dc .5 -dr 1 -sj .7 -st .1");
830 >        op = addarg(op, "-dt .1 -dc .5 -dr 1 -ss 1 -st .1");
831          if ( (overture = vint(INDIRECT)) ) {
832                  sprintf(op, " -ab %d", overture);
833                  op += strlen(op);
# Line 850 | Line 851 | medqopts(                      /* medium quality rendering options */
851          d = ambval();
852          sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
853          op += strlen(op);
854 <        op = addarg(op, "-lr 8 -lw .002");
854 >        op = addarg(op, "-lr 8 -lw 1e-4");
855   }
856  
857  
# Line 904 | Line 905 | hiqopts(                               /* high quality rendering options */
905                  op = addarg(op, "-ds .1 -dj .9");
906          else
907                  op = addarg(op, "-ds .2");
908 <        op = addarg(op, "-dt .05 -dc .75 -dr 3 -sj 1 -st .01");
908 >        op = addarg(op, "-dt .05 -dc .75 -dr 3 -ss 16 -st .01");
909          sprintf(op, " -ab %d", overture=vint(INDIRECT)+1);
910          op += strlen(op);
911          if (vdef(AMBFILE)) {
# Line 926 | Line 927 | hiqopts(                               /* high quality rendering options */
927          d = ambval();
928          sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
929          op += strlen(op);
930 <        op = addarg(op, "-lr 12 -lw .0005");
930 >        op = addarg(op, "-lr 12 -lw 1e-5");
931   }
932  
933  
# Line 1123 | Line 1124 | specview(                              /* get proper view spec from vs */
1124                  while (!isspace(*vs))           /* else skip id */
1125                          if (!*vs++)
1126                                  return(NULL);
1127 <                if (upax) {                     /* specify up vector */
1127 >                if (upax) {                     /* prepend up vector */
1128                          strcpy(cp, vup[upax+3]);
1129                          cp += strlen(cp);
1130                  }
# Line 1152 | Line 1153 | getview(                               /* get view n, or NULL if none */
1153          if (viewselect != NULL) {               /* command-line selected */
1154                  if (n)                          /* only do one */
1155                          return(NULL);
1156 +                                        
1157 +                if (isint(viewselect)) {        /* view number? */
1158 +                        n = atoi(viewselect)-1;
1159 +                        goto numview;
1160 +                }
1161                  if (viewselect[0] == '-') {     /* already specified */
1162 <                        if (vn != NULL) *vn = '\0';
1162 >                        if (vn != NULL)
1163 >                                strcpy(vn, "0");
1164                          return(viewselect);
1165                  }
1166                  if (vn != NULL) {
# Line 1162 | Line 1169 | getview(                               /* get view n, or NULL if none */
1169                                  ;
1170                          *vn = '\0';
1171                  }
1165 –                                                /* view number? */
1166 –                if (isint(viewselect))
1167 –                        return(specview(nvalue(VIEWS, atoi(viewselect)-1)));
1172                                                  /* check list */
1173                  while ((mv = nvalue(VIEWS, n++)) != NULL)
1174                          if (matchword(viewselect, mv))
1175                                  return(specview(mv));
1176 +
1177                  return(specview(viewselect));   /* standard view? */
1178          }
1179 + numview:
1180          mv = nvalue(VIEWS, n);          /* use view n */
1181 <        if ((vn != NULL) & (mv != NULL)) {
1182 <                register char   *mv2 = mv;
1183 <                if (*mv2 != '-')
1181 >        if ((vn != NULL) & (mv != NULL))
1182 >                if (*mv != '-') {
1183 >                        register char   *mv2 = mv;
1184                          while (*mv2 && !isspace(*mv2))
1185                                  *vn++ = *mv2++;
1186 <                *vn = '\0';
1187 <        }
1186 >                        *vn = '\0';
1187 >                } else
1188 >                        sprintf(vn, "%d", n+1);
1189 >
1190          return(specview(mv));
1191   }
1192  
# Line 1247 | Line 1255 | rvu(                           /* run rvu with first view */
1255   }
1256  
1257  
1258 + static int
1259 + syncf_done(                     /* check if an rpiece sync file is complete */
1260 +        char *sfname
1261 + )
1262 + {
1263 +        FILE    *fp = fopen(sfname, "r");
1264 +        int     todo = 1;
1265 +        int     x, y;
1266 +
1267 +        if (fp == NULL)
1268 +                return(0);
1269 +        if (fscanf(fp, "%d %d", &x, &y) != 2)
1270 +                goto checked;
1271 +        todo = x*y;             /* total number of tiles */
1272 +        if (fscanf(fp, "%d %d", &x, &y) != 2 || (x != 0) | (y != 0))
1273 +                goto checked;
1274 +                                /* XXX assume no redundant tiles */
1275 +        while (fscanf(fp, "%d %d", &x, &y) == 2)
1276 +                if (!--todo)
1277 +                        break;
1278 + checked:
1279 +        fclose(fp);
1280 +        return(!todo);
1281 + }
1282 +
1283 +
1284   static void
1285   rpict(                          /* run rpict and pfilt for each view */
1286          char    *opts,
1287          char    *po
1288   )
1289   {
1290 <        char    combuf[PATH_MAX];
1290 > #define do_rpiece       (sfile[0]!='\0')
1291 >        char    combuf[5*PATH_MAX+512];
1292          char    rawfile[PATH_MAX], picfile[PATH_MAX];
1293          char    zopt[PATH_MAX+4], rep[PATH_MAX+16], res[32];
1294 <        char    rppopt[128], *pfile = NULL;
1294 >        char    rppopt[32], sfile[PATH_MAX], *pfile = NULL;
1295          char    pfopts[128];
1296          char    vs[32], *vw;
1297          int     vn, mult;
1298          FILE    *fp;
1299          time_t  rfdt, pfdt;
1300 +        int     xres, yres;
1301 +        double  aspect;
1302 +        int     n;
1303                                          /* get pfilt options */
1304          pfiltopts(pfopts);
1305                                          /* get resolution, reporting */
# Line 1276 | Line 1314 | rpict(                         /* run rpict and pfilt for each view */
1314                  mult = 3;
1315                  break;
1316          }
1317 <        {
1318 <                int     xres, yres;
1319 <                double  aspect;
1320 <                int     n;
1321 <                n = sscanf(vval(RESOLUTION), "%d %d %lf", &xres, &yres, &aspect);
1322 <                if (n == 3)
1323 <                        sprintf(res, "-x %d -y %d -pa %.3f",
1324 <                                        mult*xres, mult*yres, aspect);
1325 <                else if (n) {
1326 <                        if (n == 1) yres = xres;
1289 <                        sprintf(res, "-x %d -y %d", mult*xres, mult*yres);
1290 <                } else
1291 <                        badvalue(RESOLUTION);
1292 <        }
1317 >        n = sscanf(vval(RESOLUTION), "%d %d %lf", &xres, &yres, &aspect);
1318 >        if (n == 3)
1319 >                sprintf(res, "-x %d -y %d -pa %.3f",
1320 >                                mult*xres, mult*yres, aspect);
1321 >        else if (n) {
1322 >                aspect = 1.;
1323 >                if (n == 1) yres = xres;
1324 >                sprintf(res, "-x %d -y %d", mult*xres, mult*yres);
1325 >        } else
1326 >                badvalue(RESOLUTION);
1327          rep[0] = '\0';
1328          if (vdef(REPORT)) {
1329                  double  minutes;
1296 –                int     n;
1330                  n = sscanf(vval(REPORT), "%lf %s", &minutes, rawfile);
1331                  if (n == 2)
1332                          sprintf(rep, " -t %d -e %s", (int)(minutes*60), rawfile);
# Line 1303 | Line 1336 | rpict(                         /* run rpict and pfilt for each view */
1336                          badvalue(REPORT);
1337          }
1338                                          /* set up parallel rendering */
1339 <        if ((nprocs > 1) & (!vdef(ZFILE))) {
1340 <                strcpy(rppopt, "-S 1 -PP pfXXXXXX");
1341 <                pfile = rppopt+9;
1342 <                if (mktemp(pfile) == NULL)
1339 >        sfile[0] = '\0';
1340 >        if ((nprocs > 1) & !touchonly & !vdef(ZFILE) &&
1341 >                                        getview(0, vs) != NULL) {
1342 >                if (!strcmp(c_rpict, DEF_RPICT_PATH) &&
1343 >                                getview(1, NULL) == NULL) {
1344 >                        sprintf(sfile, "%s_%s_rpsync.txt",
1345 >                                vdef(RAWFILE) ? vval(RAWFILE) : vval(PICTURE),
1346 >                                        vs);
1347 >                        strcpy(rppopt, "-PP pfXXXXXX");
1348 >                } else {
1349 >                        strcpy(rppopt, "-S 1 -PP pfXXXXXX");
1350 >                }
1351 >                pfile = rppopt + strlen(rppopt) - 8;
1352 >                if (mktemp(pfile) == NULL) {
1353 >                        if (do_rpiece) {
1354 >                                fprintf(stderr, "%s: cannot create\n", pfile);
1355 >                                quit(1);
1356 >                        }
1357                          pfile = NULL;
1358 +                }
1359          }
1360          vn = 0;                                 /* do each view */
1361          while ((vw = getview(vn++, vs)) != NULL) {
1362                  if (sayview)
1363                          myprintview(vw, stdout);
1316 –                if (!vs[0])
1317 –                        sprintf(vs, "%d", vn);
1364                  sprintf(picfile, "%s_%s.hdr", vval(PICTURE), vs);
1365                  if (vdef(ZFILE))
1366                          sprintf(zopt, " -z %s_%s.zbf", vval(ZFILE), vs);
# Line 1336 | Line 1382 | rpict(                         /* run rpict and pfilt for each view */
1382                                  touch(picfile);
1383                          continue;
1384                  }
1385 <                if (next_process()) {           /* parallel running? */
1385 >                                                /* parallel running? */
1386 >                if (do_rpiece) {
1387 >                        if (rfdt < oct1date || !fdate(sfile)) {
1388 >                                int     xdiv = 8+nprocs/3, ydiv = 8+nprocs/3;
1389 >                                if (rfdt >= oct1date) {
1390 >                                        fprintf(stderr,
1391 >                "%s: partial output not created with %s\n", rawfile, c_rpiece);
1392 >                                        quit(1);
1393 >                                }
1394 >                                if (rfdt) {     /* start fresh */
1395 >                                        rmfile(rawfile);
1396 >                                        rfdt = 0;
1397 >                                }
1398 >                                if (!silent)
1399 >                                        printf("\techo %d %d > %s\n",
1400 >                                                        xdiv, ydiv, sfile);
1401 >                                if ((fp = fopen(sfile, "w")) == NULL) {
1402 >                                        fprintf(stderr, "%s: cannot create\n",
1403 >                                                        sfile);
1404 >                                        quit(1);
1405 >                                }
1406 >                                fprintf(fp, "%d %d\n", xdiv, ydiv);
1407 >                                fclose(fp);
1408 >                        }
1409 >                } else if (next_process(0)) {
1410                          if (pfile != NULL)
1411 <                                sleep(20);
1411 >                                sleep(10);
1412                          continue;
1413 <                }
1413 >                } else if (!inchild())
1414 >                        pfile = NULL;
1415                  /* XXX Remember to call finish_process() */
1416                                                  /* build rpict command */
1417 <                if (rfdt >= oct1date) {         /* recover */
1418 <                        sprintf(combuf, "%s%s%s%s%s -ro %s %s", c_rpict,
1419 <                                rep, opts, po, zopt, rawfile, oct1name);
1420 <                        if (runcom(combuf))     /* run rpict */
1417 >                if (rfdt >= oct1date) {         /* already in progress */
1418 >                        if (do_rpiece) {
1419 >                                sprintf(combuf, "%s -R %s %s%s %s %s%s%s -o %s %s",
1420 >                                                c_rpiece, sfile, rppopt, rep, vw,
1421 >                                                res, opts, po, rawfile, oct1name);
1422 >                                while (next_process(1)) {
1423 >                                        sleep(10);
1424 >                                        combuf[strlen(c_rpiece)+2] = 'F';
1425 >                                }
1426 >                        } else
1427 >                                sprintf(combuf, "%s%s%s%s%s -ro %s %s", c_rpict,
1428 >                                        rep, opts, po, zopt, rawfile, oct1name);
1429 >                        if (runcom(combuf))     /* run rpict/rpiece */
1430                                  goto rperror;
1431                  } else {
1432                          if (overture) {         /* run overture calculation */
1433                                  sprintf(combuf,
1434 <                                "%s%s %s%s -x 64 -y 64 -ps 1 %s > %s",
1434 >                                        "%s%s %s%s -x 64 -y 64 -ps 1 %s > %s",
1435                                                  c_rpict, rep, vw, opts,
1436                                                  oct1name, overfile);
1437 <                                if (runcom(combuf)) {
1438 <                                        fprintf(stderr,
1437 >                                if (!do_rpiece || !next_process(0)) {
1438 >                                        if (runcom(combuf)) {
1439 >                                                fprintf(stderr,
1440                                          "%s: error in overture for view %s\n",
1441 <                                                progname, vs);
1442 <                                        quit(1);
1443 <                                }
1441 >                                                        progname, vs);
1442 >                                                quit(1);
1443 >                                        }
1444   #ifndef NULL_DEVICE
1445 <                                rmfile(overfile);
1445 >                                        rmfile(overfile);
1446   #endif
1447 +                                } else if (do_rpiece)
1448 +                                        sleep(20);
1449                          }
1450 <                        sprintf(combuf, "%s%s %s %s%s%s%s %s > %s",
1451 <                                        c_rpict, rep, vw, res, opts, po,
1452 <                                        zopt, oct1name, rawfile);
1453 <                        if (pfile != NULL && inchild()) {
1454 <                                                /* rpict persistent mode */
1450 >                        if (do_rpiece) {
1451 >                                sprintf(combuf, "%s -F %s %s%s %s %s%s%s -o %s %s",
1452 >                                                c_rpiece, sfile, rppopt, rep, vw,
1453 >                                                res, opts, po, rawfile, oct1name);
1454 >                                while (next_process(1))
1455 >                                        sleep(10);
1456 >                        } else {
1457 >                                sprintf(combuf, "%s%s %s %s%s%s%s %s > %s",
1458 >                                                c_rpict, rep, vw, res, opts, po,
1459 >                                                zopt, oct1name, rawfile);
1460 >                        }
1461 >                        if ((pfile != NULL) & !do_rpiece) {
1462                                  if (!silent)
1463                                          printf("\t%s\n", combuf);
1464                                  fflush(stdout);
1465                                  sprintf(combuf, "%s%s %s %s%s%s %s > %s",
1466 <                                                c_rpict, rep, rppopt, res, opts,
1467 <                                                po, oct1name, rawfile);
1466 >                                                c_rpict, rep, rppopt, res,
1467 >                                                opts, po, oct1name, rawfile);
1468                                  fp = popen(combuf, "w");
1469                                  if (fp == NULL)
1470                                          goto rperror;
1471                                  myprintview(vw, fp);
1472                                  if (pclose(fp))
1473                                          goto rperror;
1474 <                        } else {                /* rpict normal mode */
1475 <                                if (runcom(combuf))
1476 <                                        goto rperror;
1474 >                        } else if (runcom(combuf))
1475 >                                goto rperror;
1476 >                }
1477 >                if (do_rpiece) {                /* need to finish raw, first */
1478 >                        finish_process();
1479 >                        wait_process(1);
1480 >                        if (!syncf_done(sfile)) {
1481 >                                fprintf(stderr,
1482 >                        "%s: %s did not complete rendering of view %s\n",
1483 >                                                progname, c_rpiece, vs);
1484 >                                quit(1);
1485                          }
1486                  }
1487                  if (!vdef(RAWFILE) || strcmp(vval(RAWFILE),vval(PICTURE))) {
1488                                                  /* build pfilt command */
1489 <                        if (mult > 1)
1489 >                        if (do_rpiece)
1490 >                                sprintf(combuf,
1491 >                                        "%s%s -x %d -y %d -p %.3f %s > %s",
1492 >                                        c_pfilt, pfopts, xres, yres, aspect,
1493 >                                        rawfile, picfile);
1494 >                        else if (mult > 1)
1495                                  sprintf(combuf, "%s%s -x /%d -y /%d %s > %s",
1496                                          c_pfilt, pfopts, mult, mult,
1497                                          rawfile, picfile);
# Line 1409 | Line 1512 | rpict(                         /* run rpict and pfilt for each view */
1512                          mvfile(rawfile, combuf);
1513                  } else
1514                          rmfile(rawfile);
1515 <                finish_process();               /* exit if child */
1515 >                if (do_rpiece)                  /* done with sync file */
1516 >                        rmfile(sfile);
1517 >                else
1518 >                        finish_process();       /* exit if child */
1519          }
1520          wait_process(1);                /* wait for children to finish */
1521 <        if (pfile != NULL) {            /* clean up rpict persistent mode */
1521 >        if (pfile != NULL) {            /* clean up persistent rpict */
1522                  RT_PID  pid;
1523                  fp = fopen(pfile, "r");
1524                  if (fp != NULL) {
# Line 1426 | Line 1532 | rpict(                         /* run rpict and pfilt for each view */
1532   rperror:
1533          fprintf(stderr, "%s: error rendering view %s\n", progname, vs);
1534          quit(1);
1535 + #undef do_rpiece
1536   }
1537  
1538  
# Line 1438 | Line 1545 | touch(                 /* update a file */
1545                  printf("\ttouch %s\n", fn);
1546          if (!nprocs)
1547                  return(0);
1441 – #ifdef notused
1442 –        if (access(fn, F_OK) == -1)             /* create it */
1443 –                if (close(open(fn, O_WRONLY|O_CREAT, 0666)) == -1)
1444 –                        return(-1);
1445 – #endif
1548          return(setfdate(fn, time((time_t *)NULL)));
1549   }
1550  
# Line 1490 | Line 1592 | mvfile(                /* move a file */
1592  
1593   #ifdef RHAS_FORK_EXEC
1594   static int
1595 < next_process(void)                      /* fork the next process (max. nprocs) */
1595 > next_process(int reserve)               /* fork the next process */
1596   {
1597          RT_PID  child_pid;
1598  
# Line 1501 | Line 1603 | next_process(void)                     /* fork the next process (max. np
1603                                  progname);
1604                  quit(1);
1605          }
1606 +        if (reserve > 0 && children_running >= nprocs-reserve)
1607 +                return(0);              /* caller holding back process(es) */
1608          if (children_running >= nprocs)
1609                  wait_process(0);        /* wait for someone to finish */
1610          fflush(NULL);                   /* flush output */
1611          child_pid = fork();             /* split process */
1612          if (child_pid == 0) {           /* we're the child */
1613                  children_running = -1;
1614 +                nprocs = 1;
1615                  return(0);
1616          }
1617          if (child_pid > 0) {            /* we're the parent */
# Line 1548 | Line 1653 | wait_process(                  /* wait for process(es) to finish */
1653   }
1654   #else   /* ! RHAS_FORK_EXEC */
1655   static int
1656 < next_process(void)
1656 > next_process(int reserve)
1657   {
1658          return(0);                      /* cannot start new process */
1659   }

Diff Legend

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