--- ray/src/util/rmatrix.c	2015/05/04 23:27:04	2.14
+++ ray/src/util/rmatrix.c	2015/07/22 04:29:56	2.16
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: rmatrix.c,v 2.14 2015/05/04 23:27:04 greg Exp $";
+static const char RCSid[] = "$Id: rmatrix.c,v 2.16 2015/07/22 04:29:56 greg Exp $";
 #endif
 /*
  * General matrix operations.
@@ -10,6 +10,7 @@ static const char RCSid[] = "$Id: rmatrix.c,v 2.14 201
 #include <string.h>
 #include <fcntl.h>
 #include "resolu.h"
+#include "rtprocess.h"
 #include "rmatrix.h"
 
 static char	rmx_mismatch_warn[] = "WARNING: data type mismatch\n";
@@ -513,8 +514,7 @@ rmx_multiply(const RMATRIX *m1, const RMATRIX *m2)
 	        for (k = mres->ncomp; k--; ) {
 		    long double	d = 0;
 		    for (h = m1->ncols; h--; )
-			d += (long double)rmx_lval(m1,i,h,k) *
-				(long double)rmx_lval(m2,h,j,k);
+			d += rmx_lval(m1,i,h,k) * rmx_lval(m2,h,j,k);
 		    rmx_lval(mres,i,j,k) = (double)d;
 		}
 	return(mres);