Index: trunk/psLib/src/imageops/psImageStats.c
===================================================================
--- trunk/psLib/src/imageops/psImageStats.c	(revision 2294)
+++ trunk/psLib/src/imageops/psImageStats.c	(revision 2436)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.54 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-11-06 00:57:44 $
+ *  @version $Revision: 1.55 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-11-24 22:14:39 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -183,6 +183,11 @@
     sums[i][j]: This will contain the sum of 
                 input->data.F32[x][y] *
-                psPolynomial1DEval((float) x, chebPolys[i]) *
-                psPolynomial1DEval((float) y, chebPolys[j]);
+                psPolynomial1DEval(
+chebPolys[i],
+(float) x) *
+                psPolynomial1DEval(
+chebPolys[j],
+(float) y, 
+);
         over all pixels (x,y) in the image.
   *****************************************************************************/
@@ -249,6 +254,12 @@
                     }
                     sums[i][j] += pixel *
-                                  psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) *
-                                  psPolynomial1DEval(cScalingFactors[y], chebPolys[j]);
+                                  psPolynomial1DEval(
+                                      chebPolys[i],
+                                      rScalingFactors[x]
+                                  ) *
+                                  psPolynomial1DEval(
+                                      chebPolys[j],
+                                      cScalingFactors[y]
+                                  );
                 }
             }
@@ -336,6 +347,13 @@
                 for (j = 0; j < coeffs->nY; j++) {
                     polySum +=
-                        psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) *
-                        psPolynomial1DEval(cScalingFactors[y], chebPolys[j]) *
+                        psPolynomial1DEval(
+                            chebPolys[i],
+                            rScalingFactors[x]
+                        ) *
+                        psPolynomial1DEval(
+                            chebPolys[j],
+                            cScalingFactors[y]
+
+                        ) *
                         coeffs->coeff[i][j];
 
