--- ray/src/rt/func.c 1989/02/21 14:56:17 1.2 +++ ray/src/rt/func.c 1991/08/08 13:08:01 1.21 @@ -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,47 +15,45 @@ static char SCCSid[] = "$SunId$ LBL"; #include "otypes.h" -typedef struct { - double sca; /* scalefactor */ - double xfm[4][4]; /* transform matrix */ -} XF; +#define INITFILE "rayinit.cal" +#define DEFVNAME "FILE_LOADED`" -static OBJREC *fobj; /* current function object */ -static RAY *fray; /* current function ray */ -static XF fxf; /* current transformation */ +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 +}; +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 long eclock; - static char *initfile = "rayinit.cal"; - - 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); - initfile = NULL; - } + static long lastrno = -1; + /* check to see if already set */ + if (m == fobj && r->rno == lastrno) + return(0); 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 */ + return(1); } @@ -66,29 +64,28 @@ RAY *r; register XF *mxf; if ((mxf = (XF *)m->os) == NULL) { - register int n = m->oargs.nsargs; - register char **sa = m->oargs.sarg; + 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); - return; + return(setmap(m, r, mxf)); memerr: error(SYSTEM, "out of memory in setfunc"); -#undef mxf } @@ -98,7 +95,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); } @@ -126,44 +123,112 @@ l_arg() /* return nth real argument */ double +l_erf() /* error function */ +{ + extern double erf(); + + return(erf(argument(1))); +} + + +double +l_erfc() /* cumulative error function */ +{ + extern double erfc(); + + return(erfc(argument(1))); +} + + +funcfile(fname) /* set context, load file if necessary */ +char *fname; +{ + extern char *setcontext(); + static char initfile[] = INITFILE; + + if (initfile[0]) { /* initialize on first call */ + setcontext(""); + 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); + funset("erf", 1, ':', l_erf); + funset("erfc", 1, ':', l_erfc); + setnoisefuncs(); + loadfunc(initfile); + initfile[0] = '\0'; + } + setcontext(fname); + if (!vardefined(DEFVNAME)) { + loadfunc(fname); + varset(DEFVNAME, ':', 1.0); + } +} + + +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 (range [-1,1]) */ + return( fray->rod <= -1.0 ? -1.0 : + fray->rod >= 1.0 ? 1.0 : + 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"); }