src/util/rad.c

/* Copyright (c) 1993 Regents of the University of California */

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/*
 * Executive program for oconv, rpict and pfilt
 */

#include "standard.h"
#include "paths.h"
#include <ctype.h>


typedef struct {
        char    *name;          /* variable name */
        short   nick;           /* # characters required for nickname */
        short   nass;           /* # assignments made */
        char    *value;         /* assigned value(s) */
        int     (*fixval)();    /* assignment checking function */
} VARIABLE;

int     onevalue(), catvalues();

                                /* variables */
#define OBJECT          0               /* object files */
#define SCENE           1               /* scene files */
#define MATERIAL        2               /* material files */
#define RENDER          3               /* rendering options */
#define OCONV           4               /* oconv options */
#define PFILT           5               /* pfilt options */
#define VIEW            6               /* view(s) for picture(s) */
#define ZONE            7               /* simulation zone */
#define QUALITY         8               /* desired rendering quality */
#define OCTREE          9               /* octree file name */
#define PICTURE         10              /* picture file name */
#define AMBFILE         11              /* ambient file name */
#define OPTFILE         12              /* rendering options file */
#define EXPOSURE        13              /* picture exposure setting */
#define RESOLUTION      14              /* maximum picture resolution */
#define UP              15              /* view up (X, Y or Z) */
#define INDIRECT        16              /* indirection in lighting */
#define DETAIL          17              /* level of scene detail */
#define PENUMBRAS       18              /* shadow penumbras are desired */
#define VARIABILITY     19              /* level of light variability */
#define REPORT          20              /* report frequency and errfile */
                                /* total number of variables */
#define NVARS           21

VARIABLE        vv[NVARS] = {           /* variable-value pairs */
        {"objects",     3,      0,      NULL,   catvalues},
        {"scene",       3,      0,      NULL,   catvalues},
        {"materials",   3,      0,      NULL,   catvalues},
        {"render",      3,      0,      NULL,   catvalues},
        {"oconv",       3,      0,      NULL,   catvalues},
        {"pfilt",       2,      0,      NULL,   catvalues},
        {"view",        2,      0,      NULL,   NULL},
        {"ZONE",        2,      0,      NULL,   onevalue},
        {"QUALITY",     3,      0,      NULL,   onevalue},
        {"OCTREE",      3,      0,      NULL,   onevalue},
        {"PICTURE",     3,      0,      NULL,   onevalue},
        {"AMBFILE",     3,      0,      NULL,   onevalue},
        {"OPTFILE",     3,      0,      NULL,   onevalue},
        {"EXPOSURE",    3,      0,      NULL,   onevalue},
        {"RESOLUTION",  3,      0,      NULL,   onevalue},
        {"UP",          2,      0,      NULL,   onevalue},
        {"INDIRECT",    3,      0,      NULL,   onevalue},
        {"DETAIL",      3,      0,      NULL,   onevalue},
        {"PENUMBRAS",   3,      0,      NULL,   onevalue},
        {"VARIABILITY", 3,      0,      NULL,   onevalue},
        {"REPORT",      3,      0,      NULL,   onevalue},
};

VARIABLE        *matchvar();
char    *nvalue();
int     vscale();

#define UPPER(c)        ((c)&~0x20)     /* ASCII trick */

#define vnam(vc)        (vv[vc].name)
#define vdef(vc)        (vv[vc].nass)
#define vval(vc)        (vv[vc].value)
#define vint(vc)        atoi(vval(vc))
#define vlet(vc)        UPPER(vval(vc)[0])
#define vbool(vc)       (vlet(vc)=='T')

#define HIGH            2
#define MEDIUM          1
#define LOW             0

int     lowqopts(), medqopts(), hiqopts();
int     (*setqopts[3])() = {lowqopts, medqopts, hiqopts};

#define renderopts      (*setqopts[vscale(QUALITY)])

                                /* overture calculation file */
#ifdef NIX
char    overfile[] = "overture.raw";
#else
char    overfile[] = "/dev/null";
#endif

extern long     fdate(), time();

long    scenedate;              /* date of latest scene or object file */
long    octreedate;             /* date of octree */

int     explicate = 0;          /* explicate variables */
int     silent = 0;             /* do work silently */
int     noaction = 0;           /* don't do anything */
char    *rvdevice = NULL;       /* rview output device */
char    *viewselect = NULL;     /* specific view only */

