--- ray/src/rt/func.c 1989/02/02 10:41:23 1.1 +++ ray/src/rt/func.c 1991/06/17 08:28:26 1.16 @@ -1,4 +1,4 @@ -/* Copyright (c) 1986 Regents of the University of California */ +/* Copyright (c) 1991 Regents of the University of California */ #ifndef lint static char SCCSid[] = "$SunId$ LBL"; @@ -15,46 +15,57 @@ static char SCCSid[] = "$SunId$ LBL"; #include "otypes.h" -typedef struct { - double sca; /* scalefactor */ - double xfm[4][4]; /* transform matrix */ -} XF; +XF unitxf = { /* identity transform */ + 1.0, 0.0, 0.0, 0.0, + 0.0, 1.0, 0.0, 0.0, + 0.0, 0.0, 1.0, 0.0, + 0.0, 0.0, 0.0, 1.0, + 1.0 +}; -static OBJREC *fobj; /* current function object */ -static RAY *fray; /* current function ray */ -static XF fxf; /* current transformation */ +XF funcxf; /* current transformation */ +static OBJREC *fobj = NULL; /* current function object */ +static RAY *fray = NULL; /* current function ray */ -setmap(m, r, sca, xfm) /* set channels for function call */ +setmap(m, r, bx) /* set channels for function call */ OBJREC *m; register RAY *r; -double sca; -double xfm[4][4]; +XF *bx; { - extern double l_noise3(), l_noise3a(), l_noise3b(), l_noise3c(); - extern double l_hermite(), l_fnoise3(), l_arg(); + extern double l_arg(); extern long eclock; static char *initfile = "rayinit.cal"; - + static long lastrno = -1; + /* check to see if already set */ + if (m == fobj && r->rno == lastrno) + return; + /* initialize if first call */ if (initfile != NULL) { loadfunc(initfile); - scompile(NULL, "Dx=$1;Dy=$2;Dz=$3;"); - scompile(NULL, "Nx=$4;Ny=$5;Nz=$6;"); - scompile(NULL, "Px=$7;Py=$8;Pz=$9;"); - scompile(NULL, "T=$10;Rdot=$11;"); - funset("arg", 1, l_arg); - funset("noise3", 3, l_noise3); - funset("noise3a", 3, l_noise3a); - funset("noise3b", 3, l_noise3b); - funset("noise3c", 3, l_noise3c); - funset("hermite", 5, l_hermite); - funset("fnoise3", 3, l_fnoise3); + scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0); + scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0); + scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0); + scompile("T=$10;Rdot=$11;", NULL, 0); + scompile("S=$12;Tx=$13;Ty=$14;Tz=$15;", NULL, 0); + scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0); + scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0); + scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0); + funset("arg", 1, '=', l_arg); + setnoisefuncs(); initfile = NULL; } fobj = m; fray = r; - fxf.sca = r->ros * sca; - multmat4(fxf.xfm, r->rox, xfm); + if (r->rox != NULL) + if (bx != &unitxf) { + funcxf.sca = r->rox->b.sca * bx->sca; + multmat4(funcxf.xfm, r->rox->b.xfm, bx->xfm); + } else + copystruct(&funcxf, &r->rox->b); + else + copystruct(&funcxf, bx); + lastrno = r->rno; eclock++; /* notify expression evaluator */ } @@ -63,31 +74,32 @@ setfunc(m, r) /* simplified interface to setmap */ register OBJREC *m; RAY *r; { -#define mxf ((XF *)m->os) + register XF *mxf; - if (mxf == NULL) { - register int n = m->oargs.nsargs; - register char **sa = m->oargs.sarg; + if ((mxf = (XF *)m->os) == NULL) { + register int n; + register char **sa; - while (n > 0 && **sa != '-') { - n--; - sa++; + for (n = m->oargs.nsargs, sa = m->oargs.sarg; + n > 0 && **sa != '-'; n--, sa++) + ; + if (n == 0) + mxf = &unitxf; + else { + mxf = (XF *)malloc(sizeof(XF)); + if (mxf == NULL) + goto memerr; + if (invxf(mxf, n, sa) != n) + objerror(m, USER, "bad transform"); + if (mxf->sca < 0.0) + mxf->sca = -mxf->sca; } - mxf = (XF *)malloc(sizeof(XF)); - if (mxf == NULL) - goto memerr; - mxf->sca = 1.0; - setident4(mxf->xfm); - if (invxf(mxf->xfm, &mxf->sca, n, sa) != n) - objerror(m, USER, "bad transform"); - if (mxf->sca < 0.0) - mxf->sca = -mxf->sca; + m->os = (char *)mxf; } - setmap(m, r, mxf->sca, mxf->xfm); + setmap(m, r, mxf); return; memerr: error(SYSTEM, "out of memory in setfunc"); -#undef mxf } @@ -97,7 +109,7 @@ char *fname; extern char *libpath; /* library search path */ char *ffname; - if ((ffname = getpath(fname, libpath)) == NULL) { + if ((ffname = getpath(fname, libpath, R_OK)) == NULL) { sprintf(errmsg, "cannot find function file \"%s\"", fname); error(USER, errmsg); } @@ -128,41 +140,58 @@ double chanvalue(n) /* return channel n to calcomp */ register int n; { - double res; + double sum; register RAY *r; - n--; /* for convenience */ + if (--n < 0) + goto badchan; - if (n < 0 || n > 10) - error(USER, "illegal channel number"); + if (n < 3) /* ray direction */ - if (n == 9) { /* distance */ + return( ( fray->rdir[0]*funcxf.xfm[0][n] + + fray->rdir[1]*funcxf.xfm[1][n] + + fray->rdir[2]*funcxf.xfm[2][n] ) + / funcxf.sca ); - res = fray->rot; - for (r = fray->parent; r != NULL; r = r->parent) - res += r->rot; - res *= fxf.sca; + if (n < 6) /* surface normal */ - } else if (n == 10) { /* dot product */ + return( ( fray->ron[0]*funcxf.xfm[0][n-3] + + fray->ron[1]*funcxf.xfm[1][n-3] + + fray->ron[2]*funcxf.xfm[2][n-3] ) + / funcxf.sca ); - res = fray->rod; + if (n < 9) /* intersection */ - } else if (n < 3) { /* ray direction */ - res = ( fray->rdir[0]*fxf.xfm[0][n] + - fray->rdir[1]*fxf.xfm[1][n] + - fray->rdir[2]*fxf.xfm[2][n] ) - / fxf.sca ; - } else if (n < 6) { /* surface normal */ - res = ( fray->ron[0]*fxf.xfm[0][n-3] + - fray->ron[1]*fxf.xfm[1][n-3] + - fray->ron[2]*fxf.xfm[2][n-3] ) - / fxf.sca ; - } else { /* intersection */ - res = fray->rop[0]*fxf.xfm[0][n-6] + - fray->rop[1]*fxf.xfm[1][n-6] + - fray->rop[2]*fxf.xfm[2][n-6] + - fxf.xfm[3][n-6] ; + return( fray->rop[0]*funcxf.xfm[0][n-6] + + fray->rop[1]*funcxf.xfm[1][n-6] + + fray->rop[2]*funcxf.xfm[2][n-6] + + funcxf.xfm[3][n-6] ); + + if (n == 9) { /* distance */ + + sum = fray->rot; + for (r = fray->parent; r != NULL; r = r->parent) + sum += r->rot; + return(sum * funcxf.sca); + } + if (n == 10) /* dot product */ + return(fray->rod); - return(res); + if (n == 11) /* scale */ + return(funcxf.sca); + + if (n < 15) /* origin */ + return(funcxf.xfm[3][n-12]); + + if (n < 18) /* i unit vector */ + return(funcxf.xfm[0][n-15] / funcxf.sca); + + if (n < 21) /* j unit vector */ + return(funcxf.xfm[1][n-15] / funcxf.sca); + + if (n < 24) /* k unit vector */ + return(funcxf.xfm[2][n-21] / funcxf.sca); +badchan: + error(USER, "illegal channel number"); }