src/rt/rv2.c

#ifndef lint
static const char       RCSid[] = "$Id: rv2.c,v 2.40 2003/06/05 19:29:34 schorsch Exp $";
#endif
/*
 *  rv2.c - command routines used in tracing a view.
 *
 *  External symbols declared in rpaint.h
 */

#include "copyright.h"

#include  <ctype.h>
#include  <string.h>

#include  "platform.h"
#include  "ray.h"
#include  "otypes.h"
#include  "rpaint.h"

extern int  psample;                    /* pixel sample size */
extern double  maxdiff;                 /* max. sample difference */

#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  rifname[128];              /* rad input file name */

extern char  *progname;
extern char  *octname;


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


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

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


void
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);
                putc('\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, "fore and aft clipping plane (%.6g %.6g): ",
                        ourview.vfore, ourview.vaft);
        (*dev->comout)(buf);
        (*dev->comin)(buf, NULL);
        if (buf[0] == CTRL('C')) return;
        if (sscanf(buf, "%lf %lf", &nv.vfore, &nv.vaft) == 2)
                change++;
        else {
                nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
        }
        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);
}


void
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, "", R_OK)) == 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();
}


void
saveview(s)                             /* save view to rad file */
char  *s;
{
        char  view[64];
        char  *fname;
        FILE  *fp;

        if (*atos(view, sizeof(view), s)) {
                if (isint(view)) {
                        error(COMMAND, "cannot write view by number");
                        return;
                }
                s = sskip(s);
        }
        while (isspace(*s))
                s++;
        if (*s)
                atos(rifname, sizeof(rifname), s);
        else if (rifname[0] == '\0') {
                error(COMMAND, "no previous rad file");
                return;
        }
        if ((fname = getpath(rifname, NULL, 0)) == NULL ||
                        (fp = fopen(fname, "a")) == NULL) {
                sprintf(errmsg, "cannot open \"%s\"", rifname);
                error(COMMAND, errmsg);
                return;
        }
        fputs("view= ", fp);
        fputs(view, fp);
        fprintview(&ourview, fp);
        putc('\n', fp);
        fclose(fp);
}


void
loadview(s)                             /* load view from rad file */
char  *s;
{
        char  buf[512];
        char  *fname;
        FILE  *fp;
        VIEW  nv;

        strcpy(buf, "rad -n -s -V -v ");
        if (sscanf(s, "%s", buf+strlen(buf)) == 1)
                s = sskip(s);
        else
                strcat(buf, "1");
        if (*s)
                atos(rifname, sizeof(rifname), s);
        else if (rifname[0] == '\0') {
                error(COMMAND, "no previous rad file");
                return;
        }
        if ((fname = getpath(rifname, "", R_OK)) == NULL) {
                sprintf(errmsg, "cannot access \"%s\"", rifname);
                error(COMMAND, errmsg);
                return;
        }
        sprintf(buf+strlen(buf), " %s", fname);
        if ((fp = popen(buf, "r")) == NULL) {
                error(COMMAND, "cannot run rad");
                return;
        }
        buf[0] = '\0';
        fgets(buf, sizeof(buf), fp);
        pclose(fp);
        copystruct(&nv, &stdview);
        if (!sscanview(&nv, buf)) {
                error(COMMAND, "rad error -- no such view?");
                return;
        }
        newview(&nv);
}


void
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.horiz = ourview.horiz; nv.vert = ourview.vert;
        nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
        nv.hoff = ourview.hoff; nv.voff = ourview.voff;
        zoomview(&nv, zfact);
        newview(&nv);
}


void
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);
}


void
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;
        nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
        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);
}


void
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(sskip2(s,2), 0, vc, &mag) < 0)
                return;
        moveview(angle, elev, mag, vc);
}


void
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;

        for (cp = s; isspace(*cp); cp++)
                ;
        if (*cp == '@') {
                adapt++;
                while (isspace(*++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);
        } 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();
}

typedef union {int i; double d; COLOR C;}       *MyUptr;

int
getparam(str, dsc, typ, p)              /* get variable from user */
char  *str, *dsc;
int  typ;
void  *p;
{
        register MyUptr  ptr = (MyUptr)p;
        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' ||
                                        strchr("yY+1tTnN-0fF", buf[0]) == NULL)
                                return(0);
                }
                ptr->i = strchr("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);
        }
}


