Index: trunk/psLib/src/math/psPolynomial.c
===================================================================
--- trunk/psLib/src/math/psPolynomial.c	(revision 1407)
+++ trunk/psLib/src/math/psPolynomial.c	(revision 1440)
@@ -8,6 +8,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.21 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-08-07 00:06:06 $
+ *  @version $Revision: 1.22 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-08-09 23:34:57 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -49,12 +49,12 @@
 
 /*****************************************************************************/
-static void polynomial1DFree(psPolynomial1D * myPoly);
-static void polynomial2DFree(psPolynomial2D * myPoly);
-static void polynomial3DFree(psPolynomial3D * myPoly);
-static void polynomial4DFree(psPolynomial4D * myPoly);
-static void dPolynomial1DFree(psDPolynomial1D * myPoly);
-static void dPolynomial2DFree(psDPolynomial2D * myPoly);
-static void dPolynomial3DFree(psDPolynomial3D * myPoly);
-static void dPolynomial4DFree(psDPolynomial4D * myPoly);
+static void polynomial1DFree(psPolynomial1D* myPoly);
+static void polynomial2DFree(psPolynomial2D* myPoly);
+static void polynomial3DFree(psPolynomial3D* myPoly);
+static void polynomial4DFree(psPolynomial4D* myPoly);
+static void dPolynomial1DFree(psDPolynomial1D* myPoly);
+static void dPolynomial2DFree(psDPolynomial2D* myPoly);
+static void dPolynomial3DFree(psDPolynomial3D* myPoly);
+static void dPolynomial4DFree(psDPolynomial4D* myPoly);
 
 /*****************************************************************************/
@@ -120,7 +120,7 @@
  NOTE: There is no way to seed the random generator.
  *****************************************************************************/
