Index: trunk/psLib/test/math/tst_psMinimizeLMM.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimizeLMM.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psMinimizeLMM.c	(revision 6484)
@@ -61,4 +61,7 @@
     myCoords->n = NUM_DATA_POINTS;
     psVector *y = psVectorAlloc(NUM_DATA_POINTS, PS_TYPE_F32);
+    myParams->n = NUM_PARAMS;
+    myDerivs->n = NUM_PARAMS;
+    y->n = NUM_DATA_POINTS;
 
     for (psS32 i=0;i<NUM_DATA_POINTS;i++) {
@@ -66,4 +69,5 @@
         ((psVector *) (myCoords->data[i]))->data.F32[0] = (float) (i+10);
         ((psVector *) (myCoords->data[i]))->data.F32[1] = (float) (i+3);
+        ((psVector *) (myCoords->data[i]))->n = 2;
         psTrace(__func__, 5, "x[%d] is vector (%f, %f)\n", i, ((psVector *) (myCoords->data[i]))->data.F32[0],
                 ((psVector *) (myCoords->data[i]))->data.F32[1]);
Index: trunk/psLib/test/math/tst_psMinimizePowell.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimizePowell.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psMinimizePowell.c	(revision 6484)
@@ -81,4 +81,6 @@
     myParams = psVectorAlloc(NUM_PARAMS, PS_TYPE_F32);
     myParamMask = psVectorAlloc(NUM_PARAMS, PS_TYPE_U8);
+    myParams->n = NUM_PARAMS;
+    myParamMask->n = NUM_PARAMS;
     min = psMinimizationAlloc(100, 0.01);
 
@@ -164,4 +166,5 @@
         ((psVector *) (myCoords->data[i]))->data.F32[0] = (float) (i+10);
         ((psVector *) (myCoords->data[i]))->data.F32[1] = (float) (i+3);
+        ((psVector *) (myCoords->data[i]))->n = 2;
     }
     for (i=0;i<NUM_PARAMS;i++) {
@@ -169,4 +172,5 @@
         myParams->data.F32[i] = 0.0;
         myParams->data.F32[i] = (float) i;
+        myParams->n++;
     }
 
@@ -239,4 +243,5 @@
         ((psVector *) (myCoords->data[i]))->data.F32[0] = (float) (i+10);
         ((psVector *) (myCoords->data[i]))->data.F32[1] = (float) (i+3);
+        ((psVector *) (myCoords->data[i]))->n++;
     }
     for (i=0;i<NUM_PARAMS;i++) {
@@ -245,4 +250,6 @@
         myParams->data.F32[i] = (float) i;
         myParamMask->data.U8[i] = 0;
+        myParams->n++;
+        myParamMask->n++;
     }
 