int     overture = 0;           /* overture calculation needed */

char    *progname;              /* global argv[0] */

char    radname[MAXPATH];       /* root Radiance file name */


main(argc, argv)
int     argc;
char    *argv[];
{
        char    ropts[512];
        int     i;

        progname = argv[0];
                                /* get options */
        for (i = 1; i < argc && argv[i][0] == '-'; i++)
                switch (argv[i][1]) {
                case 's':
                        silent++;
                        break;
                case 'n':
                        noaction++;
                        break;
                case 'e':
                        explicate++;
                        break;
                case 'o':
                        rvdevice = argv[++i];
                        break;
                case 'v':
                        viewselect = argv[++i];
                        break;
                default:
                        goto userr;
                }
        if (i >= argc)
                goto userr;
                                /* assign Radiance root file name */
        rootname(radname, argv[i]);
                                /* load variable values */
        load(argv[i]);
                                /* get any additional assignments */
        for (i++; i < argc; i++)
                setvariable(argv[i]);
                                /* check assignments */
        checkvalues();
                                /* check files and dates */
        checkfiles();
                                /* set default values as necessary */
        setdefaults();
                                /* print all values if requested */
        if (explicate)
                printvals();
                                /* run simulation */
        oconv();
        renderopts(ropts);
        xferopts(ropts);
        if (rvdevice != NULL)
                rview(ropts);
        else
                rpict(ropts);
        exit(0);
userr:
        fprintf(stderr,
        "Usage: %s [-s][-n][-e][-v view][-o dev] rfile [VAR=value ..]\n",
                        progname);
        exit(1);
}


rootname(rn, fn)                /* remove tail from end of fn */
register char   *rn, *fn;
{
        char    *tp, *dp;

        for (tp = NULL, dp = rn; *rn = *fn++; rn++)
                if (ISDIRSEP(*rn))
                        dp = rn;
                else if (*rn == '.')
                        tp = rn;
        if (tp != NULL && tp > dp)
                *tp = '\0';
}


load(rfname)                    /* load Radiance simulation file */
char    *rfname;
{
        FILE    *fp;
        char    buf[512];
        register char   *cp;

        if (rfname == NULL)
                fp = stdin;
        else if ((fp = fopen(rfname, "r")) == NULL)
                syserr(rfname);
        while (fgetline(buf, sizeof(buf), fp) != NULL) {
                for (cp = buf; *cp; cp++) {
                        switch (*cp) {
                        case '\\':
                        case '\n':
                                *cp = ' ';
                                continue;
                        case '#':
                                *cp = '\0';
                                break;
                        default:
                                continue;
                        }
                        break;
                }
                setvariable(buf);
        }
        fclose(fp);
}


setvariable(ass)                /* assign variable according to string */
register char   *ass;
{
        char    varname[32];
        register char   *cp;
        register VARIABLE       *vp;
        register int    i;
        int     n;

        while (isspace(*ass))           /* skip leading space */
                ass++;
        cp = varname;                   /* extract name */
        while (cp < varname+sizeof(varname)-1
                        && *ass && !isspace(*ass) && *ass != '=')
                *cp++ = *ass++;
        *cp = '\0';
        if (!varname[0])
                return;         /* no variable name! */
                                        /* trim value */
        while (isspace(*ass) || *ass == '=')
                ass++;
        cp = ass + strlen(ass);
        do
                *cp-- = '\0';
        while (cp >= ass && isspace(*cp));
        n = cp - ass + 1;
        if (!n) {
                fprintf(stderr, "%s: warning - missing value for variable '%s'\n",
                                progname, varname);
                return;
        }
                                        /* match variable from list */
        vp = matchvar(varname);
        if (vp == NULL) {
                fprintf(stderr, "%s: unknown variable '%s'\n",
                                progname, varname);
                exit(1);
        }
                                        /* assign new value */
        if (i = vp->nass) {
                cp = vp->value;
                while (i--)
                        while (*cp++)
                                ;
                i = cp - vp->value;
                vp->value = realloc(vp->value, i+n+1);
        } else
                vp->value = malloc(n+1);
        if (vp->value == NULL)
                syserr(progname);
        strcpy(vp->value+i, ass);
        vp->nass++;
}


VARIABLE *
matchvar(nam)                   /* match a variable by its name */
char    *nam;
{
        int     n = strlen(nam);
        register int    i;

        for (i = 0; i < NVARS; i++)
                if (n >= vv[i].nick && !strncmp(nam, vv[i].name, n))
                        return(vv+i);
        return(NULL);
}