-psVector *psGaussianDev(float mean, float sigma, int Npts)
-{
-    psVector *gauss = NULL;
+psVector* psGaussianDev(float mean, float sigma, int Npts)
+{
+    psVector* gauss = NULL;
     const gsl_rng_type *T = NULL;
     gsl_rng *r = NULL;
@@ -144,10 +144,10 @@
     This routine must allocate memory for the polynomial structures.
  *****************************************************************************/
-psPolynomial1D *psPolynomial1DAlloc(int n)
+psPolynomial1D* psPolynomial1DAlloc(int n)
 {
     int i = 0;
-    psPolynomial1D *newPoly = NULL;
-
-    newPoly = (psPolynomial1D *) psAlloc(sizeof(psPolynomial1D));
+    psPolynomial1D* newPoly = NULL;
+
+    newPoly = (psPolynomial1D* ) psAlloc(sizeof(psPolynomial1D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial1DFree);
     newPoly->n = n;
@@ -164,11 +164,11 @@
 }
 
-psPolynomial2D *psPolynomial2DAlloc(int nX, int nY)
-{
-    int x = 0;
-    int y = 0;
-    psPolynomial2D *newPoly = NULL;
-
-    newPoly = (psPolynomial2D *) psAlloc(sizeof(psPolynomial2D));
+psPolynomial2D* psPolynomial2DAlloc(int nX, int nY)
+{
+    int x = 0;
+    int y = 0;
+    psPolynomial2D* newPoly = NULL;
+
+    newPoly = (psPolynomial2D* ) psAlloc(sizeof(psPolynomial2D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial2DFree);
     newPoly->nX = nX;
@@ -194,12 +194,12 @@
 }
 
-psPolynomial3D *psPolynomial3DAlloc(int nX, int nY, int nZ)
+psPolynomial3D* psPolynomial3DAlloc(int nX, int nY, int nZ)
 {
     int x = 0;
     int y = 0;
     int z = 0;
-    psPolynomial3D *newPoly = NULL;
-
-    newPoly = (psPolynomial3D *) psAlloc(sizeof(psPolynomial3D));
+    psPolynomial3D* newPoly = NULL;
+
+    newPoly = (psPolynomial3D* ) psAlloc(sizeof(psPolynomial3D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial3DFree);
     newPoly->nX = nX;
@@ -233,5 +233,5 @@
 }
 
-psPolynomial4D *psPolynomial4DAlloc(int nW, int nX, int nY, int nZ)
+psPolynomial4D* psPolynomial4DAlloc(int nW, int nX, int nY, int nZ)
 {
     int w = 0;
@@ -239,7 +239,7 @@
     int y = 0;
     int z = 0;
-    psPolynomial4D *newPoly = NULL;
-
-    newPoly = (psPolynomial4D *) psAlloc(sizeof(psPolynomial4D));
+    psPolynomial4D* newPoly = NULL;
+
+    newPoly = (psPolynomial4D* ) psAlloc(sizeof(psPolynomial4D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial4DFree);
     newPoly->nW = nW;
@@ -281,5 +281,5 @@
 }
 
-static void polynomial1DFree(psPolynomial1D * myPoly)
+static void polynomial1DFree(psPolynomial1D* myPoly)
 {
     psFree(myPoly->coeff);
@@ -288,5 +288,5 @@
 }
 
-static void polynomial2DFree(psPolynomial2D * myPoly)
+static void polynomial2DFree(psPolynomial2D* myPoly)
 {
     int x = 0;
@@ -302,5 +302,5 @@
 }
 
-static void polynomial3DFree(psPolynomial3D * myPoly)
+static void polynomial3DFree(psPolynomial3D* myPoly)
 {
     int x = 0;
@@ -323,5 +323,5 @@
 }
 
-static void polynomial4DFree(psPolynomial4D * myPoly)
+static void polynomial4DFree(psPolynomial4D* myPoly)
 {
     int w = 0;
@@ -353,5 +353,5 @@
     Polynomial coefficients will be accessed in [w][x][y][z] fashion.
  *****************************************************************************/
-float psPolynomial1DEval(float x, const psPolynomial1D * myPoly)
+float psPolynomial1DEval(float x, const psPolynomial1D* myPoly)
 {
     int loop_x = 0;
@@ -380,5 +380,5 @@
 }
 
-float psPolynomial2DEval(float x, float y, const psPolynomial2D * myPoly)
+float psPolynomial2DEval(float x, float y, const psPolynomial2D* myPoly)
 {
     int loop_x = 0;
@@ -400,5 +400,5 @@
 }
 
-float psPolynomial3DEval(float x, float y, float z, const psPolynomial3D * myPoly)
+float psPolynomial3DEval(float x, float y, float z, const psPolynomial3D* myPoly)
 {
     int loop_x = 0;
@@ -426,5 +426,5 @@
 }
 
-float psPolynomial4DEval(float w, float x, float y, float z, const psPolynomial4D * myPoly)
+float psPolynomial4DEval(float w, float x, float y, float z, const psPolynomial4D* myPoly)
 {
     int loop_w = 0;
@@ -458,10 +458,10 @@
 }
 
-psDPolynomial1D *psDPolynomial1DAlloc(int n)
+psDPolynomial1D* psDPolynomial1DAlloc(int n)
 {
     int i = 0;
-    psDPolynomial1D *newPoly = NULL;
-
-    newPoly = (psDPolynomial1D *) psAlloc(sizeof(psDPolynomial1D));
+    psDPolynomial1D* newPoly = NULL;
+
+    newPoly = (psDPolynomial1D* ) psAlloc(sizeof(psDPolynomial1D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial1DFree);
     newPoly->n = n;
@@ -478,11 +478,11 @@
 }
 
-psDPolynomial2D *psDPolynomial2DAlloc(int nX, int nY)
-{
-    int x = 0;
-    int y = 0;
-    psDPolynomial2D *newPoly = NULL;
-
-    newPoly = (psDPolynomial2D *) psAlloc(sizeof(psDPolynomial2D));
+psDPolynomial2D* psDPolynomial2DAlloc(int nX, int nY)
+{
+    int x = 0;
+    int y = 0;
+    psDPolynomial2D* newPoly = NULL;
+
+    newPoly = (psDPolynomial2D* ) psAlloc(sizeof(psDPolynomial2D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial2DFree);
     newPoly->nX = nX;
@@ -508,12 +508,12 @@
 }
 
-psDPolynomial3D *psDPolynomial3DAlloc(int nX, int nY, int nZ)
+psDPolynomial3D* psDPolynomial3DAlloc(int nX, int nY, int nZ)
 {
     int x = 0;
     int y = 0;
     int z = 0;
-    psDPolynomial3D *newPoly = NULL;
-
-    newPoly = (psDPolynomial3D *) psAlloc(sizeof(psDPolynomial3D));
+    psDPolynomial3D* newPoly = NULL;
+
+    newPoly = (psDPolynomial3D* ) psAlloc(sizeof(psDPolynomial3D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial3DFree);
     newPoly->nX = nX;
@@ -547,5 +547,5 @@
 }
 
-psDPolynomial4D *psDPolynomial4DAlloc(int nW, int nX, int nY, int nZ)
+psDPolynomial4D* psDPolynomial4DAlloc(int nW, int nX, int nY, int nZ)
 {
     int w = 0;
@@ -553,7 +553,7 @@
     int y = 0;
     int z = 0;
-    psDPolynomial4D *newPoly = NULL;
-
-    newPoly = (psDPolynomial4D *) psAlloc(sizeof(psDPolynomial4D));
+    psDPolynomial4D* newPoly = NULL;
+
+    newPoly = (psDPolynomial4D* ) psAlloc(sizeof(psDPolynomial4D));
     p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial4DFree);
     newPoly->nW = nW;
@@ -595,5 +595,5 @@
 }
 
-static void dPolynomial1DFree(psDPolynomial1D * myPoly)
+static void dPolynomial1DFree(psDPolynomial1D* myPoly)
 {
     psFree(myPoly->coeff);
@@ -602,5 +602,5 @@
 }
 
-static void dPolynomial2DFree(psDPolynomial2D * myPoly)
+static void dPolynomial2DFree(psDPolynomial2D* myPoly)
 {
     int x = 0;
@@ -616,5 +616,5 @@
 }
 
-static void dPolynomial3DFree(psDPolynomial3D * myPoly)
+static void dPolynomial3DFree(psDPolynomial3D* myPoly)
 {
     int x = 0;
@@ -637,5 +637,5 @@
 }
 
-static void dPolynomial4DFree(psDPolynomial4D * myPoly)
+static void dPolynomial4DFree(psDPolynomial4D* myPoly)
 {
     int w = 0;
@@ -667,5 +667,5 @@
     Polynomial coefficients will be accessed in [w][x][y][z] fashion.
  *****************************************************************************/
-double psDPolynomial1DEval(double x, const psDPolynomial1D * myPoly)
+double psDPolynomial1DEval(double x, const psDPolynomial1D* myPoly)
 {
     int loop_x = 0;
@@ -686,5 +686,5 @@
 }
 
-double psDPolynomial2DEval(double x, double y, const psDPolynomial2D * myPoly)
+double psDPolynomial2DEval(double x, double y, const psDPolynomial2D* myPoly)
 {
     int loop_x = 0;
@@ -706,5 +706,5 @@
 }
 
-double psDPolynomial3DEval(double x, double y, double z, const psDPolynomial3D * myPoly)
+double psDPolynomial3DEval(double x, double y, double z, const psDPolynomial3D* myPoly)
 {
     int loop_x = 0;
@@ -732,5 +732,5 @@
 }
 
-double psDPolynomial4DEval(double w, double x, double y, double z, const psDPolynomial4D * myPoly)
+double psDPolynomial4DEval(double w, double x, double y, double z, const psDPolynomial4D* myPoly)
 {
     int loop_w = 0;