Index: /trunk/psLib/src/dataManip/psFunctions.c
===================================================================
--- /trunk/psLib/src/dataManip/psFunctions.c	(revision 3096)
+++ /trunk/psLib/src/dataManip/psFunctions.c	(revision 3097)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.84 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-01-26 20:41:04 $
+ *  @version $Revision: 1.85 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-26 20:47:41 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -1841,4 +1841,6 @@
  
 XXX: Ensure that domain[i+1] != domain[i]
+ 
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAlloc(psS32 numSplines,
@@ -1884,4 +1886,5 @@
 
 /*****************************************************************************
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAllocGeneric(const psVector *bounds,
@@ -2142,4 +2145,7 @@
 XXX: The spline eval functions require input and output to be F32.  however
      the spline fit functions require F32 and F64.
+ 
+XXX: This only works if spline0>knots if psF32.  Must add support for psU32 and
+psF64.
  *****************************************************************************/
 psF32 psSpline1DEval(
@@ -2150,4 +2156,5 @@
     PS_PTR_CHECK_NULL(spline, NAN);
     PS_INT_CHECK_NON_NEGATIVE(spline->n, NAN);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NAN);
 
     psS32 binNum;
@@ -2186,4 +2193,5 @@
     PS_VECTOR_CHECK_NULL(x, NULL);
     PS_VECTOR_CHECK_TYPE_F32_OR_F64(x, NULL);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NULL);
 
     psS32 i;
Index: /trunk/psLib/src/dataManip/psMinimize.c
===================================================================
--- /trunk/psLib/src/dataManip/psMinimize.c	(revision 3096)
+++ /trunk/psLib/src/dataManip/psMinimize.c	(revision 3097)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.102 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-01-26 20:41:04 $
+ *  @version $Revision: 1.103 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-26 20:47:41 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -162,4 +162,6 @@
  
 XXX: This is F32 only
+ 
+XXX: spline->knots must be psF32
  *****************************************************************************/
 /*
@@ -177,4 +179,5 @@
     PS_VECTOR_CHECK_NULL(y, NAN);
     PS_VECTOR_CHECK_TYPE(y, PS_TYPE_F32, NAN);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NULL);
  
     psS32 n;
@@ -244,4 +247,6 @@
 XXX: What happens if X is NULL, then an index vector is generated for X, but
 that index vector lies outside the range vectors in mySpline?
+ 
+XXX: Assumes mySpline->knots is psF32.  Must add psU32 and psF64.
  *****************************************************************************/
 psSpline1D *psVectorFitSpline1D(psSpline1D *mySpline,              ///< The spline which will be generated.
@@ -252,4 +257,5 @@
     PS_VECTOR_CHECK_NULL(y, NULL);
     PS_VECTOR_CHECK_TYPE_F32_OR_F64(y, NULL);
+    PS_VECTOR_CHECK_TYPE(mySpline->knots, PS_TYPE_F32, NULL);
 
     psTrace(".psLib.dataManip.psVectorFitSpline1D", 4,
Index: /trunk/psLib/src/math/psMinimize.c
===================================================================
--- /trunk/psLib/src/math/psMinimize.c	(revision 3096)
+++ /trunk/psLib/src/math/psMinimize.c	(revision 3097)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.102 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-01-26 20:41:04 $
+ *  @version $Revision: 1.103 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-26 20:47:41 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -162,4 +162,6 @@
  
 XXX: This is F32 only
+ 
+XXX: spline->knots must be psF32
  *****************************************************************************/
 /*
@@ -177,4 +179,5 @@
     PS_VECTOR_CHECK_NULL(y, NAN);
     PS_VECTOR_CHECK_TYPE(y, PS_TYPE_F32, NAN);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NULL);
  
     psS32 n;
@@ -244,4 +247,6 @@
 XXX: What happens if X is NULL, then an index vector is generated for X, but
 that index vector lies outside the range vectors in mySpline?
+ 
+XXX: Assumes mySpline->knots is psF32.  Must add psU32 and psF64.
  *****************************************************************************/
 psSpline1D *psVectorFitSpline1D(psSpline1D *mySpline,              ///< The spline which will be generated.
@@ -252,4 +257,5 @@
     PS_VECTOR_CHECK_NULL(y, NULL);
     PS_VECTOR_CHECK_TYPE_F32_OR_F64(y, NULL);
+    PS_VECTOR_CHECK_TYPE(mySpline->knots, PS_TYPE_F32, NULL);
 
     psTrace(".psLib.dataManip.psVectorFitSpline1D", 4,
Index: /trunk/psLib/src/math/psPolynomial.c
===================================================================
--- /trunk/psLib/src/math/psPolynomial.c	(revision 3096)
+++ /trunk/psLib/src/math/psPolynomial.c	(revision 3097)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.84 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-01-26 20:41:04 $
+ *  @version $Revision: 1.85 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-26 20:47:41 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -1841,4 +1841,6 @@
  
 XXX: Ensure that domain[i+1] != domain[i]
+ 
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAlloc(psS32 numSplines,
@@ -1884,4 +1886,5 @@
 
 /*****************************************************************************
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAllocGeneric(const psVector *bounds,
@@ -2142,4 +2145,7 @@
 XXX: The spline eval functions require input and output to be F32.  however
      the spline fit functions require F32 and F64.
+ 
+XXX: This only works if spline0>knots if psF32.  Must add support for psU32 and
+psF64.
  *****************************************************************************/
 psF32 psSpline1DEval(
@@ -2150,4 +2156,5 @@
     PS_PTR_CHECK_NULL(spline, NAN);
     PS_INT_CHECK_NON_NEGATIVE(spline->n, NAN);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NAN);
 
     psS32 binNum;
@@ -2186,4 +2193,5 @@
     PS_VECTOR_CHECK_NULL(x, NULL);
     PS_VECTOR_CHECK_TYPE_F32_OR_F64(x, NULL);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NULL);
 
     psS32 i;
Index: /trunk/psLib/src/math/psSpline.c
===================================================================
--- /trunk/psLib/src/math/psSpline.c	(revision 3096)
+++ /trunk/psLib/src/math/psSpline.c	(revision 3097)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.84 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-01-26 20:41:04 $
+ *  @version $Revision: 1.85 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-26 20:47:41 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -1841,4 +1841,6 @@
  
 XXX: Ensure that domain[i+1] != domain[i]
+ 
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAlloc(psS32 numSplines,
@@ -1884,4 +1886,5 @@
 
 /*****************************************************************************
+XXX: What should be the defualty type for knots be?  psF32 is assumed.
  *****************************************************************************/
 psSpline1D *psSpline1DAllocGeneric(const psVector *bounds,
@@ -2142,4 +2145,7 @@
 XXX: The spline eval functions require input and output to be F32.  however
      the spline fit functions require F32 and F64.
+ 
+XXX: This only works if spline0>knots if psF32.  Must add support for psU32 and
+psF64.
  *****************************************************************************/
 psF32 psSpline1DEval(
@@ -2150,4 +2156,5 @@
     PS_PTR_CHECK_NULL(spline, NAN);
     PS_INT_CHECK_NON_NEGATIVE(spline->n, NAN);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NAN);
 
     psS32 binNum;
@@ -2186,4 +2193,5 @@
     PS_VECTOR_CHECK_NULL(x, NULL);
     PS_VECTOR_CHECK_TYPE_F32_OR_F64(x, NULL);
+    PS_VECTOR_CHECK_TYPE(spline->knots, PS_TYPE_F32, NULL);
 
     psS32 i;
