src/rt/rv2.c

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

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

/*
 *  rv2.c - command routines used in tracing a view.
 *
 *     3/24/87
 */

#include  "ray.h"

#include  "octree.h"

#include  "otypes.h"

#include  "rpaint.h"

#include  <ctype.h>

#define  CTRL(c)        ('c'-'@')

extern char  VersionID[];
extern char  *progname;
extern char  *octname;


getframe(s)                             /* get a new frame */
char  *s;
{
        if (getrect(s, &pframe) < 0)
                return;
        pdepth = 0;
}


getrepaint(s)                           /* get area and repaint */
char  *s;
{
        RECT  box;

        if (getrect(s, &box) < 0)
                return;
        paintrect(&ptrunk, 0, 0, hresolu, vresolu, &box);
}


getview(s)                              /* get/show view parameters */
char  *s;
{
        FILE  *fp;
        char  buf[128];
        char  *fname;
        int  change = 0;
        VIEW  nv;

        if (sscanf(s, "%s", buf) == 1) {        /* write parameters to a file */
                if ((fname = getpath(buf, NULL, 0)) == NULL ||
                                (fp = fopen(fname, "a")) == NULL) {
                        sprintf(errmsg, "cannot open \"%s\"", buf);
                        error(COMMAND, errmsg);
                        return;
                }
                fputs(progname, fp);
                fprintview(&ourview, fp);
                fputs(sskip(s), fp);
                fputs("\n", fp);
                fclose(fp);
                return;
        }
        sprintf(buf, "view type (%c): ", ourview.type);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (buf[0] && buf[0] != ourview.type) {
                nv.type = buf[0];
                change++;
        } else
                nv.type = ourview.type;
        sprintf(buf, "view point (%.6g %.6g %.6g): ",
                        ourview.vp[0], ourview.vp[1], ourview.vp[2]);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (sscanf(buf, "%lf %lf %lf", &nv.vp[0], &nv.vp[1], &nv.vp[2]) == 3)
                change++;
        else
                VCOPY(nv.vp, ourview.vp);
        sprintf(buf, "view direction (%.6g %.6g %.6g): ",
                        ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (sscanf(buf,"%lf %lf %lf",&nv.vdir[0],&nv.vdir[1],&nv.vdir[2]) == 3)
                change++;
        else
                VCOPY(nv.vdir, ourview.vdir);
        sprintf(buf, "view up (%.6g %.6g %.6g): ",
                        ourview.vup[0], ourview.vup[1], ourview.vup[2]);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (sscanf(buf,"%lf %lf %lf",&nv.vup[0],&nv.vup[1],&nv.vup[2]) == 3)
                change++;
        else
                VCOPY(nv.vup, ourview.vup);
        sprintf(buf, "view horiz and vert size (%.6g %.6g): ",
                        ourview.horiz, ourview.vert);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2)
                change++;
        else {
                nv.horiz = ourview.horiz; nv.vert = ourview.vert;
        }
        sprintf(buf, "view shift and lift (%.6g %.6g): ",
                        ourview.hoff, ourview.voff);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL(C)) return;
        if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2)
                change++;
        else {
                nv.hoff = ourview.hoff; nv.voff = ourview.voff;
        }
        if (change)
                newview(&nv);
}


lastview(s)                             /* return to a previous view */
char  *s;
{
        char  buf[128];
        char  *fname;
        int  success;
        VIEW  nv;

        if (sscanf(s, "%s", buf) == 1) {        /* get parameters from a file */
                copystruct(&nv, &stdview);
                if ((fname = getpath(buf, NULL, 0)) == NULL ||
                                (success = viewfile(fname, &nv, 0, 0)) == -1) {
                        sprintf(errmsg, "cannot open \"%s\"", buf);
                        error(COMMAND, errmsg);
                        return;
                }
                if (!success)
                        error(COMMAND, "wrong file format");
                else
                        newview(&nv);
                return;
        }
        if (oldview.type == 0) {        /* no old view! */
                error(COMMAND, "no previous view");
                return;
        }
        copystruct(&nv, &ourview);
        copystruct(&ourview, &oldview);
        copystruct(&oldview, &nv);
        newimage();
}


getaim(s)                               /* aim camera */
char  *s;
{
        extern double  tan(), atan();
        double  zfact;
        VIEW  nv;

        if (getinterest(s, 1, nv.vdir, &zfact) < 0)
                return;
        nv.type = ourview.type;
        VCOPY(nv.vp, ourview.vp);
        VCOPY(nv.vup, ourview.vup);
        nv.hoff = ourview.hoff; nv.voff = ourview.voff;
        nv.horiz = ourview.horiz; nv.vert = ourview.vert;
        zoomview(&nv, zfact);
        newview(&nv);
}


getmove(s)                              /* move camera */
char  *s;
{
        FVECT  vc;
        double  mag;

        if (getinterest(s, 0, vc, &mag) < 0)
                return;
        moveview(0.0, 0.0, mag, vc);
}


