--- ray/src/common/xf.c	1989/10/13 19:34:51	1.6
+++ ray/src/common/xf.c	1990/02/20 18:02:30	1.8
@@ -26,12 +26,13 @@ char  *av[];
 {
 	double  atof(), sin(), cos();
 	double  xfmat[4][4], m4[4][4];
-	double  xfsca, theta;
+	double  xfsca, dtmp;
 	int  i, icnt;
 
 	setident4(retmat);
 	*retsca = 1.0;
 
+	icnt = 1;
 	setident4(xfmat);
 	xfsca = 1.0;
 
@@ -52,21 +53,21 @@ char  *av[];
 			switch (av[i][2]) {
 			case 'x':
 				checkarg(3,1);
-				theta = d2r(atof(av[++i]));
-				m4[1][1] = m4[2][2] = cos(theta);
-				m4[2][1] = -(m4[1][2] = sin(theta));
+				dtmp = d2r(atof(av[++i]));
+				m4[1][1] = m4[2][2] = cos(dtmp);
+				m4[2][1] = -(m4[1][2] = sin(dtmp));
 				break;
 			case 'y':
 				checkarg(3,1);
-				theta = d2r(atof(av[++i]));
-				m4[0][0] = m4[2][2] = cos(theta);
-				m4[0][2] = -(m4[2][0] = sin(theta));
+				dtmp = d2r(atof(av[++i]));
+				m4[0][0] = m4[2][2] = cos(dtmp);
+				m4[0][2] = -(m4[2][0] = sin(dtmp));
 				break;
 			case 'z':
 				checkarg(3,1);
-				theta = d2r(atof(av[++i]));
-				m4[0][0] = m4[1][1] = cos(theta);
-				m4[1][0] = -(m4[0][1] = sin(theta));
+				dtmp = d2r(atof(av[++i]));
+				m4[0][0] = m4[1][1] = cos(dtmp);
+				m4[1][0] = -(m4[0][1] = sin(dtmp));
 				break;
 			default:
 				return(i);
@@ -75,10 +76,13 @@ char  *av[];
 
 		case 's':			/* scale */
 			checkarg(2,1);
+			dtmp = atof(av[i+1]);
+			if (dtmp == 0.0) goto done;
+			i++;
 			xfsca *=
 			m4[0][0] = 
 			m4[1][1] = 
-			m4[2][2] = atof(av[++i]);
+			m4[2][2] = dtmp;
 			break;
 
 		case 'm':			/* mirror */
@@ -105,14 +109,14 @@ char  *av[];
 
 		case 'i':			/* iterate */
 			checkarg(2,1);
-			icnt = atoi(av[++i]);
 			while (icnt-- > 0) {
 				multmat4(retmat, retmat, xfmat);
 				*retsca *= xfsca;
 			}
+			icnt = atoi(av[++i]);
 			setident4(xfmat);
 			xfsca = 1.0;
-			break;
+			continue;
 
 		default:
 			return(i);
@@ -121,8 +125,10 @@ char  *av[];
 		multmat4(xfmat, xfmat, m4);
 	}
 done:
-	multmat4(retmat, retmat, xfmat);
-	*retsca *= xfsca;
+	while (icnt-- > 0) {
+		multmat4(retmat, retmat, xfmat);
+		*retsca *= xfsca;
+	}
 	return(i);
 }
 
@@ -137,12 +143,13 @@ char  *av[];
 {
 	double  atof(), sin(), cos();
 	double  xfmat[4][4], m4[4][4];
-	double  xfsca, theta;
+	double  xfsca, dtmp;
 	int  i, icnt;
 
 	setident4(retmat);
 	*retsca = 1.0;
 
+	icnt = 1;
 	setident4(xfmat);
 	xfsca = 1.0;
 
@@ -163,21 +170,21 @@ char  *av[];
 			switch (av[i][2]) {
 			case 'x':
 				checkarg(3,1);
-				theta = -d2r(atof(av[++i]));
-				m4[1][1] = m4[2][2] = cos(theta);
-				m4[2][1] = -(m4[1][2] = sin(theta));
+				dtmp = -d2r(atof(av[++i]));
+				m4[1][1] = m4[2][2] = cos(dtmp);
+				m4[2][1] = -(m4[1][2] = sin(dtmp));
 				break;
 			case 'y':
 				checkarg(3,1);
-				theta = -d2r(atof(av[++i]));
-				m4[0][0] = m4[2][2] = cos(theta);
-				m4[0][2] = -(m4[2][0] = sin(theta));
+				dtmp = -d2r(atof(av[++i]));
+				m4[0][0] = m4[2][2] = cos(dtmp);
+				m4[0][2] = -(m4[2][0] = sin(dtmp));
 				break;
 			case 'z':
 				checkarg(3,1);
-				theta = -d2r(atof(av[++i]));
-				m4[0][0] = m4[1][1] = cos(theta);
-				m4[1][0] = -(m4[0][1] = sin(theta));
+				dtmp = -d2r(atof(av[++i]));
+				m4[0][0] = m4[1][1] = cos(dtmp);
+				m4[1][0] = -(m4[0][1] = sin(dtmp));
 				break;
 			default:
 				return(i);
@@ -186,10 +193,13 @@ char  *av[];
 
 		case 's':			/* scale */
 			checkarg(2,1);
+			dtmp = atof(av[i+1]);
+			if (dtmp == 0.0) goto done;
+			i++;
 			xfsca *=
 			m4[0][0] = 
 			m4[1][1] = 
-			m4[2][2] = 1.0 / atof(av[++i]);
+			m4[2][2] = 1.0 / dtmp;
 			break;
 
 		case 'm':			/* mirror */
@@ -216,11 +226,11 @@ char  *av[];
 
 		case 'i':			/* iterate */
 			checkarg(2,1);
-			icnt = atoi(av[++i]);
 			while (icnt-- > 0) {
 				multmat4(retmat, xfmat, retmat);
 				*retsca *= xfsca;
 			}
+			icnt = atoi(av[++i]);
 			setident4(xfmat);
 			xfsca = 1.0;
 			break;
@@ -232,8 +242,10 @@ char  *av[];
 		multmat4(xfmat, m4, xfmat);	/* left multiply */
 	}
 done:
-	multmat4(retmat, xfmat, retmat);
-	*retsca *= xfsca;
+	while (icnt-- > 0) {
+		multmat4(retmat, xfmat, retmat);
+		*retsca *= xfsca;
+	}
 	return(i);
 }
 #endif