char *
nvalue(vp, n)                   /* return nth variable value */
VARIABLE        *vp;
register int    n;
{
        register char   *cp;

        if (vp == NULL || n < 0 || n >= vp->nass)
                return(NULL);
        cp = vp->value;
        while (n--)
                while (*cp++)
                        ;
        return(cp);
}


int
vscale(vc)                      /* return scale for variable vc */
int     vc;
{
        switch(vlet(vc)) {
        case 'H':
                return(HIGH);
        case 'M':
                return(MEDIUM);
        case 'L':
                return(LOW);
        }
        badvalue(vc);
}


checkvalues()                   /* check assignments */
{
        register int    i;

        for (i = 0; i < NVARS; i++)
                if (vv[i].fixval != NULL)
                        (*vv[i].fixval)(vv+i);
}


onevalue(vp)                    /* only one assignment for this variable */
register VARIABLE       *vp;
{
        if (vp->nass < 2)
                return;
        fprintf(stderr, "%s: warning - multiple assignment of variable '%s'\n",
                        progname, vp->name);
        do
                vp->value += strlen(vp->value)+1;
        while (--vp->nass > 1);
}


catvalues(vp)                   /* concatenate variable values */
register VARIABLE       *vp;
{
        register char   *cp;

        if (vp->nass < 2)
                return;
        for (cp = vp->value; vp->nass > 1; vp->nass--) {
                while (*cp)
                        cp++;
                *cp++ = ' ';
        }
}


long
checklast(fnames)                       /* check files and find most recent */
register char   *fnames;
{
        char    thisfile[MAXPATH];
        long    thisdate, lastdate = -1;
        register char   *cp;

        while (*fnames) {
                while (isspace(*fnames)) fnames++;
                cp = thisfile;
                while (*fnames && !isspace(*fnames))
                        *cp++ = *fnames++;
                *cp = '\0';
                if ((thisdate = fdate(thisfile)) < 0)
                        syserr(thisfile);
                if (thisdate > lastdate)
                        lastdate = thisdate;
        }
        return(lastdate);
}


checkfiles()                    /* check for existence and modified times */
{
        char    *cp;
        long    objdate;

        if (!vdef(OCTREE)) {
                if ((cp = bmalloc(strlen(radname)+5)) == NULL)
                        syserr(progname);
                sprintf(cp, "%s.oct", radname);
                vval(OCTREE) = cp;
                vdef(OCTREE)++;
        }
        octreedate = fdate(vval(OCTREE));
        scenedate = -1;
        if (vdef(SCENE)) {
                scenedate = checklast(vval(SCENE));
                if (vdef(OBJECT)) {
                        objdate = checklast(vval(OBJECT));
                        if (objdate > scenedate)
                                scenedate = objdate;
                }
        }
        if (octreedate < 0 & scenedate < 0) {
                fprintf(stderr, "%s: need '%s' or '%s'\n", progname,
                                vnam(OCTREE), vnam(SCENE));
                exit(1);
        }
}       


getoctcube(org, sizp)           /* get octree bounding cube */
double  org[3], *sizp;
{
        extern FILE     *popen();
        static double   oorg[3], osiz = 0.;
        char    buf[MAXPATH+16];
        FILE    *fp;

        if (osiz <= FTINY) {
                oconv();                /* does nothing if done already */
                sprintf(buf, "getinfo -d < %s", vval(OCTREE));
                if ((fp = popen(buf, "r")) == NULL)
                        syserr("getinfo");
                if (fscanf(fp, "%lf %lf %lf %lf", &oorg[0], &oorg[1],
                                &oorg[2], &osiz) != 4) {
                        fprintf(stderr,
                        "%s: error reading bounding cube from getinfo\n",
                                        progname);
                        exit(1);
                }
                pclose(fp);
        }
        org[0] = oorg[0]; org[1] = oorg[1]; org[2] = oorg[2]; *sizp = osiz;
}


