Index: trunk/psLib/src/math/psPolynomial.c
===================================================================
--- trunk/psLib/src/math/psPolynomial.c	(revision 3397)
+++ trunk/psLib/src/math/psPolynomial.c	(revision 3405)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.93 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-03-10 01:11:12 $
+ *  @version $Revision: 1.94 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-03-10 22:33:57 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1746,12 +1746,7 @@
 
     psVector *tmp;
-    psVector *myW;
-    psVector *myX;
-    psVector *myY;
-    psVector *myZ;
-    psS32 i;
-    psS32 vecLen=x->n;
-
-    vecLen=w->n;
+    psS32 vecLen=w->n;
+
+    // Determine the output vector size from min of input vectors
     if (y->n < vecLen) {
         vecLen = y->n;
@@ -1764,26 +1759,17 @@
     }
 
+    // Allocate output vector
     tmp = psVectorAlloc(vecLen, PS_TYPE_F64);
-    for (i=0;i<vecLen;i++) {
+
+    // Evaluate the polynomial
+    for (psS32 i = 0; i < vecLen; i++) {
         tmp->data.F64[i] = psDPolynomial4DEval(myPoly,
-                                               myW->data.F64[i],
-                                               myX->data.F64[i],
-                                               myY->data.F64[i],
-                                               myZ->data.F64[i]);
-
-    }
-
-    if (w->type.type != PS_TYPE_F64) {
-        psFree(myW);
-    }
-    if (x->type.type != PS_TYPE_F64) {
-        psFree(myX);
-    }
-    if (y->type.type != PS_TYPE_F64) {
-        psFree(myY);
-    }
-    if (z->type.type != PS_TYPE_F64) {
-        psFree(myZ);
-    }
+                                               w->data.F64[i],
+                                               x->data.F64[i],
+                                               y->data.F64[i],
+                                               z->data.F64[i]);
+    }
+
+    // Return output vector
     return(tmp);
 }