--- ray/src/rt/func.c 1990/04/07 16:12:56 1.8 +++ 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,51 +15,45 @@ static char SCCSid[] = "$SunId$ LBL"; #include "otypes.h" -typedef struct { - double xfm[4][4]; /* transform matrix */ - double sca; /* scalefactor */ -} 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, xfm, sca) /* set channels for function call */ + +setmap(m, r, bx) /* set channels for function call */ OBJREC *m; register RAY *r; -double xfm[4][4]; -double sca; +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;"); - scompile(NULL, "S=$12;Tx=$13;Ty=$14;Tz=$15;"); - scompile(NULL, "Ix=$16;Iy=$17;Iz=$18;"); - scompile(NULL, "Jx=$19;Jy=$20;Jz=$21;"); - scompile(NULL, "Kx=$22;Ky=$23;Kz=$24;"); - 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->robs * sca; - multmat4(fxf.xfm, r->robx, 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); } @@ -76,20 +70,22 @@ RAY *r; for (n = m->oargs.nsargs, sa = m->oargs.sarg; n > 0 && **sa != '-'; n--, sa++) ; - mxf = (XF *)malloc(sizeof(XF)); - if (mxf == NULL) - goto memerr; - if (invxf(mxf->xfm, &mxf->sca, n, sa) != n) - objerror(m, USER, "bad transform"); - if (mxf->sca < 0.0) - mxf->sca = -mxf->sca; + 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; + } m->os = (char *)mxf; } - setmap(m, r, mxf->xfm, mxf->sca); - return; + return(setmap(m, r, mxf)); memerr: error(SYSTEM, "out of memory in setfunc"); -#undef mxf } @@ -127,61 +123,112 @@ l_arg() /* return nth real argument */ double -chanvalue(n) /* return channel n to calcomp */ -register int n; +l_erf() /* error function */ { - double sum; - register RAY *r; + extern double erf(); - n--; /* for convenience */ + return(erf(argument(1))); +} - if (n < 0 || n > 23) - error(USER, "illegal channel number"); - if (n == 9) { /* distance */ +double +l_erfc() /* cumulative error function */ +{ + extern double erfc(); - sum = fray->rot; - for (r = fray->parent; r != NULL; r = r->parent) - sum += r->rot; - return(sum * fxf.sca); + 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'; } - if (n == 10) /* dot product */ - return(fray->rod); + setcontext(fname); + if (!vardefined(DEFVNAME)) { + loadfunc(fname); + varset(DEFVNAME, ':', 1.0); + } +} + +double +chanvalue(n) /* return channel n to calcomp */ +register int n; +{ + double sum; + register RAY *r; + + if (--n < 0) + goto badchan; + if (n < 3) /* ray direction */ - return( ( fray->rdir[0]*fxf.xfm[0][n] + - fray->rdir[1]*fxf.xfm[1][n] + - fray->rdir[2]*fxf.xfm[2][n] ) - / fxf.sca ); + return( ( fray->rdir[0]*funcxf.xfm[0][n] + + fray->rdir[1]*funcxf.xfm[1][n] + + fray->rdir[2]*funcxf.xfm[2][n] ) + / funcxf.sca ); if (n < 6) /* surface normal */ - return( ( 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 ); + 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 ); if (n < 9) /* intersection */ - return( 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 ); + if (n == 11) /* scale */ - return(fxf.sca); + return(funcxf.sca); if (n < 15) /* origin */ - return(fxf.xfm[3][n-12]); + return(funcxf.xfm[3][n-12]); if (n < 18) /* i unit vector */ - return(fxf.xfm[0][n-15] / fxf.sca); + return(funcxf.xfm[0][n-15] / funcxf.sca); if (n < 21) /* j unit vector */ - return(fxf.xfm[1][n-15] / fxf.sca); + return(funcxf.xfm[1][n-15] / funcxf.sca); if (n < 24) /* k unit vector */ - return(fxf.xfm[2][n-21] / fxf.sca); + return(funcxf.xfm[2][n-21] / funcxf.sca); +badchan: + error(USER, "illegal channel number"); }