Changeset 2789 for trunk/psLib/src/math/psPolynomial.c

Timestamp:

Dec 22, 2004, 11:00:36 AM (22 years ago)

Author:

gusciora

Message:

Fixed bug 254.

File:

• trunk/psLib/src/math/psPolynomial.c (modified) (9 diffs)

trunk/psLib/src/math/psPolynomial.c

@@ -7 +7 @@
  *  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 +297 @@
                     "polysum+= sum*coeff [%f+= (%f * %f)\n", polySum, xSum, myPoly->coeff[loop_x]);
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -370 +370 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -431 +431 @@
                 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 +504 @@
                     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 +577 @@
         if (myPoly->mask[loop_x] == 0) {
             polySum += xSum * myPoly->coeff[loop_x];
-            xSum *= x;
-        }
+        }
+        xSum *= x;
     }
 
@@ -624 +624 @@
             if (myPoly->mask[loop_x][loop_y] == 0) {
                 polySum += ySum * myPoly->coeff[loop_x][loop_y];
-                ySum *= y;
             }
+            ySum *= y;
         }
         xSum *= x;
@@ -684 +684 @@
                 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 +757 @@
                     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;