--- ray/src/rt/t_func.c	1989/02/21 14:56:22	1.2
+++ ray/src/rt/t_func.c	1989/10/05 07:23:20	1.7
@@ -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;
@@ -55,15 +56,11 @@ register RAY  *r;
 		mxf = (XFORM *)malloc(sizeof(XFORM));
 		if (mxf == NULL)
 			goto memerr;
-		mxf->fore.sca = 1.0;
-		setident4(mxf->fore.xfm);
 		if (xf(mxf->fore.xfm, &mxf->fore.sca,
 			m->oargs.nsargs-4, sa+4) != m->oargs.nsargs-4)
 			objerror(m, USER, "bad transform");
 		if (mxf->fore.sca < 0.0)
 			mxf->fore.sca = -mxf->fore.sca;
-		mxf->back.sca = 1.0;
-		setident4(mxf->back.xfm);
 		invxf(mxf->back.xfm, &mxf->back.sca,
 				m->oargs.nsargs-4, sa+4);
 		if (mxf->back.sca < 0.0)
@@ -71,7 +68,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 +79,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");