setdefaults()                   /* set default values for unassigned var's */
{
        double  org[3], size;
        char    buf[128];

        if (!vdef(ZONE)) {
                getoctcube(org, &size);
                sprintf(buf, "E %g %g %g %g %g %g", org[0], org[0]+size,
                                org[1], org[1]+size, org[2], org[2]+size);
                vval(ZONE) = savqstr(buf);
                vdef(ZONE)++;
        }
        if (!vdef(INDIRECT)) {
                vval(INDIRECT) = "0";
                vdef(INDIRECT)++;
        }
        if (!vdef(QUALITY)) {
                vval(QUALITY) = "L";
                vdef(QUALITY)++;
        }
        if (!vdef(RESOLUTION)) {
                vval(RESOLUTION) = "512";
                vdef(RESOLUTION)++;
        }
        if (!vdef(PICTURE)) {
                vval(PICTURE) = radname;
                vdef(PICTURE)++;
        }
        if (!vdef(VIEW)) {
                vval(VIEW) = "X";
                vdef(VIEW)++;
        }
        if (!vdef(DETAIL)) {
                vval(DETAIL) = "M";
                vdef(DETAIL)++;
        }
        if (!vdef(PENUMBRAS)) {
                vval(PENUMBRAS) = "F";
                vdef(PENUMBRAS)++;
        }
        if (!vdef(VARIABILITY)) {
                vval(VARIABILITY) = "L";
                vdef(VARIABILITY)++;
        }
}


printvals()                     /* print variable values */
{
        register int    i, j;

        for (i = 0; i < NVARS; i++)
                for (j = 0; j < vdef(i); j++)
                        printf("%s= %s\n", vnam(i), nvalue(vv+i, j));
        fflush(stdout);
}


oconv()                         /* run oconv if necessary */
{
        char    combuf[512], ocopts[64];

        if (octreedate >= scenedate)    /* check dates */
                return;
                                        /* build command */
        oconvopts(ocopts);
        if (vdef(MATERIAL))
                sprintf(combuf, "oconv%s %s %s > %s", ocopts,
                                vval(MATERIAL), vval(SCENE), vval(OCTREE));
        else
                sprintf(combuf, "oconv%s %s > %s", ocopts,
                                vval(SCENE), vval(OCTREE));
        
        if (runcom(combuf)) {           /* run it */
                fprintf(stderr, "%s: error generating octree\n\t%s removed\n",
                                progname, vval(OCTREE));
                unlink(vval(OCTREE));
                exit(1);
        }
        octreedate = time(0);
}


char *
addarg(op, arg)                         /* add argument and advance pointer */
register char   *op, *arg;
{
        *op = ' ';
        while (*++op = *arg++)
                ;
        return(op);
}


oconvopts(oo)                           /* get oconv options */
register char   *oo;
{
        /* BEWARE:  This may be called via setdefaults(), so no assumptions */

        *oo = '\0';
        if (vdef(OCONV))
                addarg(oo, vval(OCONV));
}


double
ambval()                                /* compute ambient value */
{
        if (vdef(EXPOSURE)) {
                if (vval(EXPOSURE)[0] == '+' || vval(EXPOSURE)[0] == '-')
                        return(.5/pow(2.,atof(vval(EXPOSURE))));
                if (isdigit(vval(EXPOSURE)[0]) || vval(EXPOSURE)[0] == '.')
                        return(.5/atof(vval(EXPOSURE)));
                badvalue(EXPOSURE);
        }
        if (vlet(ZONE) == 'E')
                return(10.);
        if (vlet(ZONE) == 'I')
                return(.01);
        badvalue(ZONE);
}


lowqopts(op)                            /* low quality rendering options */
register char   *op;
{
        double  d, org[3], siz[3];

        *op = '\0';
        if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
                        &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
                badvalue(ZONE);
        siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
        getoctcube(org, &d);
        d *= 3./(siz[0]+siz[1]+siz[2]);
        switch (vscale(DETAIL)) {
        case LOW:
                op = addarg(op, "-ps 16 -dp 16");
                sprintf(op, " -ar %d", (int)(4*d));
                op += strlen(op);
                break;
        case MEDIUM:
                op = addarg(op, "-ps 8 -dp 32");
                sprintf(op, " -ar %d", (int)(8*d));
                op += strlen(op);
                break;
        case HIGH:
                op = addarg(op, "-ps 4 -dp 64");
                sprintf(op, " -ar %d", (int)(16*d));
                op += strlen(op);
                break;
        }
        op = addarg(op, "-pt .16");
        if (vbool(PENUMBRAS))
                op = addarg(op, "-ds .4");
        else
                op = addarg(op, "-ds 0");
        op = addarg(op, "-dt .2 -dc .25 -dr 0 -sj 0 -st .7");
        if (vdef(AMBFILE)) {
                sprintf(op, " -af %s", vval(AMBFILE));
                op += strlen(op);
        } else
                overture = 0;
        switch (vscale(VARIABILITY)) {
        case LOW:
                op = addarg(op, "-aa .4 -ad 32");
                break;
        case MEDIUM:
                op = addarg(op, "-aa .3 -ad 64");
                break;
        case HIGH:
                op = addarg(op, "-aa .25 -ad 128");
                break;
        }
        op = addarg(op, "-as 0");
        d = ambval();
        sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
        op += strlen(op);
        op = addarg(op, "-lr 3 -lw .02");
        if (vdef(RENDER))
                op = addarg(op, vval(RENDER));
}


