src/rt/rv2.c

/* Copyright (c) 1995 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  *atos();

extern FILE  *popen();

extern char  rifname[128];              /* rad input file name */

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


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


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


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


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


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


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


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


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 int  ambvwt;
        extern double  ambacc;
        extern int  ambres;
        extern int  ambdiv;
        extern int  ambssamp;
        extern int  ambounce;
        extern COLOR  cextinction;
        extern COLOR  salbedo;
        extern double  seccg;
        extern double  ssampdist;
        extern int  directvis;
        extern double  srcsizerat;
        extern int  do_irrad;
        extern double  specjitter;
        extern double  specthresh;
        extern int  backvis;
        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', &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':                       /* back faces or black and white */
                switch (s[1]) {
                case 'v':                       /* back face visibility */
                        getparam(s+2, "back face visibility", 'b', &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', &greyscale);
                        break;
                default:
                        goto badparam;
                }
                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 'w':                       /* weight */
                        getparam(s+2, "ambient value weight", 'i', &ambvwt);
                        break;
                case 'a':                       /* accuracy */
                        if (getparam(s+2, "ambient accuracy", 'r', &ambacc))
                                setambacc(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 'm':                       /* medium */
                switch (s[1]) {
                case 'e':                       /* extinction coefficient */
                        getparam(s+2, "extinction coefficient", 'C',
                                        (COLOR *)cextinction);
                        break;
                case 'a':                       /* scattering albedo */
                        getparam(s+2, "scattering albedo", 'C',
                                        (COLOR *)salbedo);
                        break;
                case 'g':                       /* scattering eccentricity */
                        getparam(s+2, "scattering eccentricity", 'r', &seccg);
                        break;
                case 's':                       /* sampling distance */
                        getparam(s+2, "mist sampling distance", 'r',
                                        &ssampdist);
                        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;

        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);
                (*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) (%.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;

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