--- ray/src/rt/t_data.c 1989/10/05 07:23:19 1.5 +++ ray/src/rt/t_data.c 1994/01/12 16:46:53 2.4 @@ -1,4 +1,4 @@ -/* Copyright (c) 1989 Regents of the University of California */ +/* Copyright (c) 1992 Regents of the University of California */ #ifndef lint static char SCCSid[] = "$SunId$ LBL"; @@ -14,6 +14,8 @@ static char SCCSid[] = "$SunId$ LBL"; #include "data.h" +#include "func.h" + /* * A stored texture is specified as follows: * @@ -31,86 +33,53 @@ static char SCCSid[] = "$SunId$ LBL"; * to get from the original coordinates to the current coordinates. */ -typedef struct { - struct { - double sca; /* scale factor */ - double xfm[4][4]; /* transformation matrix */ - } fore, back; -} XFORM; - t_data(m, r) /* interpolate texture data */ register OBJREC *m; RAY *r; { - extern double varvalue(), funvalue(), datavalue(); - extern int errno; int nv; - FVECT dval, disp; - double pt[MAXDIM]; + FVECT disp; + double dval[3], pt[MAXDIM]; double d; DATARRAY *dp; - register XFORM *mxf; - register char **sa; + register MFUNC *mf; register int i; if (m->oargs.nsargs < 8) objerror(m, USER, "bad # arguments"); - sa = m->oargs.sarg; - - for (i = 7; i < m->oargs.nsargs && sa[i][0] != '-'; i++) - ; - nv = i-7; - if ((mxf = (XFORM *)m->os) == NULL) { - mxf = (XFORM *)malloc(sizeof(XFORM)); - if (mxf == NULL) - goto memerr; - if (xf(mxf->fore.xfm, &mxf->fore.sca, - m->oargs.nsargs-i, sa+i) != m->oargs.nsargs-i) - objerror(m, USER, "bad transform"); - if (mxf->fore.sca < 0.0) - mxf->fore.sca = -mxf->fore.sca; - invxf(mxf->back.xfm, &mxf->back.sca, - m->oargs.nsargs-i, sa+i); - if (mxf->back.sca < 0.0) - mxf->back.sca = -mxf->back.sca; - m->os = (char *)mxf; - } - - setmap(m, r, mxf->back.xfm, mxf->back.sca); - - if (nv > MAXDIM) - goto dimerr; - if (!vardefined(sa[7])) - loadfunc(sa[6]); + dp = getdata(m->oargs.sarg[3]); + i = (1 << (nv = dp->nd)) - 1; + mf = getfunc(m, 6, i<<7, 1); + setfunc(m, r); errno = 0; for (i = 0; i < nv; i++) - pt[i] = varvalue(sa[i+7]); + pt[i] = evalue(mf->ep[i]); if (errno) goto computerr; - for (i = 0; i < 3; i++) { - dp = getdata(sa[i+3]); + dval[0] = datavalue(dp, pt); + for (i = 1; i < 3; i++) { + dp = getdata(m->oargs.sarg[i+3]); if (dp->nd != nv) - goto dimerr; + objerror(m, USER, "dimension error"); dval[i] = datavalue(dp, pt); } errno = 0; for (i = 0; i < 3; i++) - disp[i] = funvalue(sa[i], 3, dval); + disp[i] = funvalue(m->oargs.sarg[i], 3, dval); if (errno) goto computerr; - - multv3(disp, disp, mxf->fore.xfm); - multv3(disp, disp, r->rofx); - d = 1.0 / (mxf->fore.sca * r->rofs); + if (mf->f != &unitxf) + multv3(disp, disp, mf->f->xfm); + if (r->rox != NULL) { + multv3(disp, disp, r->rox->f.xfm); + d = 1.0 / (mf->f->sca * r->rox->f.sca); + } else + d = 1.0 / mf->f->sca; for (i = 0; i < 3; i++) r->pert[i] += disp[i] * d; - return; -dimerr: - objerror(m, USER, "dimension error"); -memerr: - error(SYSTEM, "out of memory in t_data"); + return(0); computerr: objerror(m, WARNING, "compute error"); - return; + return(0); }