medqopts(op)                            /* medium quality rendering options */
register char   *op;
{
        double  d, org[3], siz[3];

        *op = '\0';
        if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
                        &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
                badvalue(ZONE);
        siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
        getoctcube(org, &d);
        d *= 3./(siz[0]+siz[1]+siz[2]);
        switch (vscale(DETAIL)) {
        case LOW:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 4" : "-ps 8");
                op = addarg(op, "-dp 64");
                sprintf(op, " -ar %d", (int)(8*d));
                op += strlen(op);
                break;
        case MEDIUM:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 3" : "-ps 6");
                op = addarg(op, "-dp 128");
                sprintf(op, " -ar %d", (int)(16*d));
                op += strlen(op);
                break;
        case HIGH:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 2" : "-ps 4");
                op = addarg(op, "-dp 256");
                sprintf(op, " -ar %d", (int)(32*d));
                op += strlen(op);
                break;
        }
        op = addarg(op, "-pt .08");
        if (vbool(PENUMBRAS))
                op = addarg(op, "-ds .2 -dj .35");
        else
                op = addarg(op, "-ds .3");
        op = addarg(op, "-dt .1 -dc .5 -dr 1 -sj .7 -st .15");
        if (overture = vint(INDIRECT)) {
                sprintf(op, " -ab %d", overture);
                op += strlen(op);
        }
        if (vdef(AMBFILE)) {
                sprintf(op, " -af %s", vval(AMBFILE));
                op += strlen(op);
        } else
                overture = 0;
        switch (vscale(VARIABILITY)) {
        case LOW:
                op = addarg(op, "-aa .25 -ad 128 -as 0");
                break;
        case MEDIUM:
                op = addarg(op, "-aa .2 -ad 300 -as 64");
                break;
        case HIGH:
                op = addarg(op, "-aa .15 -ad 500 -as 128");
                break;
        }
        d = ambval();
        sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
        op += strlen(op);
        op = addarg(op, "-lr 6 -lw .002");
        if (vdef(RENDER))
                op = addarg(op, vval(RENDER));
}


hiqopts(op)                             /* high quality rendering options */
register char   *op;
{
        double  d, org[3], siz[3];

        *op = '\0';
        if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf", &org[0],
                        &siz[0], &org[1], &siz[1], &org[2], &siz[2]) != 6)
                badvalue(ZONE);
        siz[0] -= org[0]; siz[1] -= org[1]; siz[2] -= org[2];
        getoctcube(org, &d);
        d *= 3./(siz[0]+siz[1]+siz[2]);
        switch (vscale(DETAIL)) {
        case LOW:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 8");
                op = addarg(op, "-dp 256");
                sprintf(op, " -ar %d", (int)(16*d));
                op += strlen(op);
                break;
        case MEDIUM:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 5");
                op = addarg(op, "-dp 512");
                sprintf(op, " -ar %d", (int)(32*d));
                op += strlen(op);
                break;
        case HIGH:
                op = addarg(op, vbool(PENUMBRAS) ? "-ps 1" : "-ps 3");
                op = addarg(op, "-dp 1024");
                sprintf(op, " -ar %d", (int)(64*d));
                op += strlen(op);
                break;
        }
        op = addarg(op, "-pt .04");
        if (vbool(PENUMBRAS))
                op = addarg(op, "-ds .1 -dj .7");
        else
                op = addarg(op, "-ds .2");
        op = addarg(op, "-dt .05 -dc .75 -dr 3 -sj 1 -st .03");
        sprintf(op, " -ab %d", overture=vint(INDIRECT)+1);
        op += strlen(op);
        if (vdef(AMBFILE)) {
                sprintf(op, " -af %s", vval(AMBFILE));
                op += strlen(op);
        } else
                overture = 0;
        switch (vscale(VARIABILITY)) {
        case LOW:
                op = addarg(op, "-aa .15 -ad 200 -as 0");
                break;
        case MEDIUM:
                op = addarg(op, "-aa .125 -ad 512 -as 128");
                break;
        case HIGH:
                op = addarg(op, "-aa .08 -ad 850 -as 256");
                break;
        }
        d = ambval();
        sprintf(op, " -av %.2g %.2g %.2g", d, d, d);
        op += strlen(op);
        op = addarg(op, "-lr 12 -lw .0005");
        if (vdef(RENDER))
                op = addarg(op, vval(RENDER));
}


