/* 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);
}