--- ray/src/rt/t_func.c 1989/02/21 14:56:22 1.2 +++ ray/src/rt/t_func.c 1989/04/19 21:20:29 1.6 @@ -43,6 +43,7 @@ register RAY *r; extern double varvalue(); extern int errno; FVECT disp; + double d; register XFORM *mxf; register int i; register char **sa; @@ -71,7 +72,7 @@ register RAY *r; m->os = (char *)mxf; } - setmap(m, r, mxf->back.sca, mxf->back.xfm); + setmap(m, r, mxf->back.xfm, mxf->back.sca); if (!vardefined(sa[0])) loadfunc(sa[3]); @@ -82,15 +83,11 @@ register RAY *r; objerror(m, WARNING, "compute error"); return; } - if (mxf->fore.xfm == NULL) - for (i = 0; i < 3; i++) - r->pert[i] += disp[i]; - else - for (i = 0; i < 3; i++) - r->pert[i] += ( disp[0]*mxf->fore.xfm[0][i] + - disp[1]*mxf->fore.xfm[1][i] + - disp[2]*mxf->fore.xfm[2][i] ) - / mxf->fore.sca; + multv3(disp, disp, mxf->fore.xfm); + multv3(disp, disp, r->rofx); + d = 1.0 / (mxf->fore.sca * r->rofs); + for (i = 0; i < 3; i++) + r->pert[i] += disp[i] * d; return; memerr: error(SYSTEM, "out of memory in t_func");