Changeset 5102 for trunk/psLib/src/math/psSpline.c

Timestamp:

Sep 22, 2005, 3:56:54 PM (21 years ago)

Author:

gusciora

Message:

....

File:

• trunk/psLib/src/math/psSpline.c (modified) (6 diffs)

trunk/psLib/src/math/psSpline.c

@@ -7 +7 @@
 *  splines.
 *
-*  @version $Revision: 1.127 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-09-22 02:47:16 $
+*  @version $Revision: 1.128 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-09-23 01:56:54 $
 *
 *
@@ -207 +207 @@
 }
 
-
 /*****************************************************************************
 CalculateSecondDerivs(): Given a set of x/y vectors corresponding to a
@@ -223 +222 @@
 XXX: do an F64 version?
  *****************************************************************************/
-static psF32 *calculateSecondDerivs(const psVector* x,        ///< Ordinates (or NULL to just use the indices)
+static psF32 *calculateSecondDerivs(const psVector* x,        ///< Ordinates (or NULL to just use the in!
                                     const psVector* y)        ///< Coordinates
 {
@@ -229 +228 @@
             "---- calculateSecondDerivs() begin ----\n");
 
-    unsigned int i;
-    unsigned int k;
+    psS32 i;
+    psS32 k;
     psF32 sig;
     psF32 p;
-    unsigned int n = y->n;
+    psS32 n = y->n;
     psF32 *u = (psF32 *) psAlloc(n * sizeof(psF32));
     psF32 *derivs2 = (psF32 *) psAlloc(n * sizeof(psF32));
@@ -253 +252 @@
 
         psTrace(".psLib.dataManip.calculateSecondDerivs", 6,
-                "X[%u] is %f\n", i, X[i]);
+                "X[%d] is %f\n", i, X[i]);
         psTrace(".psLib.dataManip.calculateSecondDerivs", 6,
-                "Y[%u] is %f\n", i, Y[i]);
+                "Y[%d] is %f\n", i, Y[i]);
         psTrace(".psLib.dataManip.calculateSecondDerivs", 6,
-                "u[%u] is %f\n", i, u[i]);
+                "u[%d] is %f\n", i, u[i]);
     }
 
@@ -264 +263 @@
     derivs2[n-1] = (u[n-1] - (qn * u[n-2])) / ((qn * derivs2[n-2]) + 1.0);
 
-    //XXX: Changed from k=n-2;k>=0; because k is now unsigned//
-    for (k=(n-1);k>0;k--) {
+    for (k=(n-2);k>=0;k--) {
         derivs2[k] = derivs2[k] * derivs2[k+1] + u[k];
 
         psTrace(".psLib.dataManip.calculateSecondDerivs", 6,
-                "derivs2[%u] is %f\n", k, derivs2[k]);
-    }
-
+                "derivs2[%d] is %f\n", k, derivs2[k]);
+    }
     psFree(u);
     psTrace(".psLib.dataManip.calculateSecondDerivs", 4,
             "---- calculateSecondDerivs() end ----\n");
+
     return(derivs2);
 }
+
+
 
 /*****************************************************************************/