Index: /trunk/psLib/src/dataManip/psFunctions.c
===================================================================
--- /trunk/psLib/src/dataManip/psFunctions.c	(revision 2788)
+++ /trunk/psLib/src/dataManip/psFunctions.c	(revision 2789)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.75 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-12-22 05:09:32 $
+ *  @version $Revision: 1.76 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-22 21:00:36 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -297,6 +297,6 @@
                     "polysum+= sum*coeff [%f+= (%f * %f)\n", polySum, xSum, myPoly->coeff[loop_x]);
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -370,6 +370,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -431,6 +431,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -504,6 +504,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;
@@ -577,6 +577,6 @@
         if (myPoly->mask[loop_x] == 0) {
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -624,6 +624,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -684,6 +684,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -757,6 +757,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;
Index: /trunk/psLib/src/math/psPolynomial.c
===================================================================
--- /trunk/psLib/src/math/psPolynomial.c	(revision 2788)
+++ /trunk/psLib/src/math/psPolynomial.c	(revision 2789)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.75 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-12-22 05:09:32 $
+ *  @version $Revision: 1.76 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-22 21:00:36 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -297,6 +297,6 @@
                     "polysum+= sum*coeff [%f+= (%f * %f)\n", polySum, xSum, myPoly->coeff[loop_x]);
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -370,6 +370,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -431,6 +431,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -504,6 +504,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;
@@ -577,6 +577,6 @@
         if (myPoly->mask[loop_x] == 0) {
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -624,6 +624,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -684,6 +684,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -757,6 +757,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;
Index: /trunk/psLib/src/math/psSpline.c
===================================================================
--- /trunk/psLib/src/math/psSpline.c	(revision 2788)
+++ /trunk/psLib/src/math/psSpline.c	(revision 2789)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.75 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-12-22 05:09:32 $
+ *  @version $Revision: 1.76 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-22 21:00:36 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -297,6 +297,6 @@
                     "polysum+= sum*coeff [%f+= (%f * %f)\n", polySum, xSum, myPoly->coeff[loop_x]);
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -370,6 +370,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -431,6 +431,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -504,6 +504,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;
@@ -577,6 +577,6 @@
         if (myPoly->mask[loop_x] == 0) {
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -624,6 +624,6 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -684,6 +684,6 @@
                 if (myPoly->mask[loop_x][loop_y][loop_z] == 0) {
                     polySum += zSum * myPoly->coeff[loop_x][loop_y][loop_z];
-                    zSum *= z;
                 }
+                zSum *= z;
             }
             ySum *= y;
@@ -757,6 +757,6 @@
                     if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) {
                         polySum += zSum * myPoly->coeff[loop_w][loop_x][loop_y][loop_z];
-                        zSum *= z;
                     }
+                    zSum *= z;
                 }
                 ySum *= y;