xferopts(ro)                            /* transfer options if indicated */
char    *ro;
{
        int     fd, n;
        
        n = strlen(ro);
        if (n < 2)
                return;
        if (vdef(OPTFILE)) {
                if ((fd = open(vval(OPTFILE), O_WRONLY|O_CREAT|O_TRUNC, 0666)) == -1)
                        syserr(vval(OPTFILE));
                if (write(fd, ro+1, n-1) != n-1)
                        syserr(vval(OPTFILE));
                write(fd, "\n", 1);
                close(fd);
                sprintf(ro, " \"^%s\"", vval(OPTFILE));
        }
#ifdef MSDOS
        else if (n > 50) {
                register char   *evp = bmalloc(n+6);
                if (evp == NULL)
                        syserr(progname);
                strcpy(evp, "ROPT=");
                strcat(evp, ro);
                if (putenv(evp) != 0) {
                        fprintf(stderr, "%s: out of environment space\n",
                                        progname);
                        exit(1);
                }
                strcpy(ro, " $ROPT");
        }
#endif
}


pfiltopts(po)                           /* get pfilt options */
register char   *po;
{
        *po = '\0';
        if (vdef(EXPOSURE)) {
                po = addarg(po, "-1 -e");
                po = addarg(po, vval(EXPOSURE));
        }
        if (vscale(QUALITY) == HIGH)
                po = addarg(po, "-r .65");
        if (vdef(PFILT))
                po = addarg(po, vval(PFILT));
}


matchword(s1, s2)                       /* match white-delimited words */
register char   *s1, *s2;
{
        while (isspace(*s1)) s1++;
        while (isspace(*s2)) s2++;
        while (*s1 && !isspace(*s1))
                if (*s1++ != *s2++)
                        return(0);
        return(!*s2 || isspace(*s2));
}