getrotate(s)                            /* rotate camera */
char  *s;
{
        extern double  normalize(), tan(), atan();
        VIEW  nv;
        FVECT  v1;
        double  angle, elev, zfact;
        
        elev = 0.0; zfact = 1.0;
        if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) {
                error(COMMAND, "missing angle");
                return;
        }
        nv.type = ourview.type;
        VCOPY(nv.vp, ourview.vp);
        VCOPY(nv.vup, ourview.vup);
        nv.hoff = ourview.hoff; nv.voff = ourview.voff;
        spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
        if (elev != 0.0) {
                fcross(v1, nv.vdir, ourview.vup);
                normalize(v1);
                spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
        }
        nv.horiz = ourview.horiz; nv.vert = ourview.vert;
        zoomview(&nv, zfact);
        newview(&nv);
}


getpivot(s)                             /* pivot viewpoint */
register char  *s;
{
        FVECT  vc;
        double  angle, elev, mag;

        elev = 0.0;
        if (sscanf(s, "%lf %lf", &angle, &elev) < 1) {
                error(COMMAND, "missing angle");
                return;
        }
        if (getinterest(sskip(sskip(s)), 0, vc, &mag) < 0)
                return;
        moveview(angle, elev, mag, vc);
}


getexposure(s)                          /* get new exposure */
char  *s;
{
        double  atof(), pow(), fabs();
        char  buf[128];
        register char  *cp;
        register PNODE  *p;
        RECT  r;
        int  x, y;
        double  e;

        for (cp = s; isspace(*cp); cp++)
                ;
        if (*cp == '\0') {              /* normalize to point */
                if (dev->getcur == NULL)
                        return;
                (*dev->comout)("Pick point for exposure\n");
                if ((*dev->getcur)(&x, &y) == ABORT)
                        return;
                r.l = r.d = 0;
                r.r = hresolu; r.u = vresolu;
                p = findrect(x, y, &ptrunk, &r, -1);
                e = 1.0;
        } else {
                if (*cp == '=') {       /* absolute setting */
                        p = NULL;
                        e = 1.0/exposure;
                        for (cp++; isspace(*cp); cp++)
                                ;
                        if (*cp == '\0') {      /* interactive */
                                sprintf(buf, "exposure (%lf): ", exposure);
                                (*dev->comout)(buf);
                                (*dev->comin)(buf, NULL);
                                for (cp = buf; isspace(*cp); cp++)
                                        ;
                                if (*cp == '\0')
                                        return;
                        }
                } else {                /* normalize to average */
                        p = &ptrunk;
                        e = 1.0;
                }
                if (*cp == '+' || *cp == '-')   /* f-stops */
                        e *= pow(2.0, atof(cp));
                else                            /* multiplier */
                        e *= atof(cp);
        }
        if (p != NULL) {                /* relative setting */
                if (bright(p->v) <= FTINY) {
                        error(COMMAND, "cannot normalize to zero");
                        return;
                }
                e *= 0.5 / bright(p->v);
        }
        if (e <= FTINY || fabs(1.0 - e) <= FTINY)
                return;
        scalepict(&ptrunk, e);
        exposure *= e;
        redraw();
}


getparam(str, dsc, typ, ptr)            /* get variable from user */
char  *str, *dsc;
int  typ;
register union {int i; double d; COLOR C;}  *ptr;
{
        extern char  *index();
        int  i0;
        double  d0, d1, d2;
        char  buf[48];

        switch (typ) {
        case 'i':                       /* integer */
                if (sscanf(str, "%d", &i0) != 1) {
                        (*dev->comout)(dsc);
                        sprintf(buf, " (%d): ", ptr->i);
                        (*dev->comout)(buf);
                        (*dev->comin)(buf, NULL);
                        if (sscanf(buf, "%d", &i0) != 1)
                                break;
                }
                ptr->i = i0;
                break;
        case 'r':                       /* real */
                if (sscanf(str, "%lf", &d0) != 1) {
                        (*dev->comout)(dsc);
                        sprintf(buf, " (%.6g): ", ptr->d);
                        (*dev->comout)(buf);
                        (*dev->comin)(buf, NULL);
                        if (sscanf(buf, "%lf", &d0) != 1)
                                break;
                }
                ptr->d = d0;
                break;
        case 'b':                       /* boolean */
                if (sscanf(str, "%1s", buf) != 1) {
                        (*dev->comout)(dsc);
                        sprintf(buf, " (%c): ", ptr->i ? 'y' : 'n');
                        (*dev->comout)(buf);
                        (*dev->comin)(buf, NULL);
                        if (buf[0] == '\0' ||
                                        index("yY+1tTnN-0fF", buf[0]) == NULL)
                                break;
                }
                ptr->i = index("yY+1tT", buf[0]) != NULL;
                break;
        case 'C':                       /* color */
                if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
                        (*dev->comout)(dsc);
                        sprintf(buf, " (%.6g %.6g %.6g): ",
                                        colval(ptr->C,RED),
                                        colval(ptr->C,GRN),
                                        colval(ptr->C,BLU));
                        (*dev->comout)(buf);
                        (*dev->comin)(buf, NULL);
                        if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
                                break;
                }
                setcolor(ptr->C, d0, d1, d2);
                break;
        }
}