void
setparam(s)                             /* get/set program parameter */
register char  *s;
{
        char  buf[128];
        
        if (s[0] == '\0') {
                (*dev->comout)(
                "aa ab ad ar as av aw b dc dv dj ds dt i lr lw me ma mg ms ps pt sj st bv: ");
                (*dev->comin)(buf, NULL);
                s = buf;
        }
        switch (s[0]) {
        case 'l':                       /* limit */
                switch (s[1]) {
                case 'w':                       /* weight */
                        getparam(s+2, "limit weight", 'r',
                                        (void *)&minweight);
                        break;
                case 'r':                       /* reflection */
                        getparam(s+2, "limit reflection", 'i',
                                        (void *)&maxdepth);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'd':                       /* direct */
                switch (s[1]) {
                case 'j':                       /* jitter */
                        getparam(s+2, "direct jitter", 'r',
                                        (void *)&dstrsrc);
                        break;
                case 'c':                       /* certainty */
                        getparam(s+2, "direct certainty", 'r',
                                        (void *)&shadcert);
                        break;
                case 't':                       /* threshold */
                        getparam(s+2, "direct threshold", 'r',
                                        (void *)&shadthresh);
                        break;
                case 'v':                       /* visibility */
                        getparam(s+2, "direct visibility", 'b',
                                        (void *)&directvis);
                        break;
                case 's':                       /* sampling */
                        getparam(s+2, "direct sampling", 'r',
                                        (void *)&srcsizerat);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'b':                       /* back faces or black and white */
                switch (s[1]) {
                case 'v':                       /* back face visibility */
                        getparam(s+2, "back face visibility", 'b',
                                        (void *)&backvis);
                        break;
                case '\0':                      /* black and white */
                case ' ':
                case 'y': case 'Y': case 't': case 'T': case '1': case '+':
                case 'n': case 'N': case 'f': case 'F': case '0': case '-':
                        getparam(s+1, "black and white", 'b',
                                        (void *)&greyscale);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'i':                       /* irradiance */
                getparam(s+1, "irradiance", 'b',
                                (void *)&do_irrad);
                break;
        case 'a':                       /* ambient */
                switch (s[1]) {
                case 'v':                       /* value */
                        getparam(s+2, "ambient value", 'C',
                                        (void *)ambval);
                        break;
                case 'w':                       /* weight */
                        getparam(s+2, "ambient value weight", 'i',
                                        (void *)&ambvwt);
                        break;
                case 'a':                       /* accuracy */
                        if (getparam(s+2, "ambient accuracy", 'r',
                                        (void *)&ambacc))
                                setambacc(ambacc);
                        break;
                case 'd':                       /* divisions */
                        getparam(s+2, "ambient divisions", 'i',
                                        (void *)&ambdiv);
                        break;
                case 's':                       /* samples */
                        getparam(s+2, "ambient super-samples", 'i',
                                        (void *)&ambssamp);
                        break;
                case 'b':                       /* bounces */
                        getparam(s+2, "ambient bounces", 'i',
                                        (void *)&ambounce);
                        break;
                case 'r':
                        if (getparam(s+2, "ambient resolution", 'i',
                                        (void *)&ambres))
                                setambres(ambres);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'm':                       /* medium */
                switch (s[1]) {
                case 'e':                       /* extinction coefficient */
                        getparam(s+2, "extinction coefficient", 'C',
                                        (void *)cextinction);
                        break;
                case 'a':                       /* scattering albedo */
                        getparam(s+2, "scattering albedo", 'C',
                                        (void *)salbedo);
                        break;
                case 'g':                       /* scattering eccentricity */
                        getparam(s+2, "scattering eccentricity", 'r',
                                        (void *)&seccg);
                        break;
                case 's':                       /* sampling distance */
                        getparam(s+2, "mist sampling distance", 'r',
                                        (void *)&ssampdist);
                        break;
                default:
                        goto badparam;
                }
                break;
        case 'p':                       /* pixel */
                switch (s[1]) {
                case 's':                       /* sample */
                        if (getparam(s+2, "pixel sample", 'i',
                                        (void *)&psample))
                                pdepth = 0;
                        break;
                case 't':                       /* threshold */
                        if (getparam(s+2, "pixel threshold", 'r',
                                        (void *)&maxdiff))
                                pdepth = 0;
                        break;
                default:
                        goto badparam;
                }
                break;
        case 's':                       /* specular */
                switch (s[1]) {
                case 'j':                       /* jitter */
                        getparam(s+2, "specular jitter", 'r',
                                        (void *)&specjitter);
                        break;
                case 't':                       /* threshold */
                        getparam(s+2, "specular threshold", 'r',
                                        (void *)&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;
        }
}


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

        thisray.rmax = 0.0;

        if (!sscanvec(s, thisray.rorg) ||
                        !sscanvec(sskip2(s,3), thisray.rdir)) {

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

                if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir,
                        &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) {
                        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 \"%s\"",
                                thisray.rod < 0.0 ? "back of " : "",
                                thisray.ro->omod == OVOID ? VOIDID :
                                        objptr(thisray.ro->omod)->oname,
                                ofun[thisray.ro->otype].funame,
                                thisray.ro->oname);
                if ((ino = objptr(thisray.robj)) != thisray.ro)
                        sprintf(buf+strlen(buf), " in %s \"%s\"",
                                        ofun[ino->otype].funame, ino->oname);
                (*dev->comout)(buf);
                (*dev->comin)(buf, NULL);
                if (thisray.rot >= FHUGE)
                        (*dev->comout)("at infinity");
                else {
                        sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)",
                                        thisray.rop[0], thisray.rop[1],
                                        thisray.rop[2], thisray.rt);
                        (*dev->comout)(buf);
                }
                (*dev->comin)(buf, NULL);
                sprintf(buf, "value (%.5g %.5g %.5g) (%.3gL)",
                                colval(thisray.rcol,RED),
                                colval(thisray.rcol,GRN),
                                colval(thisray.rcol,BLU),
                                luminance(thisray.rcol));
                (*dev->comout)(buf);
        }
        (*dev->comin)(buf, NULL);
}


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

        while (isspace(*s))
                s++;
        if (*s)
                atos(buf, sizeof(buf), s);
        else if (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;
        }
        SET_FILE_BINARY(fp);
        (*dev->comout)("writing \"");
        (*dev->comout)(fname);
        (*dev->comout)("\"...\n");
                                                /* write header */
        newheader("RADIANCE", fp);
        fputs(progname, fp);
        fprintview(&ourview, fp);
        if (octname != NULL)
                fprintf(fp, " %s\n", octname);
        else
                putc('\n', fp);
        fprintf(fp, "SOFTWARE= %s\n", VersionID);
        fputnow(fp);
        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((void *)scanline);
        if (fclose(fp) < 0)
                error(COMMAND, "write error");
}
Revision:	2.41
Committed:	Mon Jun 30 14:59:13 2003 UTC (20 years, 10 months ago) by schorsch
Content type:	text/plain
Branch:	MAIN
Changes since 2.40:	+6 -5 lines
Log Message:	Replaced most outdated BSD function calls with their posix equivalents, and cleaned up a few other platform dependencies.
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	schorsch	2.41	static const char RCSid[] = "$Id: rv2.c,v 2.40 2003/06/05 19:29:34 schorsch Exp $";
3	greg	1.1	#endif
4			/*
5			* rv2.c - command routines used in tracing a view.
6			*
7	greg	2.38	* External symbols declared in rpaint.h
8			*/
9
10	schorsch	2.41	#include "copyright.h"
11
12	schorsch	2.40	#include <ctype.h>
13	schorsch	2.41	#include <string.h>
14	greg	1.1
15	schorsch	2.40	#include "platform.h"
16	greg	1.1	#include "ray.h"
17			#include "otypes.h"
18			#include "rpaint.h"
19
20	greg	2.38	extern int psample; /* pixel sample size */
21			extern double maxdiff; /* max. sample difference */
22	greg	1.1
23	greg	2.4	#define CTRL(c) ((c)-'@')
24	greg	1.1
25	greg	2.6	#ifdef SMLFLT
26			#define sscanvec(s,v) (sscanf(s,"%f %f %f",v,v+1,v+2)==3)
27			#else
28			#define sscanvec(s,v) (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3)
29			#endif
30
31	greg	2.20	extern char rifname[128]; /* rad input file name */
32	greg	2.18
33	greg	1.1	extern char *progname;
34	greg	1.29	extern char *octname;
35	greg	1.1
36
37	greg	2.38	void
38	greg	1.1	getframe(s) /* get a new frame */
39			char *s;
40			{
41	greg	1.12	if (getrect(s, &pframe) < 0)
42	greg	1.1	return;
43			pdepth = 0;
44			}
45
46
47	greg	2.38	void
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	2.38	void
60	greg	1.1	getview(s) /* get/show view parameters */
61			char *s;
62			{
63			FILE *fp;
64			char buf[128];
65	greg	1.10	char *fname;
66	greg	1.3	int change = 0;
67	greg	1.1	VIEW nv;
68
69	greg	2.9	while (isspace(*s))
70			s++;
71			if (s == '-') { / command line parameters */
72			copystruct(&nv, &ourview);
73			if (sscanview(&nv, s))
74			newview(&nv);
75			else
76			error(COMMAND, "bad view option(s)");
77			return;
78			}
79	greg	1.1	if (sscanf(s, "%s", buf) == 1) { /* write parameters to a file */
80	greg	1.10	if ((fname = getpath(buf, NULL, 0)) == NULL \|\|
81	greg	1.1	(fp = fopen(fname, "a")) == NULL) {
82			sprintf(errmsg, "cannot open \"%s\"", buf);
83			error(COMMAND, errmsg);
84			return;
85			}
86	greg	1.2	fputs(progname, fp);
87	greg	1.1	fprintview(&ourview, fp);
88	greg	1.19	fputs(sskip(s), fp);
89	greg	2.18	putc('\n', fp);
90	greg	1.1	fclose(fp);
91			return;
92			}
93			sprintf(buf, "view type (%c): ", ourview.type);
94			(*dev->comout)(buf);
95	greg	1.18	(*dev->comin)(buf, NULL);
96	greg	2.4	if (buf[0] == CTRL('C')) return;
97	greg	1.3	if (buf[0] && buf[0] != ourview.type) {
98	greg	1.1	nv.type = buf[0];
99	greg	1.3	change++;
100			} else
101	greg	1.1	nv.type = ourview.type;
102			sprintf(buf, "view point (%.6g %.6g %.6g): ",
103			ourview.vp[0], ourview.vp[1], ourview.vp[2]);
104			(*dev->comout)(buf);
105	greg	1.18	(*dev->comin)(buf, NULL);
106	greg	2.4	if (buf[0] == CTRL('C')) return;
107	greg	2.6	if (sscanvec(buf, nv.vp))
108	greg	1.3	change++;
109			else
110	greg	1.1	VCOPY(nv.vp, ourview.vp);
111			sprintf(buf, "view direction (%.6g %.6g %.6g): ",
112			ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
113			(*dev->comout)(buf);
114	greg	1.18	(*dev->comin)(buf, NULL);
115	greg	2.4	if (buf[0] == CTRL('C')) return;
116	greg	2.6	if (sscanvec(buf, nv.vdir))
117	greg	1.3	change++;
118			else
119	greg	1.1	VCOPY(nv.vdir, ourview.vdir);
120			sprintf(buf, "view up (%.6g %.6g %.6g): ",
121			ourview.vup[0], ourview.vup[1], ourview.vup[2]);
122			(*dev->comout)(buf);
123	greg	1.18	(*dev->comin)(buf, NULL);
124	greg	2.4	if (buf[0] == CTRL('C')) return;
125	greg	2.6	if (sscanvec(buf, nv.vup))
126	greg	1.3	change++;
127			else
128	greg	1.1	VCOPY(nv.vup, ourview.vup);
129			sprintf(buf, "view horiz and vert size (%.6g %.6g): ",
130			ourview.horiz, ourview.vert);
131			(*dev->comout)(buf);
132	greg	1.18	(*dev->comin)(buf, NULL);
133	greg	2.4	if (buf[0] == CTRL('C')) return;
134	greg	1.3	if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2)
135			change++;
136			else {
137	greg	1.1	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
138			}
139	greg	2.26	sprintf(buf, "fore and aft clipping plane (%.6g %.6g): ",
140			ourview.vfore, ourview.vaft);
141			(*dev->comout)(buf);
142			(*dev->comin)(buf, NULL);
143			if (buf[0] == CTRL('C')) return;
144			if (sscanf(buf, "%lf %lf", &nv.vfore, &nv.vaft) == 2)
145			change++;
146			else {
147			nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
148			}
149	greg	1.14	sprintf(buf, "view shift and lift (%.6g %.6g): ",
150			ourview.hoff, ourview.voff);
151	greg	1.1	(*dev->comout)(buf);
152	greg	1.18	(*dev->comin)(buf, NULL);
153	greg	2.4	if (buf[0] == CTRL('C')) return;
154	greg	1.14	if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2)
155	greg	1.3	change++;
156			else {
157	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
158	greg	1.1	}
159	greg	1.3	if (change)
160			newview(&nv);
161	greg	1.1	}
162
163
164	greg	2.38	void
165	greg	1.1	lastview(s) /* return to a previous view */
166			char *s;
167			{
168			char buf[128];
169	greg	1.10	char *fname;
170	greg	1.1	int success;
171			VIEW nv;
172
173			if (sscanf(s, "%s", buf) == 1) { /* get parameters from a file */
174	greg	1.17	copystruct(&nv, &stdview);
175	greg	2.19	if ((fname = getpath(buf, "", R_OK)) == NULL \|\|
176	greg	1.33	(success = viewfile(fname, &nv, NULL)) == -1) {
177	greg	1.1	sprintf(errmsg, "cannot open \"%s\"", buf);
178			error(COMMAND, errmsg);
179			return;
180			}
181			if (!success)
182			error(COMMAND, "wrong file format");
183			else
184			newview(&nv);
185			return;
186			}
187	greg	1.20	if (oldview.type == 0) { /* no old view! */
188	greg	1.1	error(COMMAND, "no previous view");
189			return;
190			}
191	greg	1.17	copystruct(&nv, &ourview);
192			copystruct(&ourview, &oldview);
193			copystruct(&oldview, &nv);
194	greg	1.1	newimage();
195	greg	2.18	}
196
197
198	greg	2.38	void
199	greg	2.18	saveview(s) /* save view to rad file */
200			char *s;
201			{
202			char view[64];
203			char *fname;
204			FILE *fp;
205
206	greg	2.19	if (*atos(view, sizeof(view), s)) {
207			if (isint(view)) {
208			error(COMMAND, "cannot write view by number");
209			return;
210			}
211	greg	2.18	s = sskip(s);
212	greg	2.19	}
213	greg	2.20	while (isspace(*s))
214			s++;
215			if (*s)
216			atos(rifname, sizeof(rifname), s);
217			else if (rifname[0] == '\0') {
218	greg	2.18	error(COMMAND, "no previous rad file");
219			return;
220			}
221			if ((fname = getpath(rifname, NULL, 0)) == NULL \|\|
222			(fp = fopen(fname, "a")) == NULL) {
223			sprintf(errmsg, "cannot open \"%s\"", rifname);
224			error(COMMAND, errmsg);
225			return;
226			}
227			fputs("view= ", fp);
228			fputs(view, fp);
229			fprintview(&ourview, fp);
230			putc('\n', fp);
231			fclose(fp);
232			}
233
234
235	greg	2.38	void
236	greg	2.18	loadview(s) /* load view from rad file */
237			char *s;
238			{
239			char buf[512];
240			char *fname;
241			FILE *fp;
242			VIEW nv;
243
244			strcpy(buf, "rad -n -s -V -v ");
245			if (sscanf(s, "%s", buf+strlen(buf)) == 1)
246			s = sskip(s);
247			else
248			strcat(buf, "1");
249	greg	2.20	if (*s)
250			atos(rifname, sizeof(rifname), s);
251			else if (rifname[0] == '\0') {
252	greg	2.18	error(COMMAND, "no previous rad file");
253			return;
254			}
255	greg	2.19	if ((fname = getpath(rifname, "", R_OK)) == NULL) {
256	greg	2.18	sprintf(errmsg, "cannot access \"%s\"", rifname);
257			error(COMMAND, errmsg);
258			return;
259			}
260			sprintf(buf+strlen(buf), " %s", fname);
261			if ((fp = popen(buf, "r")) == NULL) {
262			error(COMMAND, "cannot run rad");
263			return;
264			}
265			buf[0] = '\0';
266			fgets(buf, sizeof(buf), fp);
267			pclose(fp);
268			copystruct(&nv, &stdview);
269			if (!sscanview(&nv, buf)) {
270			error(COMMAND, "rad error -- no such view?");
271			return;
272			}
273			newview(&nv);
274	greg	1.1	}
275
276
277	greg	2.38	void
278	greg	1.1	getaim(s) /* aim camera */
279			char *s;
280			{
281			double zfact;
282			VIEW nv;
283
284			if (getinterest(s, 1, nv.vdir, &zfact) < 0)
285			return;
286	greg	1.21	nv.type = ourview.type;
287	greg	1.1	VCOPY(nv.vp, ourview.vp);
288			VCOPY(nv.vup, ourview.vup);
289	greg	2.28	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
290			nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
291	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
292	greg	1.21	zoomview(&nv, zfact);
293	greg	1.1	newview(&nv);
294			}
295
296
297	greg	2.38	void
298	greg	1.1	getmove(s) /* move camera */
299			char *s;
300			{
301			FVECT vc;
302			double mag;
303
304			if (getinterest(s, 0, vc, &mag) < 0)
305			return;
306	greg	1.4	moveview(0.0, 0.0, mag, vc);
307	greg	1.1	}
308
309
310	greg	2.38	void
311	greg	1.1	getrotate(s) /* rotate camera */
312			char *s;
313			{
314			VIEW nv;
315			FVECT v1;
316	greg	1.4	double angle, elev, zfact;
317	greg	1.1
318	greg	1.4	elev = 0.0; zfact = 1.0;
319			if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) {
320	greg	1.1	error(COMMAND, "missing angle");
321			return;
322			}
323	greg	1.21	nv.type = ourview.type;
324	greg	1.1	VCOPY(nv.vp, ourview.vp);
325			VCOPY(nv.vup, ourview.vup);
326	greg	1.14	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
327	greg	2.28	nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
328	greg	1.1	spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
329			if (elev != 0.0) {
330			fcross(v1, nv.vdir, ourview.vup);
331			normalize(v1);
332			spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
333			}
334	greg	1.21	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
335			zoomview(&nv, zfact);
336	greg	1.1	newview(&nv);
337			}
338
339
340	greg	2.38	void
341	greg	1.1	getpivot(s) /* pivot viewpoint */
342			register char *s;
343			{
344			FVECT vc;
345	greg	1.4	double angle, elev, mag;
346	greg	1.1
347	greg	1.4	elev = 0.0;
348			if (sscanf(s, "%lf %lf", &angle, &elev) < 1) {
349	greg	1.1	error(COMMAND, "missing angle");
350			return;
351			}
352	greg	2.33	if (getinterest(sskip2(s,2), 0, vc, &mag) < 0)
353	greg	1.1	return;
354	greg	1.4	moveview(angle, elev, mag, vc);
355	greg	1.1	}
356
357
358	greg	2.38	void
359	greg	1.1	getexposure(s) /* get new exposure */
360			char *s;
361			{
362			char buf[128];
363			register char *cp;
364			RECT r;
365			int x, y;
366	greg	2.16	register PNODE *p = &ptrunk;
367			int adapt = 0;
368			double e = 1.0;
369	greg	2.25
370	greg	1.1	for (cp = s; isspace(*cp); cp++)
371			;
372	greg	2.16	if (*cp == '@') {
373			adapt++;
374	greg	2.25	while (isspace(*++cp))
375			;
376	greg	2.16	}
377	greg	1.1	if (cp == '\0') { / normalize to point */
378			if (dev->getcur == NULL)
379			return;
380			(*dev->comout)("Pick point for exposure\n");
381			if ((*dev->getcur)(&x, &y) == ABORT)
382			return;
383			r.l = r.d = 0;
384	greg	1.14	r.r = hresolu; r.u = vresolu;
385	greg	1.1	p = findrect(x, y, &ptrunk, &r, -1);
386			} else {
387			if (cp == '=') { / absolute setting */
388			p = NULL;
389			e = 1.0/exposure;
390			for (cp++; isspace(*cp); cp++)
391			;
392			if (cp == '\0') { / interactive */
393	greg	2.15	sprintf(buf, "exposure (%f): ", exposure);
394	greg	1.1	(*dev->comout)(buf);
395	greg	1.18	(*dev->comin)(buf, NULL);
396	greg	1.1	for (cp = buf; isspace(*cp); cp++)
397			;
398			if (*cp == '\0')
399			return;
400			}
401			}
402			if (cp == '+' \|\| cp == '-') /* f-stops */
403			e *= pow(2.0, atof(cp));
404			else /* multiplier */
405			e *= atof(cp);
406			}
407			if (p != NULL) { /* relative setting */
408	greg	2.14	if (bright(p->v) < 1e-15) {
409	greg	1.1	error(COMMAND, "cannot normalize to zero");
410			return;
411			}
412	greg	2.16	if (adapt)
413			e *= 106./pow(1.219+pow(luminance(p->v)/exposure,.4),2.5)/exposure;
414			else
415			e *= 0.5 / bright(p->v);
416	greg	1.1	}
417			if (e <= FTINY \|\| fabs(1.0 - e) <= FTINY)
418			return;
419			scalepict(&ptrunk, e);
420			exposure *= e;
421			redraw();
422			}
423
424	greg	2.38	typedef union {int i; double d; COLOR C;} *MyUptr;
425	greg	1.1
426	greg	2.38	int
427			getparam(str, dsc, typ, p) /* get variable from user */
428	greg	1.24	char str, dsc;
429			int typ;
430	greg	2.38	void *p;
431	greg	1.24	{
432	greg	2.38	register MyUptr ptr = (MyUptr)p;
433	greg	1.24	int i0;
434			double d0, d1, d2;
435			char buf[48];
436
437			switch (typ) {
438			case 'i': /* integer */
439			if (sscanf(str, "%d", &i0) != 1) {
440			(*dev->comout)(dsc);
441			sprintf(buf, " (%d): ", ptr->i);
442			(*dev->comout)(buf);
443			(*dev->comin)(buf, NULL);
444			if (sscanf(buf, "%d", &i0) != 1)
445	greg	1.31	return(0);
446	greg	1.24	}
447			ptr->i = i0;
448	greg	1.31	return(1);
449	greg	1.24	case 'r': /* real */
450			if (sscanf(str, "%lf", &d0) != 1) {
451			(*dev->comout)(dsc);
452			sprintf(buf, " (%.6g): ", ptr->d);
453			(*dev->comout)(buf);
454			(*dev->comin)(buf, NULL);
455			if (sscanf(buf, "%lf", &d0) != 1)
456	greg	1.31	return(0);
457	greg	1.24	}
458			ptr->d = d0;
459	greg	1.31	return(1);
460	greg	1.24	case 'b': /* boolean */
461			if (sscanf(str, "%1s", buf) != 1) {
462			(*dev->comout)(dsc);
463	greg	2.2	sprintf(buf, "? (%c): ", ptr->i ? 'y' : 'n');
464	greg	1.24	(*dev->comout)(buf);
465			(*dev->comin)(buf, NULL);
466	greg	1.25	if (buf[0] == '\0' \|\|
467	schorsch	2.41	strchr("yY+1tTnN-0fF", buf[0]) == NULL)
468	greg	1.31	return(0);
469	greg	1.24	}
470	schorsch	2.41	ptr->i = strchr("yY+1tT", buf[0]) != NULL;
471	greg	1.31	return(1);
472	greg	1.24	case 'C': /* color */
473			if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
474			(*dev->comout)(dsc);
475			sprintf(buf, " (%.6g %.6g %.6g): ",
476			colval(ptr->C,RED),
477			colval(ptr->C,GRN),
478			colval(ptr->C,BLU));
479			(*dev->comout)(buf);
480			(*dev->comin)(buf, NULL);
481			if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
482	greg	1.31	return(0);
483	greg	1.24	}
484			setcolor(ptr->C, d0, d1, d2);
485	greg	1.31	return(1);
486	greg	1.24	}
487			}
488
489
490	greg	2.38	void
491	greg	1.1	setparam(s) /* get/set program parameter */
492			register char *s;
493			{
494	greg	1.24	char buf[128];
495	greg	1.1
496			if (s[0] == '\0') {
497	greg	1.23	(*dev->comout)(
498	greg	2.34	"aa ab ad ar as av aw b dc dv dj ds dt i lr lw me ma mg ms ps pt sj st bv: ");
499	greg	1.18	(*dev->comin)(buf, NULL);
500	greg	1.1	s = buf;
501			}
502			switch (s[0]) {
503			case 'l': /* limit */
504			switch (s[1]) {
505			case 'w': /* weight */
506	greg	2.38	getparam(s+2, "limit weight", 'r',
507			(void *)&minweight);
508	greg	1.1	break;
509			case 'r': /* reflection */
510	greg	2.38	getparam(s+2, "limit reflection", 'i',
511			(void *)&maxdepth);
512	greg	1.1	break;
513			default:
514			goto badparam;
515			}
516			break;
517	greg	1.5	case 'd': /* direct */
518			switch (s[1]) {
519			case 'j': /* jitter */
520	greg	2.38	getparam(s+2, "direct jitter", 'r',
521			(void *)&dstrsrc);
522	greg	1.9	break;
523			case 'c': /* certainty */
524	greg	2.38	getparam(s+2, "direct certainty", 'r',
525			(void *)&shadcert);
526	greg	1.5	break;
527	greg	1.8	case 't': /* threshold */
528	greg	2.38	getparam(s+2, "direct threshold", 'r',
529			(void *)&shadthresh);
530	greg	1.5	break;
531	greg	2.13	case 'v': /* visibility */
532	greg	2.38	getparam(s+2, "direct visibility", 'b',
533			(void *)&directvis);
534	greg	1.30	break;
535			case 's': /* sampling */
536	greg	2.38	getparam(s+2, "direct sampling", 'r',
537			(void *)&srcsizerat);
538	greg	1.23	break;
539	greg	1.5	default:
540	greg	1.1	goto badparam;
541			}
542	greg	1.23	break;
543	greg	2.27	case 'b': /* back faces or black and white */
544			switch (s[1]) {
545			case 'v': /* back face visibility */
546	greg	2.38	getparam(s+2, "back face visibility", 'b',
547			(void *)&backvis);
548	greg	2.27	break;
549			case '\0': /* black and white */
550	greg	2.29	case ' ':
551	greg	2.27	case 'y': case 'Y': case 't': case 'T': case '1': case '+':
552			case 'n': case 'N': case 'f': case 'F': case '0': case '-':
553	greg	2.38	getparam(s+1, "black and white", 'b',
554			(void *)&greyscale);
555	greg	2.27	break;
556			default:
557			goto badparam;
558			}
559	greg	1.23	break;
560			case 'i': /* irradiance */
561	greg	2.38	getparam(s+1, "irradiance", 'b',
562			(void *)&do_irrad);
563	greg	1.1	break;
564			case 'a': /* ambient */
565			switch (s[1]) {
566			case 'v': /* value */
567	greg	2.38	getparam(s+2, "ambient value", 'C',
568			(void *)ambval);
569	greg	2.34	break;
570			case 'w': /* weight */
571	greg	2.38	getparam(s+2, "ambient value weight", 'i',
572			(void *)&ambvwt);
573	greg	1.1	break;
574			case 'a': /* accuracy */
575	greg	2.38	if (getparam(s+2, "ambient accuracy", 'r',
576			(void *)&ambacc))
577	greg	2.17	setambacc(ambacc);
578	greg	1.1	break;
579			case 'd': /* divisions */
580	greg	2.38	getparam(s+2, "ambient divisions", 'i',
581			(void *)&ambdiv);
582	greg	1.1	break;
583			case 's': /* samples */
584	greg	2.38	getparam(s+2, "ambient super-samples", 'i',
585			(void *)&ambssamp);
586	greg	1.1	break;
587			case 'b': /* bounces */
588	greg	2.38	getparam(s+2, "ambient bounces", 'i',
589			(void *)&ambounce);
590	greg	1.1	break;
591			case 'r':
592	greg	2.38	if (getparam(s+2, "ambient resolution", 'i',
593			(void *)&ambres))
594	greg	2.8	setambres(ambres);
595	greg	2.32	break;
596			default:
597			goto badparam;
598			}
599			break;
600			case 'm': /* medium */
601			switch (s[1]) {
602			case 'e': /* extinction coefficient */
603			getparam(s+2, "extinction coefficient", 'C',
604	greg	2.38	(void *)cextinction);
605	greg	2.32	break;
606			case 'a': /* scattering albedo */
607	greg	2.35	getparam(s+2, "scattering albedo", 'C',
608	greg	2.38	(void *)salbedo);
609	greg	2.32	break;
610			case 'g': /* scattering eccentricity */
611	greg	2.38	getparam(s+2, "scattering eccentricity", 'r',
612			(void *)&seccg);
613	greg	2.32	break;
614			case 's': /* sampling distance */
615			getparam(s+2, "mist sampling distance", 'r',
616	greg	2.38	(void *)&ssampdist);
617	greg	1.1	break;
618			default:
619			goto badparam;
620			}
621			break;
622	greg	2.5	case 'p': /* pixel */
623	greg	1.1	switch (s[1]) {
624	greg	2.5	case 's': /* sample */
625	greg	2.38	if (getparam(s+2, "pixel sample", 'i',
626			(void *)&psample))
627	greg	1.31	pdepth = 0;
628	greg	1.1	break;
629	greg	1.8	case 't': /* threshold */
630	greg	2.38	if (getparam(s+2, "pixel threshold", 'r',
631			(void *)&maxdiff))
632	greg	1.31	pdepth = 0;
633	greg	2.5	break;
634			default:
635			goto badparam;
636			}
637			break;
638			case 's': /* specular */
639			switch (s[1]) {
640			case 'j': /* jitter */
641	greg	2.38	getparam(s+2, "specular jitter", 'r',
642			(void *)&specjitter);
643	greg	2.5	break;
644			case 't': /* threshold */
645	greg	2.38	getparam(s+2, "specular threshold", 'r',
646			(void *)&specthresh);
647	greg	1.1	break;
648			default:
649			goto badparam;
650			}
651			break;
652			case '\0': /* nothing */
653			break;
654			default:;
655			badparam:
656	greg	1.26	*sskip(s) = '\0';
657	greg	1.1	sprintf(errmsg, "%s: unknown variable", s);
658			error(COMMAND, errmsg);
659			break;
660			}
661			}
662
663
664	greg	2.38	void
665	greg	1.1	traceray(s) /* trace a single ray */
666			char *s;
667			{
668			char buf[128];
669			int x, y;
670	greg	2.37	OBJREC *ino;
671	greg	1.1	RAY thisray;
672
673	greg	2.26	thisray.rmax = 0.0;
674
675	greg	2.6	if (!sscanvec(s, thisray.rorg) \|\|
676	greg	2.33	!sscanvec(sskip2(s,3), thisray.rdir)) {
677	greg	1.1
678			if (dev->getcur == NULL)
679			return;
680			(*dev->comout)("Pick ray\n");
681			if ((*dev->getcur)(&x, &y) == ABORT)
682			return;
683
684	greg	2.26	if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir,
685			&ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) {
686	greg	1.21	error(COMMAND, "not on image");
687			return;
688			}
689
690	greg	1.1	} else if (normalize(thisray.rdir) == 0.0) {
691			error(COMMAND, "zero ray direction");
692			return;
693			}
694
695			rayorigin(&thisray, NULL, PRIMARY, 1.0);
696
697			rayvalue(&thisray);
698
699			if (thisray.ro == NULL)
700			(*dev->comout)("ray hit nothing");
701			else {
702	greg	2.23	sprintf(buf, "ray hit %s%s %s \"%s\"",
703			thisray.rod < 0.0 ? "back of " : "",
704	greg	2.22	thisray.ro->omod == OVOID ? VOIDID :
705			objptr(thisray.ro->omod)->oname,
706	greg	1.1	ofun[thisray.ro->otype].funame,
707			thisray.ro->oname);
708	greg	2.37	if ((ino = objptr(thisray.robj)) != thisray.ro)
709			sprintf(buf+strlen(buf), " in %s \"%s\"",
710			ofun[ino->otype].funame, ino->oname);
711	greg	1.1	(*dev->comout)(buf);
712	greg	1.18	(*dev->comin)(buf, NULL);
713	greg	1.1	if (thisray.rot >= FHUGE)
714			(*dev->comout)("at infinity");
715			else {
716	greg	2.30	sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)",
717			thisray.rop[0], thisray.rop[1],
718			thisray.rop[2], thisray.rt);
719	greg	1.1	(*dev->comout)(buf);
720			}
721	greg	1.18	(*dev->comin)(buf, NULL);
722	greg	2.38	sprintf(buf, "value (%.5g %.5g %.5g) (%.3gL)",
723	greg	1.1	colval(thisray.rcol,RED),
724			colval(thisray.rcol,GRN),
725	greg	1.32	colval(thisray.rcol,BLU),
726			luminance(thisray.rcol));
727	greg	1.1	(*dev->comout)(buf);
728			}
729	greg	1.18	(*dev->comin)(buf, NULL);
730	greg	1.1	}
731
732
733	greg	2.38	void
734	greg	1.1	writepict(s) /* write the picture to a file */
735			char *s;
736			{
737			static char buf[128];
738	greg	1.10	char *fname;
739	greg	1.1	FILE *fp;
740			COLR *scanline;
741			int y;
742
743	greg	2.20	while (isspace(*s))
744			s++;
745			if (*s)
746			atos(buf, sizeof(buf), s);
747			else if (buf[0] == '\0') {
748	greg	1.1	error(COMMAND, "no file");
749			return;
750			}
751	greg	1.10	if ((fname = getpath(buf, NULL, 0)) == NULL \|\|
752	greg	1.1	(fp = fopen(fname, "w")) == NULL) {
753			sprintf(errmsg, "cannot open \"%s\"", buf);
754			error(COMMAND, errmsg);
755			return;
756			}
757	schorsch	2.40	SET_FILE_BINARY(fp);
758	greg	1.1	(*dev->comout)("writing \"");
759			(*dev->comout)(fname);
760			(*dev->comout)("\"...\n");
761			/* write header */
762	greg	2.24	newheader("RADIANCE", fp);
763	greg	1.1	fputs(progname, fp);
764			fprintview(&ourview, fp);
765	greg	1.29	if (octname != NULL)
766			fprintf(fp, " %s\n", octname);
767			else
768			putc('\n', fp);
769	greg	1.28	fprintf(fp, "SOFTWARE= %s\n", VersionID);
770	greg	2.38	fputnow(fp);
771	greg	1.1	if (exposure != 1.0)
772	greg	1.13	fputexpos(exposure, fp);
773	greg	1.15	if (dev->pixaspect != 1.0)
774			fputaspect(dev->pixaspect, fp);
775	greg	1.22	fputformat(COLRFMT, fp);
776	greg	1.11	putc('\n', fp);
777	greg	1.33	fprtresolu(hresolu, vresolu, fp);
778	greg	1.1
779	greg	1.14	scanline = (COLR )malloc(hresolusizeof(COLR));
780	greg	1.27	if (scanline == NULL) {
781			error(COMMAND, "not enough memory!");
782			fclose(fp);
783			unlink(fname);
784			return;
785			}
786	greg	1.14	for (y = vresolu-1; y >= 0; y--) {
787			getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu);
788			if (fwritecolrs(scanline, hresolu, fp) < 0)
789	greg	1.1	break;
790			}
791	greg	2.38	free((void *)scanline);
792	greg	1.1	if (fclose(fp) < 0)
793			error(COMMAND, "write error");
794			}