char *
specview(vs)                            /* get proper view spec from vs */
register char   *vs;
{
        static char     vup[7][12] = {"-vu 0 0 -1","-vu 0 -1 0","-vu -1 0 0",
                        "-vu 0 0 1", "-vu 1 0 0","-vu 0 1 0","-vu 0 0 1"};
        static char     viewopts[128];
        register char   *cp;
        int     xpos, ypos, zpos, viewtype, upax;
        register int    i;
        double  cent[3], dim[3], mult, d;

        if (vs == NULL || *vs == '-')
                return(vs);
        upax = 0;                       /* get the up vector */
        if (vdef(UP)) {
                if (vval(UP)[0] == '-' || vval(UP)[0] == '+')
                        upax = 1-'X'+UPPER(vval(UP)[1]);
                else
                        upax = 1-'X'+vlet(UP);
                if (upax < 1 | upax > 3)
                        badvalue(UP);
                if (vval(UP)[0] == '-')
                        upax = -upax;
        }
                                        /* check standard view names */
        xpos = ypos = zpos = 0;
        if (*vs == 'X') {
                xpos = 1; vs++;
        } else if (*vs == 'x') {
                xpos = -1; vs++;
        }
        if (*vs == 'Y') {
                ypos = 1; vs++;
        } else if (*vs == 'y') {
                ypos = -1; vs++;
        }
        if (*vs == 'Z') {
                zpos = 1; vs++;
        } else if (*vs == 'z') {
                zpos = -1; vs++;
        }
        viewtype = 'v';
        if (*vs == 'v' | *vs == 'l' | *vs == 'a' | *vs == 'h')
                viewtype = *vs++;
        cp = viewopts;
        if ((!*vs || isspace(*vs)) && (xpos|ypos|zpos)) {       /* got one! */
                *cp++ = '-'; *cp++ = 'v'; *cp++ = 't'; *cp++ = viewtype;
                if (sscanf(vval(ZONE), "%*s %lf %lf %lf %lf %lf %lf",
                                &cent[0], &dim[0], &cent[1], &dim[1],
                                &cent[2], &dim[2]) != 6)
                        badvalue(ZONE);
                for (i = 0; i < 3; i++) {
                        dim[i] -= cent[i];
                        cent[i] += .5*dim[i];
                }
                mult = vlet(ZONE)=='E' ? 2. : .45 ;
                sprintf(cp, " -vp %.2g %.2g %.2g -vd %.2g %.2g %.2g",
                                cent[0]+xpos*mult*dim[0],
                                cent[1]+ypos*mult*dim[1],
                                cent[2]+zpos*mult*dim[2],
                                -xpos*dim[0], -ypos*dim[1], -zpos*dim[2]);
                cp += strlen(cp);
                                        /* redirect up axis if necessary */
                switch (upax) {
                case 3:                 /* plus or minus Z axis */
                case -3:
                case 0:
                        if (!(xpos|ypos))
                                upax = 2;
                        break;
                case 2:                 /* plus or minus Y axis */
                case -2:
                        if (!(xpos|zpos))
                                upax = 1;
                        break;
                case 1:                 /* plus or minus X axis */
                case -1:
                        if (!(ypos|zpos))
                                upax = 3;
                        break;
                }
                cp = addarg(cp, vup[upax+3]);
                switch (viewtype) {
                case 'v':
                        cp = addarg(cp, "-vh 45 -vv 45");
                        break;
                case 'l':
                        d = sqrt(dim[0]*dim[0]+dim[1]*dim[1]+dim[2]*dim[2]);
                        sprintf(cp, " -vh %.2g -vv %.2g", d, d);
                        cp += strlen(cp);
                        break;
                case 'a':
                case 'h':
                        cp = addarg(cp, "-vh 180 -vv 180");
                        break;
                }
        } else {
                while (!isspace(*vs))           /* else skip id */
                        if (!*vs++)
                                return(NULL);
                if (upax) {                     /* specify up vector */
                        strcpy(cp, vup[upax+3]);
                        cp += strlen(cp);
                }
        }
                                        /* append any additional options */
        strcpy(cp, vs);
        return(viewopts);
}


char *
getview(n, vn)                          /* get view n, or NULL if none */
int     n;
char    *vn;
{
        register char   *mv;

        if (viewselect != NULL) {
                if (n)                          /* only do one */
                        return(NULL);
                if (viewselect[0] == '-') {     /* already specified */
                        if (vn != NULL) *vn = '\0';
                        return(viewselect);
                }
                if (vn != NULL) {
                        for (mv = viewselect; *mv && !isspace(*mv);
                                        *vn++ = *mv++)
                                ;
                        *vn = '\0';
                }
                                                /* check list */
                while ((mv = nvalue(vv+VIEW, n++)) != NULL)
                        if (matchword(viewselect, mv))
                                return(specview(mv));
                return(specview(viewselect));   /* standard view? */
        }
        if (vn != NULL && (mv = nvalue(vv+VIEW, n)) != NULL) {
                if (*mv != '-')
                        while (*mv && !isspace(*mv))
                                *vn++ = *mv++;
                *vn = '\0';
        }
        return(specview(nvalue(vv+VIEW, n)));   /* use view n */
}


rview(opts)                             /* run rview with first view */
char    *opts;
{
        char    combuf[512];
                                        /* build command */
        sprintf(combuf, "rview %s%s ", getview(0, NULL), opts);
        if (rvdevice != NULL)
                sprintf(combuf+strlen(combuf), "-o %s ", rvdevice);
        strcat(combuf, vval(OCTREE));
        if (runcom(combuf)) {           /* run it */
                fprintf(stderr, "%s: error running rview\n", progname);
                exit(1);
        }
}