Index: trunk/psLib/test/math/tst_psPolyFit1D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolyFit1D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolyFit1D.c	(revision 6484)
@@ -99,4 +99,5 @@
         printf(" using a psF32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F32[i] = (psF32) i;
@@ -111,4 +112,5 @@
         printf(" using a psS32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_S32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.S32[i] = (psS32) i;
@@ -123,4 +125,5 @@
         printf(" using a psF64 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F64);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F64[i] = (psF64) i;
@@ -139,4 +142,5 @@
         printf(" using a psF32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F32[i] = setData((psF32) i);
@@ -159,4 +163,5 @@
         printf(" using a psS32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_S32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.S32[i] = (psS32) setData((psF32) i);
@@ -179,4 +184,5 @@
         printf(" using a psF64 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F64);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F64[i] = (psF64) setData((psF32) i);
@@ -203,4 +209,5 @@
         printf(" using a psF32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F32[i] = YERR;
@@ -211,4 +218,5 @@
         printf(" using a psS32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_S32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.S32[i] = (psS32) YERR;
@@ -219,4 +227,5 @@
         printf(" using a psF64 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F64);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F64[i] = YERR;
@@ -231,4 +240,5 @@
         printf(" using a psU8 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_U8);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = 0;
@@ -239,4 +249,5 @@
         printf(" using a psS32 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_S32);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = 0;
Index: trunk/psLib/test/math/tst_psPolyFit2D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolyFit2D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolyFit2D.c	(revision 6484)
@@ -100,4 +100,5 @@
         printf(" using a psF32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F32[i] = (psF32) i;
@@ -108,4 +109,5 @@
         printf(" using a psS32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_S32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.S32[i] = (psS32) i;
@@ -116,4 +118,5 @@
         printf(" using a psF64 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F64);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F64[i] = (psF64) i;
@@ -129,4 +132,5 @@
         printf(" using a psF32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F32[i] = (psF32) i;
@@ -137,4 +141,5 @@
         printf(" using a psS32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_S32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.S32[i] = (psS32) i;
@@ -145,4 +150,5 @@
         printf(" using a psF64 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F64);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F64[i] = (psF64) i;
@@ -157,4 +163,5 @@
         printf(" using a psF32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F32[i] = setData((psF32) i, (psF32) i);
@@ -177,4 +184,5 @@
         printf(" using a psS32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_S32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.S32[i] = (psS32) setData((psF32) i, (psF32) i);
@@ -197,4 +205,5 @@
         printf(" using a psF64 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F64);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F64[i] = (psF64) setData((psF32) i, (psF32) i);
@@ -221,4 +230,5 @@
         printf(" using a psF32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F32[i] = YERR;
@@ -237,4 +247,5 @@
         printf(" using a psF64 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F64);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F64[i] = YERR;
@@ -249,4 +260,5 @@
         printf(" using a psU8 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_U8);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = 0;
@@ -257,4 +269,5 @@
         printf(" using a psS32 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_S32);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = 0;
Index: trunk/psLib/test/math/tst_psPolyFit3D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolyFit3D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolyFit3D.c	(revision 6484)
@@ -117,4 +117,5 @@
         printf(" using a psF32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F32[i] = (psF32) i;
@@ -125,4 +126,5 @@
         printf(" using a psS32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_S32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.S32[i] = (psS32) i;
@@ -133,4 +135,5 @@
         printf(" using a psF64 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F64);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F64[i] = (psF64) i;
@@ -146,4 +149,5 @@
         printf(" using a psF32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F32[i] = (psF32) i;
@@ -154,4 +158,5 @@
         printf(" using a psS32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_S32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.S32[i] = (psS32) i;
@@ -162,4 +167,5 @@
         printf(" using a psF64 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F64);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F64[i] = (psF64) i;
@@ -174,4 +180,5 @@
         printf(" using a psF32 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_F32);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.F32[i] = (psF32) i;
@@ -182,4 +189,5 @@
         printf(" using a psS32 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_S32);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.S32[i] = (psS32) i;
@@ -190,4 +198,5 @@
         printf(" using a psF64 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_F64);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.F64[i] = (psF64) i;
@@ -203,4 +212,5 @@
         printf(" using a psF32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F32[i] = setData((psF32) i, (psF32) i, (psF32) i);
@@ -223,4 +233,5 @@
         printf(" using a psS32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_S32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.S32[i] = (psS32) setData((psF32) i, (psF32) i, (psF32) i);
@@ -243,4 +254,5 @@
         printf(" using a psF64 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F64);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F64[i] = (psF64) setData((psF32) i, (psF32) i, (psF32) i);
@@ -267,4 +279,5 @@
         printf(" using a psF32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F32[i] = YERR;
@@ -275,4 +288,5 @@
         printf(" using a psS32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_S32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.S32[i] = (psS32) YERR;
@@ -283,4 +297,5 @@
         printf(" using a psF64 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F64);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F64[i] = YERR;
@@ -295,4 +310,5 @@
         printf(" using a psU8 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_U8);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = 0;
@@ -303,4 +319,5 @@
         printf(" using a psS32 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_S32);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = 0;
Index: trunk/psLib/test/math/tst_psPolyFit4D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolyFit4D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolyFit4D.c	(revision 6484)
@@ -130,4 +130,5 @@
         printf(" using a psF32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F32[i] = (psF32) i;
@@ -138,4 +139,5 @@
         printf(" using a psS32 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_S32);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.S32[i] = (psS32) i;
@@ -146,4 +148,5 @@
         printf(" using a psF64 x vector\n");
         x = psVectorAlloc(numData, PS_TYPE_F64);
+        x->n = numData;
         for (psS32 i=0;i<numData;i++) {
             x->data.F64[i] = (psF64) i;
@@ -159,4 +162,5 @@
         printf(" using a psF32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F32[i] = (psF32) i;
@@ -167,4 +171,5 @@
         printf(" using a psS32 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_S32);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.S32[i] = (psS32) i;
@@ -175,4 +180,5 @@
         printf(" using a psF64 y vector\n");
         y = psVectorAlloc(numData, PS_TYPE_F64);
+        y->n = numData;
         for (psS32 i=0;i<numData;i++) {
             y->data.F64[i] = (psF64) i;
@@ -187,4 +193,5 @@
         printf(" using a psF32 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_F32);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.F32[i] = (psF32) i;
@@ -195,4 +202,5 @@
         printf(" using a psS32 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_S32);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.S32[i] = (psS32) i;
@@ -203,4 +211,5 @@
         printf(" using a psF64 z vector\n");
         z = psVectorAlloc(numData, PS_TYPE_F64);
+        z->n = numData;
         for (psS32 i=0;i<numData;i++) {
             z->data.F64[i] = (psF64) i;
@@ -215,4 +224,5 @@
         printf(" using a psF32 t vector\n");
         t = psVectorAlloc(numData, PS_TYPE_F32);
+        t->n = numData;
         for (psS32 i=0;i<numData;i++) {
             t->data.F32[i] = (psF32) i;
@@ -223,4 +233,5 @@
         printf(" using a psS32 t vector\n");
         t = psVectorAlloc(numData, PS_TYPE_S32);
+        t->n = numData;
         for (psS32 i=0;i<numData;i++) {
             t->data.S32[i] = (psS32) i;
@@ -231,4 +242,5 @@
         printf(" using a psF64 t vector\n");
         t = psVectorAlloc(numData, PS_TYPE_F64);
+        t->n = numData;
         for (psS32 i=0;i<numData;i++) {
             t->data.F64[i] = (psF64) i;
@@ -244,4 +256,5 @@
         printf(" using a psF32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F32[i] = setData((psF32) i, (psF32) i, (psF32) i, (psF32) i);
@@ -264,4 +277,5 @@
         printf(" using a psS32 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_S32);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.S32[i] = (psS32) setData((psF32) i, (psF32) i, (psF32) i, (psF32) i);
@@ -284,4 +298,5 @@
         printf(" using a psF64 f vector\n");
         f = psVectorAlloc(numData, PS_TYPE_F64);
+        f->n = numData;
         for (psS32 i=0;i<numData;i++) {
             f->data.F64[i] = (psF64) setData((psF32) i, (psF32) i, (psF32) i, (psF32) i);
@@ -308,4 +323,5 @@
         printf(" using a psF32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F32[i] = YERR;
@@ -316,4 +332,5 @@
         printf(" using a psS32 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_S32);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.S32[i] = (psS32) YERR;
@@ -324,4 +341,5 @@
         printf(" using a psF64 fErr vector\n");
         fErr = psVectorAlloc(numData, PS_TYPE_F64);
+        fErr->n = numData;
         for (psS32 i=0;i<numData;i++) {
             fErr->data.F64[i] = YERR;
@@ -336,4 +354,5 @@
         printf(" using a psU8 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_U8);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = 0;
@@ -344,4 +363,5 @@
         printf(" using a psS32 mask vector\n");
         mask = psVectorAlloc(numData, PS_TYPE_S32);
+        mask->n = numData;
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = 0;
Index: trunk/psLib/test/math/tst_psPolynomialEval1D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolynomialEval1D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolynomialEval1D.c	(revision 6484)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.1 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2006-02-02 21:05:51 $
+*  @version  $Revision: 1.2 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2006-02-24 23:43:15 $
 *
 *  XXX: Probably should test single- and multi-dimensional polynomials in
@@ -11,5 +11,5 @@
 *
 *  XXX: define ORDERS, not TERMS
-* 
+*
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -124,4 +124,6 @@
         inputOrd->data.F64[i] = poly1DXValue[i];
         inputCheb->data.F64[i] = poly1DXChebValue[i];
+        inputOrd->n++;
+        inputCheb->n++;
     }
 
Index: trunk/psLib/test/math/tst_psPolynomialEval2D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolynomialEval2D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolynomialEval2D.c	(revision 6484)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.1 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2006-02-02 21:05:51 $
+*  @version  $Revision: 1.2 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2006-02-24 23:43:15 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -157,4 +157,8 @@
         inputChebX->data.F64[i] = poly2DXYChebValue[i][0];
         inputChebY->data.F64[i] = poly2DXYChebValue[i][1];
+        inputOrdX->n++;
+        inputOrdY->n++;
+        inputChebX->n++;
+        inputChebY->n++;
     }
 
Index: trunk/psLib/test/math/tst_psPolynomialEval3D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolynomialEval3D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolynomialEval3D.c	(revision 6484)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.1 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2006-02-02 21:05:51 $
+*  @version  $Revision: 1.2 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2006-02-24 23:43:15 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -220,4 +220,10 @@
         inputChebY->data.F64[i] = Dpoly3DXYZChebValue[i][1];
         inputChebZ->data.F64[i] = Dpoly3DXYZChebValue[i][2];
+        inputOrdX->n++;
+        inputOrdY->n++;
+        inputOrdZ->n++;
+        inputChebX->n++;
+        inputChebY->n++;
+        inputChebZ->n++;
     }
 
Index: trunk/psLib/test/math/tst_psPolynomialEval4D.c
===================================================================
--- trunk/psLib/test/math/tst_psPolynomialEval4D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psPolynomialEval4D.c	(revision 6484)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.1 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2006-02-02 21:05:51 $
+*  @version  $Revision: 1.2 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2006-02-24 23:43:15 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -470,4 +470,12 @@
         inputChebY->data.F64[i] = Dpoly4DWXYZChebValue[i][2];
         inputChebZ->data.F64[i] = Dpoly4DWXYZChebValue[i][3];
+        inputOrdW->n++;
+        inputOrdX->n++;
+        inputOrdY->n++;
+        inputOrdZ->n++;
+        inputChebW->n++;
+        inputChebX->n++;
+        inputChebY->n++;
+        inputChebZ->n++;
     }
 
Index: trunk/psLib/test/math/tst_psRandom.c
===================================================================
--- trunk/psLib/test/math/tst_psRandom.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psRandom.c	(revision 6484)
@@ -122,4 +122,5 @@
     psRandom *myRNG = NULL;
     psVector *rans = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans->n = rans->nalloc;
     psStats *stats = psStatsAlloc(PS_STAT_SAMPLE_MEAN);
 
@@ -167,4 +168,5 @@
     psRandom *myRNG = NULL;
     psVector *rans = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans->n = rans->nalloc;
     psStats *stats = psStatsAlloc(PS_STAT_SAMPLE_MEAN);
 
@@ -214,4 +216,5 @@
     psRandom *myRNG = NULL;
     psVector *rans = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans->n = rans->nalloc;
     psStats *stats = psStatsAlloc(PS_STAT_SAMPLE_MEAN);
 
@@ -259,6 +262,9 @@
     psRandom *myRNG = NULL;
     psVector *rans00 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans00->n = rans00->nalloc;
     psVector *rans01 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans01->n = rans01->nalloc;
     psVector *rans02 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans02->n = rans02->nalloc;
 
     myRNG = psRandomAlloc(PS_RANDOM_TAUS, SEED);
@@ -314,6 +320,9 @@
     psRandom *myRNG = NULL;
     psVector *rans00 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans00->n = rans00->nalloc;
     psVector *rans01 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans01->n = rans01->nalloc;
     psVector *rans02 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans02->n = rans02->nalloc;
 
     myRNG = psRandomAlloc(PS_RANDOM_TAUS, SEED);
@@ -357,6 +366,9 @@
     psRandom *myRNG = NULL;
     psVector *rans00 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans00->n = rans00->nalloc;
     psVector *rans01 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans01->n = rans01->nalloc;
     psVector *rans02 = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
+    rans02->n = rans02->nalloc;
 
     myRNG = psRandomAlloc(PS_RANDOM_TAUS, SEED);
Index: trunk/psLib/test/math/tst_psSpline1D.c
===================================================================
--- trunk/psLib/test/math/tst_psSpline1D.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psSpline1D.c	(revision 6484)
@@ -13,6 +13,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2006-01-27 20:08:58 $
+*  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2006-02-24 23:43:15 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -28,7 +28,7 @@
 // {testFunction, testpointNumber, description, expected rc, boolean-ignore this test}
 testDescription tests[] = {
-                              {psSplineAllocTest,      0000, "(TEST A) psSplineAllocTest",      true, false},
-                              {psSplineEvalTest,       0000, "(TEST B) psSplineEvalTest",       true, false},
-                              {psSplineEvalVectorTest, 0000, "(TEST C) psSplineEvalVectorTest", true, false},
+                              {psSplineAllocTest, 665, "(TEST A) psSplineAllocTest",true, false},
+                              {psSplineEvalTest, 666, "(TEST B) psSplineEvalTest",true, false},
+                              {psSplineEvalVectorTest,667,"(TEST C) psSplineEvalVectorTest",true,false},
                               {NULL}
                           };
@@ -199,4 +199,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -274,4 +278,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -384,4 +392,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -459,4 +471,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -559,4 +575,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -634,4 +654,8 @@
     yF32 = psVectorAlloc(NumSplines+1, PS_TYPE_F32);
     yF64 = psVectorAlloc(NumSplines+1, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    xF64->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    yF64->n = NumSplines+1;
 
     for (psS32 i=0;i<NumSplines+1;i++) {
@@ -757,4 +781,9 @@
     psVector *yF32Test = NULL;
     psVector *xF64Test = psVectorAlloc(NumSplines, PS_TYPE_F64);
+    xF32->n = NumSplines+1;
+    yF32->n = NumSplines+1;
+    xF32Test->n = NumSplines;
+    xF64Test->n = NumSplines;
+
 
     for (psS32 i=0;i<NumSplines+1;i++) {
Index: trunk/psLib/test/math/tst_psStats07.c
===================================================================
--- trunk/psLib/test/math/tst_psStats07.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psStats07.c	(revision 6484)
@@ -72,4 +72,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.F32[i] = gaussVector->data.F32[i];
+            in->n++;
         }
 
@@ -86,4 +87,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.F64[i] = (psF64) gaussVector->data.F32[i];
+            in->n++;
         }
 
@@ -100,4 +102,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.S8[i] = (psS8) gaussVector->data.F32[i];
+            in->n++;
         }
 
@@ -114,4 +117,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.U16[i] = (psU16) gaussVector->data.F32[i];
+            in->n++;
         }
 
@@ -128,4 +132,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.S32[i] = (psS32) gaussVector->data.F32[i];
+            in->n++;
         }
 
@@ -138,4 +143,5 @@
     psFree(gaussVector);
 
+    //    in->n = in->nalloc;
     if (flags & TST_ERRORS_NULL) {
         printf("        using a NULL errors vector\n");
@@ -147,4 +153,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.F32[i] = ERRORS;
+            errors->n++;
         }
 
@@ -161,4 +168,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.F64[i] = ERRORS;
+            errors->n++;
         }
 
@@ -175,4 +183,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.S8[i] = (psS8) ERRORS;
+            errors->n++;
         }
 
@@ -189,4 +198,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.U16[i] = (psU16) ERRORS;
+            errors->n++;
         }
 
@@ -203,4 +213,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.S32[i] = (psS32) ERRORS;
+            errors->n++;
         }
 
@@ -222,4 +233,5 @@
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = (psU8) 0;
+            mask->n++;
         }
 
@@ -236,4 +248,5 @@
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = (psS32) 0;
+            mask->n++;
         }
 
Index: trunk/psLib/test/math/tst_psStats09.c
===================================================================
--- trunk/psLib/test/math/tst_psStats09.c	(revision 6449)
+++ trunk/psLib/test/math/tst_psStats09.c	(revision 6484)
@@ -68,4 +68,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.F32[i] = (psF32) i;
+            in->n++;
         }
 
@@ -82,4 +83,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.F64[i] = (psF64) i;
+            in->n++;
         }
 
@@ -96,4 +98,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.S8[i] = (psS8) i;
+            in->n++;
         }
 
@@ -110,4 +113,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.U16[i] = (psU16) i;
+            in->n++;
         }
 
@@ -124,4 +128,5 @@
         for (psS32 i=0;i<numData;i++) {
             in->data.S32[i] = (psS32) i;
+            in->n++;
         }
 
@@ -142,4 +147,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.F32[i] = ERRORS;
+            errors->n++;
         }
 
@@ -156,4 +162,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.F64[i] = ERRORS;
+            errors->n++;
         }
 
@@ -170,4 +177,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.S8[i] = (psS8) ERRORS;
+            errors->n++;
         }
 
@@ -184,4 +192,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.U16[i] = (psU16) ERRORS;
+            errors->n++;
         }
 
@@ -198,4 +207,5 @@
         for (psS32 i=0;i<numData;i++) {
             errors->data.S32[i] = (psS32) ERRORS;
+            errors->n++;
         }
 
@@ -217,4 +227,5 @@
         for (psS32 i=0;i<numData;i++) {
             mask->data.U8[i] = (psU8) 0;
+            mask->n++;
         }
 
@@ -231,4 +242,5 @@
         for (psS32 i=0;i<numData;i++) {
             mask->data.S32[i] = (psS32) 0;
+            mask->n++;
         }