src/rt/rv2.c

/* Copyright (c) 1993 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  "resolu.h"

#include  <ctype.h>

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

#ifdef  SMLFLT
#define  sscanvec(s,v)  (sscanf(s,"%f %f %f",v,v+1,v+2)==3)
#else
#define  sscanvec(s,v)  (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3)
#endif


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;

        while (isspace(*s))
                s++;
        if (*s == '-') {                        /* command line parameters */
                copystruct(&nv, &ourview);
                if (sscanview(&nv, s))
                        newview(&nv);
                else
                        error(COMMAND, "bad view option(s)");
                return;
        }
        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 (sscanvec(buf, nv.vp))
                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 (sscanvec(buf, nv.vdir))
                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 (sscanvec(buf, nv.vup))
                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, NULL)) == -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;
{
        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;
{
        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;
{
        char  buf[128];
        register char  *cp;
        RECT  r;
        int  x, y;
        register PNODE  *p = &ptrunk;
        int  adapt = 0;
        double  e = 1.0;
start:
        for (cp = s; isspace(*cp); cp++)
                ;
        if (*cp == '@') {
                adapt++;
                goto start;
        }
        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);
        } else {
                if (*cp == '=') {       /* absolute setting */
                        p = NULL;
                        e = 1.0/exposure;
                        for (cp++; isspace(*cp); cp++)
                                ;
                        if (*cp == '\0') {      /* interactive */
                                sprintf(buf, "exposure (%f): ", exposure);
                                (*dev->comout)(buf);
                                (*dev->comin)(buf, NULL);
                                for (cp = buf; isspace(*cp); cp++)
                                        ;
                                if (*cp == '\0')
                                        return;
                        }
                }
                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) < 1e-15) {
                        error(COMMAND, "cannot normalize to zero");
                        return;
                }
                if (adapt)
                        e *= 106./pow(1.219+pow(luminance(p->v)/exposure,.4),2.5)/exposure;
                else
                        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)
                                return(0);
                }
                ptr->i = i0;
                return(1);
        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)
                                return(0);
                }
                ptr->d = d0;
                return(1);
        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)
                                return(0);
                }
                ptr->i = index("yY+1tT", buf[0]) != NULL;
                return(1);
        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)
                                return(0);
                }
                setcolor(ptr->C, d0, d1, d2);
                return(1);
        }
}


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 int  ambres;
        extern int  ambdiv;
        extern int  ambssamp;
        extern int  ambounce;
        extern int  directvis;
        extern double  srcsizerat;
        extern int  do_irrad;
        extern double  specjitter;
        extern double  specthresh;
        char  buf[128];
        
        if (s[0] == '\0') {
                (*dev->comout)(
                "aa ab ad ar as av b dc di dj ds dt i lr lw ps pt sj 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 'v':                       /* visibility */
                        getparam(s+2, "direct visibility",
                                        'b', &directvis);
                        break;
                case 's':                       /* sampling */
                        getparam(s+2, "direct sampling", 'r', &srcsizerat);
                        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', (COLOR *)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':
                        if (getparam(s+2, "ambient resolution", 'i', &ambres))
                                setambres(ambres);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'p':                       /* pixel */
                switch (s[1]) {
                case 's':                       /* sample */
                        if (getparam(s+2, "pixel sample", 'i', &psample))
                                pdepth = 0;
                        break;
                case 't':                       /* threshold */
                        if (getparam(s+2, "pixel threshold", 'r', &maxdiff))
                                pdepth = 0;
                        break;
                default:
                        goto badparam;
                }
                break;
        case 's':                       /* specular */
                switch (s[1]) {
                case 'j':                       /* jitter */
                        getparam(s+2, "specular jitter", 'r', &specjitter);
                        break;
                case 't':                       /* threshold */
                        getparam(s+2, "specular threshold", 'r', &specthresh);
                        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 (!sscanvec(s, thisray.rorg) ||
                        !sscanvec(sskip(sskip(sskip(s))), thisray.rdir)) {

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