setparam(s)                             /* get/set program parameter */
register char  *s;
{
        extern int  psample;
        extern double  maxdiff;
        extern double  minweight;
        extern int  maxdepth;
        extern double  dstrsrc;
        extern double  shadthresh;
        extern double  shadcert;
        extern COLOR  ambval;
        extern double  ambacc;
        extern double  minarad;
        extern int  ambres;
        extern int  ambdiv;
        extern int  ambssamp;
        extern int  ambounce;
        extern int  directinvis;
        extern int  do_irrad;
        char  buf[128];
        
        if (s[0] == '\0') {
                (*dev->comout)(
                        "aa ab ad ar as av b dc di dj dt i lr lw sp st: ");
                (*dev->comin)(buf, NULL);
                s = buf;
        }
        switch (s[0]) {
        case 'l':                       /* limit */
                switch (s[1]) {
                case 'w':                       /* weight */
                        getparam(s+2, "limit weight", 'r', &minweight);
                        break;
                case 'r':                       /* reflection */
                        getparam(s+2, "limit reflection", 'i', &maxdepth);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'd':                       /* direct */
                switch (s[1]) {
                case 'j':                       /* jitter */
                        getparam(s+2, "direct jitter", 'r', &dstrsrc);
                        break;
                case 'c':                       /* certainty */
                        getparam(s+2, "direct certainty", 'r', &shadcert);
                        break;
                case 't':                       /* threshold */
                        getparam(s+2, "direct threshold", 'r', &shadthresh);
                        break;
                case 'i':                       /* invisibility */
                        getparam(s+2, "direct invisibility",
                                        'b', &directinvis);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'b':                       /* black and white */
                getparam(s+1, "black and white", 'b', &greyscale);
                break;
        case 'i':                       /* irradiance */
                getparam(s+1, "irradiance", 'b', &do_irrad);
                break;
        case 'a':                       /* ambient */
                switch (s[1]) {
                case 'v':                       /* value */
                        getparam(s+2, "ambient value", 'C', ambval);
                        break;
                case 'a':                       /* accuracy */
                        getparam(s+2, "ambient accuracy", 'r', &ambacc);
                        break;
                case 'd':                       /* divisions */
                        getparam(s+2, "ambient divisions", 'i', &ambdiv);
                        break;
                case 's':                       /* samples */
                        getparam(s+2, "ambient super-samples", 'i', &ambssamp);
                        break;
                case 'b':                       /* bounces */
                        getparam(s+2, "ambient bounces", 'i', &ambounce);
                        break;
                case 'r':
                        getparam(s+2, "ambient resolution", 'i', &ambres);
                        minarad = ambres > 0 ? thescene.cusize/ambres : 0.0;
                        break;
                default:
                        goto badparam;
                }
                break;
        case 's':                       /* sample */
                switch (s[1]) {
                case 'p':                       /* pixel */
                        getparam(s+2, "sample pixel", 'i', &psample);
                        pdepth = 0;
                        break;
                case 't':                       /* threshold */
                        getparam(s+2, "sample threshold", 'r', &maxdiff);
                        pdepth = 0;
                        break;
                default:
                        goto badparam;
                }
                break;
        case '\0':                      /* nothing */
                break;
        default:;
badparam:
                *sskip(s) = '\0';
                sprintf(errmsg, "%s: unknown variable", s);
                error(COMMAND, errmsg);
                break;
        }
}


traceray(s)                             /* trace a single ray */
char  *s;
{
        char  buf[128];
        int  x, y;
        RAY  thisray;

        if (sscanf(s, "%lf %lf %lf %lf %lf %lf",
                &thisray.rorg[0], &thisray.rorg[1], &thisray.rorg[2],
                &thisray.rdir[0], &thisray.rdir[1], &thisray.rdir[2]) != 6) {

                if (dev->getcur == NULL)
                        return;
                (*dev->comout)("Pick ray\n");
                if ((*dev->getcur)(&x, &y) == ABORT)
                        return;

                if (viewray(thisray.rorg, thisray.rdir, &ourview,
                                (x+.5)/hresolu, (y+.5)/vresolu) < 0) {
                        error(COMMAND, "not on image");
                        return;
                }

        } else if (normalize(thisray.rdir) == 0.0) {
                error(COMMAND, "zero ray direction");
                return;
        }

        rayorigin(&thisray, NULL, PRIMARY, 1.0);
        
        rayvalue(&thisray);

        if (thisray.ro == NULL)
                (*dev->comout)("ray hit nothing");
        else {
                sprintf(buf, "ray hit %s %s \"%s\"",
                                objptr(thisray.ro->omod)->oname,
                                ofun[thisray.ro->otype].funame,
                                thisray.ro->oname);
                (*dev->comout)(buf);
                (*dev->comin)(buf, NULL);
                if (thisray.rot >= FHUGE)
                        (*dev->comout)("at infinity");
                else {
                        sprintf(buf, "at (%.6g %.6g %.6g)", thisray.rop[0],
                                        thisray.rop[1], thisray.rop[2]);
                        (*dev->comout)(buf);
                }
                (*dev->comin)(buf, NULL);
                sprintf(buf, "with value (%.6g %.6g %.6g)",
                                colval(thisray.rcol,RED),
                                colval(thisray.rcol,GRN),
                                colval(thisray.rcol,BLU));
                (*dev->comout)(buf);
        }
        (*dev->comin)(buf, NULL);
}


writepict(s)                            /* write the picture to a file */
char  *s;
{
        static char  buf[128];
        char  *fname;
        FILE  *fp;
        COLR  *scanline;
        int  y;

        if (sscanf(s, "%s", buf) != 1 && buf[0] == '\0') {
                error(COMMAND, "no file");
                return;
        }
        if ((fname = getpath(buf, NULL, 0)) == NULL ||
                        (fp = fopen(fname, "w")) == NULL) {
                sprintf(errmsg, "cannot open \"%s\"", buf);
                error(COMMAND, errmsg);
                return;
        }
        (*dev->comout)("writing \"");
        (*dev->comout)(fname);
        (*dev->comout)("\"...\n");
                                                /* write header */
        fputs(progname, fp);
        fprintview(&ourview, fp);
        if (octname != NULL)
                fprintf(fp, " %s\n", octname);
        else
                putc('\n', fp);
        fprintf(fp, "SOFTWARE= %s\n", VersionID);
        if (exposure != 1.0)
                fputexpos(exposure, fp);
        if (dev->pixaspect != 1.0)
                fputaspect(dev->pixaspect, fp);
        fputformat(COLRFMT, fp);
        putc('\n', fp);
        fputresolu(YMAJOR|YDECR, hresolu, vresolu, fp);

        scanline = (COLR *)malloc(hresolu*sizeof(COLR));
        if (scanline == NULL) {
                error(COMMAND, "not enough memory!");
                fclose(fp);
                unlink(fname);
                return;
        }
        for (y = vresolu-1; y >= 0; y--) {
                getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu);
                if (fwritecolrs(scanline, hresolu, fp) < 0)
                        break;
        }
        free((char *)scanline);
        if (fclose(fp) < 0)
                error(COMMAND, "write error");
}
Revision:	1.29
Committed:	Fri Oct 11 12:12:40 1991 UTC (32 years, 6 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 1.28:	+6 -2 lines
Log Message:	added octree name to picture header
#	User	Rev	Content
1	greg	1.26	/* Copyright (c) 1991 Regents of the University of California */
2	greg	1.1
3			#ifndef lint
4			static char SCCSid[] = "$SunId$ LBL";
5			#endif
6
7			/*
8			* rv2.c - command routines used in tracing a view.
9			*
10			* 3/24/87
11			*/
12
13			#include "ray.h"
14
15			#include "octree.h"
16
17			#include "otypes.h"
18
19			#include "rpaint.h"
20
21			#include <ctype.h>
22
23			#define CTRL(c) ('c'-'@')
24
25	greg	1.29	extern char VersionID[];
26	greg	1.1	extern char *progname;
27	greg	1.29	extern char *octname;
28	greg	1.1
29
30			getframe(s) /* get a new frame */
31			char *s;
32			{
33	greg	1.12	if (getrect(s, &pframe) < 0)
34	greg	1.1	return;
35			pdepth = 0;
36			}
37
38
39	greg	1.12	getrepaint(s) /* get area and repaint */
40			char *s;
41			{
42			RECT box;
43
44			if (getrect(s, &box) < 0)
45			return;
46	greg	1.14	paintrect(&ptrunk, 0, 0, hresolu, vresolu, &box);
47	greg	1.12	}
48
49
50	greg	1.1	getview(s) /* get/show view parameters */
51			char *s;
52			{
53			FILE *fp;
54			char buf[128];
55	greg	1.10	char *fname;
56	greg	1.3	int change = 0;
57	greg	1.1	VIEW nv;
58
59			if (sscanf(s, "%s", buf) == 1) { /* write parameters to a file */
60	greg	1.10	if ((fname = getpath(buf, NULL, 0)) == NULL \|\|
61	greg	1.1	(fp = fopen(fname, "a")) == NULL) {
62			sprintf(errmsg, "cannot open \"%s\"", buf);
63			error(COMMAND, errmsg);
64			return;
65			}
66	greg	1.2	fputs(progname, fp);
67	greg	1.1	fprintview(&ourview, fp);
68	greg	1.19	fputs(sskip(s), fp);
69	greg	1.1	fputs("\n", fp);
70			fclose(fp);
71			return;
72			}
73			sprintf(buf, "view type (%c): ", ourview.type);
74			(*dev->comout)(buf);
75	greg	1.18	(*dev->comin)(buf, NULL);
76	greg	1.1	if (buf[0] == CTRL(C)) return;
77	greg	1.3	if (buf[0] && buf[0] != ourview.type) {
78	greg	1.1	nv.type = buf[0];
79	greg	1.3	change++;
80			} else
81	greg	1.1	nv.type = ourview.type;
82			sprintf(buf, "view point (%.6g %.6g %.6g): ",
83			ourview.vp[0], ourview.vp[1], ourview.vp[2]);
84			(*dev->comout)(buf);
85	greg	1.18	(*dev->comin)(buf, NULL);
86	greg	1.1	if (buf[0] == CTRL(C)) return;
87	greg	1.3	if (sscanf(buf, "%lf %lf %lf", &nv.vp[0], &nv.vp[1], &nv.vp[2]) == 3)
88			change++;
89			else
90	greg	1.1	VCOPY(nv.vp, ourview.vp);
91			sprintf(buf, "view direction (%.6g %.6g %.6g): ",
92			ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
93			(*dev->comout)(buf);
94	greg	1.18	(*dev->comin)(buf, NULL);
95	greg	1.1	if (buf[0] == CTRL(C)) return;
96	greg	1.3	if (sscanf(buf,"%lf %lf %lf",&nv.vdir[0],&nv.vdir[1],&nv.vdir[2]) == 3)
97			change++;
98			else
99	greg	1.1	VCOPY(nv.vdir, ourview.vdir);
100			sprintf(buf, "view up (%.6g %.6g %.6g): ",
101			ourview.vup[0], ourview.vup[1], ourview.vup[2]);
102			(*dev->comout)(buf);
103	greg	1.18	(*dev->comin)(buf, NULL);
104	greg	1.1	if (buf[0] == CTRL(C)) return;
105	greg	1.3	if (sscanf(buf,"%lf %lf %lf",&nv.vup[0],&nv.vup[1],&nv.vup[2]) == 3)
106			change++;
107			else
108	greg	1.1	VCOPY(nv.vup, ourview.vup);
109			sprintf(buf, "view horiz and vert size (%.6g %.6g): ",
110			ourview.horiz, ourview.vert);
111			(*dev->comout)(buf);
112	greg	1.18	(*dev->comin)(buf, NULL);
113	greg	1.1	if (buf[0] == CTRL(C)) return;
114	greg	1.3	if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2)
115			change++;
116			else {
117	greg	1.1	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
118			}
119	greg	1.14	sprintf(buf, "view shift and lift (%.6g %.6g): ",
120			ourview.hoff, ourview.voff);
121	greg	1.1	(*dev->comout)(buf);
122	greg	1.18	(*dev->comin)(buf, NULL);
123	greg	1.1	if (buf[0] == CTRL(C)) return;
124	greg	1.14	if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2)
125	greg	1.3	change++;
126			else {
127	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
128	greg	1.1	}
129	greg	1.3	if (change)
130			newview(&nv);
131	greg	1.1	}
132
133
134			lastview(s) /* return to a previous view */
135			char *s;
136			{
137			char buf[128];
138	greg	1.10	char *fname;
139	greg	1.1	int success;
140			VIEW nv;
141
142			if (sscanf(s, "%s", buf) == 1) { /* get parameters from a file */
143	greg	1.17	copystruct(&nv, &stdview);
144	greg	1.10	if ((fname = getpath(buf, NULL, 0)) == NULL \|\|
145	greg	1.16	(success = viewfile(fname, &nv, 0, 0)) == -1) {
146	greg	1.1	sprintf(errmsg, "cannot open \"%s\"", buf);
147			error(COMMAND, errmsg);
148			return;
149			}
150			if (!success)
151			error(COMMAND, "wrong file format");
152			else
153			newview(&nv);
154			return;
155			}
156	greg	1.20	if (oldview.type == 0) { /* no old view! */
157	greg	1.1	error(COMMAND, "no previous view");
158			return;
159			}
160	greg	1.17	copystruct(&nv, &ourview);
161			copystruct(&ourview, &oldview);
162			copystruct(&oldview, &nv);
163	greg	1.1	newimage();
164			}
165
166
167			getaim(s) /* aim camera */
168			char *s;
169			{
170			extern double tan(), atan();
171			double zfact;
172			VIEW nv;
173
174			if (getinterest(s, 1, nv.vdir, &zfact) < 0)
175			return;
176	greg	1.21	nv.type = ourview.type;
177	greg	1.1	VCOPY(nv.vp, ourview.vp);
178			VCOPY(nv.vup, ourview.vup);
179	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
180	greg	1.21	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
181			zoomview(&nv, zfact);
182	greg	1.1	newview(&nv);
183			}
184
185
186			getmove(s) /* move camera */
187			char *s;
188			{
189			FVECT vc;
190			double mag;
191
192			if (getinterest(s, 0, vc, &mag) < 0)
193			return;
194	greg	1.4	moveview(0.0, 0.0, mag, vc);
195	greg	1.1	}
196
197
198			getrotate(s) /* rotate camera */
199			char *s;
200			{
201	greg	1.4	extern double normalize(), tan(), atan();
202	greg	1.1	VIEW nv;
203			FVECT v1;
204	greg	1.4	double angle, elev, zfact;
205	greg	1.1
206	greg	1.4	elev = 0.0; zfact = 1.0;
207			if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) {
208	greg	1.1	error(COMMAND, "missing angle");
209			return;
210			}
211	greg	1.21	nv.type = ourview.type;
212	greg	1.1	VCOPY(nv.vp, ourview.vp);
213			VCOPY(nv.vup, ourview.vup);
214	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
215	greg	1.1	spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
216			if (elev != 0.0) {
217			fcross(v1, nv.vdir, ourview.vup);
218			normalize(v1);
219			spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
220			}
221	greg	1.21	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
222			zoomview(&nv, zfact);
223	greg	1.1	newview(&nv);
224			}
225
226
227			getpivot(s) /* pivot viewpoint */
228			register char *s;
229			{
230			FVECT vc;
231	greg	1.4	double angle, elev, mag;
232	greg	1.1
233	greg	1.4	elev = 0.0;
234			if (sscanf(s, "%lf %lf", &angle, &elev) < 1) {
235	greg	1.1	error(COMMAND, "missing angle");
236			return;
237			}
238	greg	1.4	if (getinterest(sskip(sskip(s)), 0, vc, &mag) < 0)
239	greg	1.1	return;
240	greg	1.4	moveview(angle, elev, mag, vc);
241	greg	1.1	}
242
243
244			getexposure(s) /* get new exposure */
245			char *s;
246			{
247			double atof(), pow(), fabs();
248			char buf[128];
249			register char *cp;
250			register PNODE *p;
251			RECT r;
252			int x, y;
253			double e;
254
255			for (cp = s; isspace(*cp); cp++)
256			;
257			if (cp == '\0') { / normalize to point */
258			if (dev->getcur == NULL)
259			return;
260			(*dev->comout)("Pick point for exposure\n");
261			if ((*dev->getcur)(&x, &y) == ABORT)
262			return;
263			r.l = r.d = 0;
264	greg	1.14	r.r = hresolu; r.u = vresolu;
265	greg	1.1	p = findrect(x, y, &ptrunk, &r, -1);
266			e = 1.0;
267			} else {
268			if (cp == '=') { / absolute setting */
269			p = NULL;
270			e = 1.0/exposure;
271			for (cp++; isspace(*cp); cp++)
272			;
273			if (cp == '\0') { / interactive */
274			sprintf(buf, "exposure (%lf): ", exposure);
275			(*dev->comout)(buf);
276	greg	1.18	(*dev->comin)(buf, NULL);
277	greg	1.1	for (cp = buf; isspace(*cp); cp++)
278			;
279			if (*cp == '\0')
280			return;
281			}
282			} else { /* normalize to average */
283			p = &ptrunk;
284			e = 1.0;
285			}
286			if (cp == '+' \|\| cp == '-') /* f-stops */
287			e *= pow(2.0, atof(cp));
288			else /* multiplier */
289			e *= atof(cp);
290			}
291			if (p != NULL) { /* relative setting */
292	greg	1.6	if (bright(p->v) <= FTINY) {
293	greg	1.1	error(COMMAND, "cannot normalize to zero");
294			return;
295			}
296	greg	1.6	e *= 0.5 / bright(p->v);
297	greg	1.1	}
298			if (e <= FTINY \|\| fabs(1.0 - e) <= FTINY)
299			return;
300			scalepict(&ptrunk, e);
301			exposure *= e;
302			redraw();
303			}
304
305
306	greg	1.24	getparam(str, dsc, typ, ptr) /* get variable from user */
307			char str, dsc;
308			int typ;
309			register union {int i; double d; COLOR C;} *ptr;
310			{
311	greg	1.25	extern char *index();
312	greg	1.24	int i0;
313			double d0, d1, d2;
314			char buf[48];
315
316			switch (typ) {
317			case 'i': /* integer */
318			if (sscanf(str, "%d", &i0) != 1) {
319			(*dev->comout)(dsc);
320			sprintf(buf, " (%d): ", ptr->i);
321			(*dev->comout)(buf);
322			(*dev->comin)(buf, NULL);
323			if (sscanf(buf, "%d", &i0) != 1)
324			break;
325			}
326			ptr->i = i0;
327			break;
328			case 'r': /* real */
329			if (sscanf(str, "%lf", &d0) != 1) {
330			(*dev->comout)(dsc);
331			sprintf(buf, " (%.6g): ", ptr->d);
332			(*dev->comout)(buf);
333			(*dev->comin)(buf, NULL);
334			if (sscanf(buf, "%lf", &d0) != 1)
335			break;
336			}
337			ptr->d = d0;
338			break;
339			case 'b': /* boolean */
340			if (sscanf(str, "%1s", buf) != 1) {
341			(*dev->comout)(dsc);
342			sprintf(buf, " (%c): ", ptr->i ? 'y' : 'n');
343			(*dev->comout)(buf);
344			(*dev->comin)(buf, NULL);
345	greg	1.25	if (buf[0] == '\0' \|\|
346			index("yY+1tTnN-0fF", buf[0]) == NULL)
347			break;
348	greg	1.24	}
349	greg	1.25	ptr->i = index("yY+1tT", buf[0]) != NULL;
350	greg	1.24	break;
351			case 'C': /* color */
352			if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
353			(*dev->comout)(dsc);
354			sprintf(buf, " (%.6g %.6g %.6g): ",
355			colval(ptr->C,RED),
356			colval(ptr->C,GRN),
357			colval(ptr->C,BLU));
358			(*dev->comout)(buf);
359			(*dev->comin)(buf, NULL);
360			if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
361			break;
362			}
363			setcolor(ptr->C, d0, d1, d2);
364			break;
365			}
366			}
367
368
369	greg	1.1	setparam(s) /* get/set program parameter */
370			register char *s;
371			{
372			extern int psample;
373			extern double maxdiff;
374			extern double minweight;
375			extern int maxdepth;
376			extern double dstrsrc;
377	greg	1.5	extern double shadthresh;
378	greg	1.9	extern double shadcert;
379	greg	1.1	extern COLOR ambval;
380			extern double ambacc;
381			extern double minarad;
382			extern int ambres;
383			extern int ambdiv;
384			extern int ambssamp;
385			extern int ambounce;
386	greg	1.23	extern int directinvis;
387			extern int do_irrad;
388	greg	1.24	char buf[128];
389	greg	1.1
390			if (s[0] == '\0') {
391	greg	1.23	(*dev->comout)(
392	greg	1.24	"aa ab ad ar as av b dc di dj dt i lr lw sp st: ");
393	greg	1.18	(*dev->comin)(buf, NULL);
394	greg	1.1	s = buf;
395			}
396			switch (s[0]) {
397			case 'l': /* limit */
398			switch (s[1]) {
399			case 'w': /* weight */
400	greg	1.24	getparam(s+2, "limit weight", 'r', &minweight);
401	greg	1.1	break;
402			case 'r': /* reflection */
403	greg	1.24	getparam(s+2, "limit reflection", 'i', &maxdepth);
404	greg	1.1	break;
405			default:
406			goto badparam;
407			}
408			break;
409	greg	1.5	case 'd': /* direct */
410			switch (s[1]) {
411			case 'j': /* jitter */
412	greg	1.24	getparam(s+2, "direct jitter", 'r', &dstrsrc);
413	greg	1.9	break;
414			case 'c': /* certainty */
415	greg	1.24	getparam(s+2, "direct certainty", 'r', &shadcert);
416	greg	1.5	break;
417	greg	1.8	case 't': /* threshold */
418	greg	1.24	getparam(s+2, "direct threshold", 'r', &shadthresh);
419	greg	1.5	break;
420	greg	1.23	case 'i': /* invisibility */
421	greg	1.24	getparam(s+2, "direct invisibility",
422			'b', &directinvis);
423	greg	1.23	break;
424	greg	1.5	default:
425	greg	1.1	goto badparam;
426			}
427	greg	1.23	break;
428			case 'b': /* black and white */
429	greg	1.24	getparam(s+1, "black and white", 'b', &greyscale);
430	greg	1.23	break;
431			case 'i': /* irradiance */
432	greg	1.24	getparam(s+1, "irradiance", 'b', &do_irrad);
433	greg	1.1	break;
434			case 'a': /* ambient */
435			switch (s[1]) {
436			case 'v': /* value */
437	greg	1.24	getparam(s+2, "ambient value", 'C', ambval);
438	greg	1.1	break;
439			case 'a': /* accuracy */
440	greg	1.24	getparam(s+2, "ambient accuracy", 'r', &ambacc);
441	greg	1.1	break;
442			case 'd': /* divisions */
443	greg	1.24	getparam(s+2, "ambient divisions", 'i', &ambdiv);
444	greg	1.1	break;
445			case 's': /* samples */
446	greg	1.24	getparam(s+2, "ambient super-samples", 'i', &ambssamp);
447	greg	1.1	break;
448			case 'b': /* bounces */
449	greg	1.24	getparam(s+2, "ambient bounces", 'i', &ambounce);
450	greg	1.1	break;
451			case 'r':
452	greg	1.24	getparam(s+2, "ambient resolution", 'i', &ambres);
453	greg	1.1	minarad = ambres > 0 ? thescene.cusize/ambres : 0.0;
454			break;
455			default:
456			goto badparam;
457			}
458			break;
459			case 's': /* sample */
460			switch (s[1]) {
461			case 'p': /* pixel */
462	greg	1.24	getparam(s+2, "sample pixel", 'i', &psample);
463	greg	1.1	pdepth = 0;
464			break;
465	greg	1.8	case 't': /* threshold */
466	greg	1.24	getparam(s+2, "sample threshold", 'r', &maxdiff);
467	greg	1.1	pdepth = 0;
468			break;
469			default:
470			goto badparam;
471			}
472			break;
473			case '\0': /* nothing */
474			break;
475			default:;
476			badparam:
477	greg	1.26	*sskip(s) = '\0';
478	greg	1.1	sprintf(errmsg, "%s: unknown variable", s);
479			error(COMMAND, errmsg);
480			break;
481			}
482			}
483
484
485			traceray(s) /* trace a single ray */
486			char *s;
487			{
488			char buf[128];
489			int x, y;
490			RAY thisray;
491
492			if (sscanf(s, "%lf %lf %lf %lf %lf %lf",
493			&thisray.rorg[0], &thisray.rorg[1], &thisray.rorg[2],
494			&thisray.rdir[0], &thisray.rdir[1], &thisray.rdir[2]) != 6) {
495
496			if (dev->getcur == NULL)
497			return;
498			(*dev->comout)("Pick ray\n");
499			if ((*dev->getcur)(&x, &y) == ABORT)
500			return;
501
502	greg	1.21	if (viewray(thisray.rorg, thisray.rdir, &ourview,
503			(x+.5)/hresolu, (y+.5)/vresolu) < 0) {
504			error(COMMAND, "not on image");
505			return;
506			}
507
508	greg	1.1	} else if (normalize(thisray.rdir) == 0.0) {
509			error(COMMAND, "zero ray direction");
510			return;
511			}
512
513			rayorigin(&thisray, NULL, PRIMARY, 1.0);
514
515			rayvalue(&thisray);
516
517			if (thisray.ro == NULL)
518			(*dev->comout)("ray hit nothing");
519			else {
520			sprintf(buf, "ray hit %s %s \"%s\"",
521			objptr(thisray.ro->omod)->oname,
522			ofun[thisray.ro->otype].funame,
523			thisray.ro->oname);
524			(*dev->comout)(buf);
525	greg	1.18	(*dev->comin)(buf, NULL);
526	greg	1.1	if (thisray.rot >= FHUGE)
527			(*dev->comout)("at infinity");
528			else {
529			sprintf(buf, "at (%.6g %.6g %.6g)", thisray.rop[0],
530			thisray.rop[1], thisray.rop[2]);
531			(*dev->comout)(buf);
532			}
533	greg	1.18	(*dev->comin)(buf, NULL);
534	greg	1.1	sprintf(buf, "with value (%.6g %.6g %.6g)",
535			colval(thisray.rcol,RED),
536			colval(thisray.rcol,GRN),
537			colval(thisray.rcol,BLU));
538			(*dev->comout)(buf);
539			}
540	greg	1.18	(*dev->comin)(buf, NULL);
541	greg	1.1	}
542
543
544			writepict(s) /* write the picture to a file */
545			char *s;
546			{
547			static char buf[128];
548	greg	1.10	char *fname;
549	greg	1.1	FILE *fp;
550			COLR *scanline;
551			int y;
552
553			if (sscanf(s, "%s", buf) != 1 && buf[0] == '\0') {
554			error(COMMAND, "no file");
555			return;
556			}
557	greg	1.10	if ((fname = getpath(buf, NULL, 0)) == NULL \|\|
558	greg	1.1	(fp = fopen(fname, "w")) == NULL) {
559			sprintf(errmsg, "cannot open \"%s\"", buf);
560			error(COMMAND, errmsg);
561			return;
562			}
563			(*dev->comout)("writing \"");
564			(*dev->comout)(fname);
565			(*dev->comout)("\"...\n");
566			/* write header */
567			fputs(progname, fp);
568			fprintview(&ourview, fp);
569	greg	1.29	if (octname != NULL)
570			fprintf(fp, " %s\n", octname);
571			else
572			putc('\n', fp);
573	greg	1.28	fprintf(fp, "SOFTWARE= %s\n", VersionID);
574	greg	1.1	if (exposure != 1.0)
575	greg	1.13	fputexpos(exposure, fp);
576	greg	1.15	if (dev->pixaspect != 1.0)
577			fputaspect(dev->pixaspect, fp);
578	greg	1.22	fputformat(COLRFMT, fp);
579	greg	1.11	putc('\n', fp);
580	greg	1.14	fputresolu(YMAJOR\|YDECR, hresolu, vresolu, fp);
581	greg	1.1
582	greg	1.14	scanline = (COLR )malloc(hresolusizeof(COLR));
583	greg	1.27	if (scanline == NULL) {
584			error(COMMAND, "not enough memory!");
585			fclose(fp);
586			unlink(fname);
587			return;
588			}
589	greg	1.14	for (y = vresolu-1; y >= 0; y--) {
590			getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu);
591			if (fwritecolrs(scanline, hresolu, fp) < 0)
592	greg	1.1	break;
593			}
594	greg	1.27	free((char *)scanline);
595	greg	1.1	if (fclose(fp) < 0)
596			error(COMMAND, "write error");
597			}