rpict(opts)                             /* run rpict and pfilt for each view */
char    *opts;
{
        char    combuf[1024];
        char    rawfile[MAXPATH], picfile[MAXPATH], rep[MAXPATH+16], res[32];
        char    pfopts[128];
        char    vs[32], *vw;
        int     vn, mult;
                                        /* get pfilt options */
        pfiltopts(pfopts);
                                        /* get resolution, reporting */
        mult = vscale(QUALITY)+1;
        {
                int     xres, yres;
                double  aspect;
                int     n;
                n = sscanf(vval(RESOLUTION), "%d %d %lf", &xres, &yres, &aspect);
                if (n == 3)
                        sprintf(res, "-x %d -y %d -pa %.3f",
                                        mult*xres, mult*yres, aspect);
                else if (n) {
                        if (n == 1) yres = xres;
                        sprintf(res, "-x %d -y %d", mult*xres, mult*yres);
                } else
                        badvalue(RESOLUTION);
        }
        rep[0] = '\0';
        if (vdef(REPORT)) {
                double  minutes;
                int     n;
                n = sscanf(vval(REPORT), "%lf %s", &minutes, rawfile);
                if (n == 2)
                        sprintf(rep, " -t %d -e %s", (int)(minutes*60), rawfile);
                else if (n == 1)
                        sprintf(rep, " -t %d", (int)(minutes*60));
                else
                        badvalue(REPORT);
        }
                                        /* check date on ambient file */
        if (vdef(AMBFILE)) {
                long    afdate = fdate(vval(AMBFILE));
                if (afdate >= 0 & octreedate > afdate)
                        rmfile(vval(AMBFILE));
        }
                                        /* do each view */
        vn = 0;
        while ((vw = getview(vn++, vs)) != NULL) {
                if (!vs[0])
                        sprintf(vs, "%d", vn);
                sprintf(picfile, "%s_%s.pic", vval(PICTURE), vs);
                                                /* check date on picture */
                if (fdate(picfile) > octreedate)
                        continue;
                                                /* build rpict command */
                sprintf(rawfile, "%s_%s.raw", vval(PICTURE), vs);
                if (fdate(rawfile) > octreedate)        /* recover */
                        sprintf(combuf, "rpict%s%s -ro %s %s",
                                        rep, opts, rawfile, vval(OCTREE));
                else {
                        if (overture) {         /* run overture calculation */
                                sprintf(combuf,
                                "rpict%s %s%s -x 64 -y 64 -ps 1 %s > %s",
                                                rep, vw, opts,
                                                vval(OCTREE), overfile);
                                if (runcom(combuf)) {
                                        fprintf(stderr,
                                        "%s: error in overture for view %s\n",
                                                progname, vs);
                                        exit(1);
                                }
#ifdef NIX
                                rmfile(overfile);
#endif
                        }
                        sprintf(combuf, "rpict%s %s %s%s %s > %s",
                                        rep, vw, res, opts,
                                        vval(OCTREE), rawfile);
                }
                if (runcom(combuf)) {           /* run rpict */
                        fprintf(stderr, "%s: error rendering view %s\n",
                                        progname, vs);
                        exit(1);
                }
                                                /* build pfilt command */
                if (mult > 1)
                        sprintf(combuf, "pfilt%s -x /%d -y /%d %s > %s",
                                        pfopts, mult, mult, rawfile, picfile);
                else
                        sprintf(combuf, "pfilt%s %s > %s", pfopts,
                                        rawfile, picfile);
                if (runcom(combuf)) {           /* run pfilt */
                        fprintf(stderr,
                        "%s: error filtering view %s\n\t%s removed\n",
                                        progname, vs, picfile);
                        unlink(picfile);
                        exit(1);
                }
                                                /* remove raw file */
                rmfile(rawfile);
        }
}


runcom(cs)                      /* run command */
char    *cs;
{
        if (!silent)            /* echo it */
                printf("\t%s\n", cs);
        if (noaction)
                return(0);
        fflush(stdout);         /* flush output and pass to shell */
        return(system(cs));
}


rmfile(fn)                      /* remove a file */
char    *fn;
{
        if (!silent)
#ifdef MSDOS
                printf("\tdel %s\n", fn);
#else
                printf("\trm -f %s\n", fn);
#endif
        if (noaction)
                return(0);
        return(unlink(fn));
}


badvalue(vc)                    /* report bad variable value and exit */
int     vc;
{
        fprintf(stderr, "%s: bad value for variable '%s'\n",
                        progname, vnam(vc));
        exit(1);
}


syserr(s)                       /* report a system error and exit */
char    *s;
{
        perror(s);
        exit(1);
}
Revision:	2.6
Committed:	Tue Mar 30 13:16:26 1993 UTC (31 years, 1 month ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.5:	+13 -4 lines
Log Message:	fixed potential problem in failed overture calculation (recovery!)