Index: trunk/psModules/src/imsubtract/pmImageSubtract.c
===================================================================
--- trunk/psModules/src/imsubtract/pmImageSubtract.c	(revision 6506)
+++ trunk/psModules/src/imsubtract/pmImageSubtract.c	(revision 6511)
@@ -7,6 +7,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-28 20:43:52 $
+ *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-04 01:01:33 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -119,7 +119,11 @@
     tmp->u = psVectorAlloc(nBF, PS_TYPE_F32);
     tmp->v = psVectorAlloc(nBF, PS_TYPE_F32);
+    tmp->u->n = tmp->u->nalloc;
+    tmp->v->n = tmp->v->nalloc;
     tmp->sigma = NULL;
     tmp->xOrder = psVectorAlloc(nBF, PS_TYPE_F32);
+    tmp->xOrder->n = tmp->xOrder->nalloc;
     tmp->yOrder = psVectorAlloc(nBF, PS_TYPE_F32);
+    tmp->yOrder->n = tmp->yOrder->nalloc;
     tmp->subIndex = 0;
     tmp->preCalc = NULL;
@@ -216,4 +220,6 @@
     tmp->u = psVectorAlloc(nBF, PS_TYPE_F32);
     tmp->v = psVectorAlloc(nBF, PS_TYPE_F32);
+    tmp->u->n = tmp->u->nalloc;
+    tmp->v->n = tmp->v->nalloc;
     tmp->sigma = psVectorAlloc(nBF, PS_TYPE_F32);
     tmp->xOrder = psVectorAlloc(nBF, PS_TYPE_F32);
@@ -223,4 +229,10 @@
     tmp->spatialOrder = spatialOrder;
     tmp->preCalc = psArrayAlloc(nBF);
+    tmp->u->n = tmp->u->nalloc;
+    tmp->v->n = tmp->v->nalloc;
+    tmp->sigma->n = tmp->sigma->nalloc;
+    tmp->xOrder->n = tmp->xOrder->nalloc;
+    tmp->yOrder->n = tmp->yOrder->nalloc;
+    tmp->preCalc->n = tmp->preCalc->nalloc;
 
     //
@@ -352,4 +364,5 @@
     } else {
         stamps = (psArray *) psArrayAlloc(xNum * yNum);
+        stamps->n = stamps->nalloc;
         for (psS32 s = 0 ; s < (xNum * yNum) ; s++) {
             stamps->data[s] = (psPtr *) p_pmStampAlloc(PM_STAMP_REJECTED);
@@ -810,4 +823,5 @@
             if (stampVector == NULL) {
                 stampVector = psVectorAlloc(numSolveParams, PS_TYPE_F64);
+                stampVector->n = stampVector->nalloc;
                 stamp->vector = stampVector;
             } else {
@@ -1028,4 +1042,5 @@
     } else {
         solution = psVectorAlloc(size, PS_TYPE_F64);
+        solution->n = solution->nalloc;
     }
 
@@ -1039,4 +1054,5 @@
     psImage *sumMatrix = psImageAlloc(size, size, PS_TYPE_F64);
     psVector *sumVector = psVectorAlloc(size, PS_TYPE_F64);
+    sumVector->n = sumVector->nalloc;
     PS_VECTOR_SET_F64(sumVector, 0.0);
     PS_IMAGE_SET_F64(sumMatrix, 0.0);
@@ -1131,4 +1147,5 @@
     } else {
         deviations = psVectorAlloc(stamps->n, PS_TYPE_F32);
+        deviations->n = deviations->nalloc;
         // XXX: Probably not necessary.
         PS_VECTOR_SET_F32(deviations, 0.0);
Index: trunk/psModules/src/imsubtract/pmSubtractBias.c
===================================================================
--- trunk/psModules/src/imsubtract/pmSubtractBias.c	(revision 6506)
+++ trunk/psModules/src/imsubtract/pmSubtractBias.c	(revision 6511)
@@ -6,6 +6,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.10 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-03-02 23:16:05 $
+ *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-04 01:01:33 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -307,4 +307,5 @@
     // XXX: use psVectorCopy if they get it working.
     out->knots = psVectorAlloc(in->knots->n, in->knots->type.type);
+    out->knots->n = out->knots->nalloc;
     for (psS32 i = 0 ; i < in->knots->n ; i++) {
         out->knots->data.F32[i] = in->knots->data.F32[i];
@@ -346,5 +347,5 @@
     // Allocate the new vector.
     psVector *newVec = psVectorAlloc(n, PS_TYPE_F32);
-
+    newVec->n = newVec->nalloc;
     //
     // If the new vector is the same size as the old, simply copy the data.
@@ -491,4 +492,5 @@
     psS32 numOverscanImages = my_psListLength(bias);
     psVector *statsAll = psVectorAlloc(numOverscanImages, PS_TYPE_F32);
+    statsAll->n = statsAll->nalloc;
     psListElem *tmpOverscan = (psListElem *) bias->head;
     psS32 i = 0;
@@ -554,4 +556,6 @@
     psVector *tmpRow = psVectorAlloc(overscanImage->numCols, PS_TYPE_F32);
     psVector *tmpCol = psVectorAlloc(overscanImage->numRows, PS_TYPE_F32);
+    tmpRow->n = tmpRow->nalloc;
+    tmpCol->n = tmpCol->nalloc;
 
     //
@@ -594,4 +598,6 @@
     psVector *tmpRow = psVectorAlloc(overscanImage->numCols, PS_TYPE_F32);
     psVector *tmpCol = psVectorAlloc(overscanImage->numRows, PS_TYPE_F32);
+    tmpRow->n = tmpRow->nalloc;
+    tmpCol->n = tmpCol->nalloc;
 
     //
@@ -649,4 +655,5 @@
     psS32 vecSize = GetOverscanSize(inImg, overScanAxis);
     psVector *overscanVector = psVectorAlloc(vecSize, PS_TYPE_F32);
+    overscanVector->n = overscanVector->nalloc;
 
     //
@@ -658,4 +665,5 @@
     psS32 numOverscanImages = my_psListLength(bias);
     psVector **overscanVectors = (psVector **) psAlloc(numOverscanImages * sizeof(psVector *));
+    //    (*overscanVectors)->n = (*overscanVectors)->nalloc;
     for (psS32 i = 0 ; i < numOverscanImages ; i++) {
         overscanVectors[i] = NULL;
@@ -711,4 +719,5 @@
     //
     psVector *tmpVec = psVectorAlloc(numOverscanImages, PS_TYPE_F32);
+    tmpVec->n = tmpVec->nalloc;
     psF64 statValue;
     for (psS32 i = 0 ; i < vecSize ; i++) {
@@ -772,4 +781,6 @@
         psVector *myBin = psVectorAlloc(numBins, PS_TYPE_F32);
         psVector *binVec = psVectorAlloc(nBinOrig, PS_TYPE_F32);
+        myBin->n = myBin->nalloc;
+        binVec->n = binVec->nalloc;
 
         for (psS32 i=0;i<numBins;i++) {
@@ -859,4 +870,5 @@
     psS32 vecSize = GetOverscanSize(inImg, overScanAxis);
     psVector *newVec = psVectorAlloc(vecSize, PS_TYPE_F32);
+    newVec->n = newVec->nalloc;
     if ((nBin > 1) && (nBin < overscanVector->n)) {
         for (psS32 i = 0 ; i < vecSize ; i++) {
@@ -905,4 +917,5 @@
 
     psVector *newVec = psVectorAlloc(vecSize, PS_TYPE_F32);
+    newVec->n = newVec->nalloc;
     for (psS32 i = 0 ; i < vecSize ; i++) {
         newVec->data.F32[i] = overscanVector->data.F32[i/nBin];
Index: trunk/psModules/src/imsubtract/pmSubtractSky.c
===================================================================
--- trunk/psModules/src/imsubtract/pmSubtractSky.c	(revision 6506)
+++ trunk/psModules/src/imsubtract/pmSubtractSky.c	(revision 6511)
@@ -6,10 +6,10 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-02-06 21:03:25 $
+ *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-04 01:01:33 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
  *
- *      
+ *
  *
  */
@@ -341,4 +341,5 @@
     psImage *Aout = psImageAlloc(localPolyTerms, localPolyTerms, PS_TYPE_F64);
     psVector *B = psVectorAlloc(localPolyTerms, PS_TYPE_F64);
+    B->n = B->nalloc;
     psVector *outPerm = NULL;
 
@@ -407,4 +408,5 @@
     PS_ASSERT_IMAGE_NON_EMPTY(Aout, NULL);
     psVector *C = psVectorAlloc(localPolyTerms, PS_TYPE_F64);
+    C->n = C->nalloc;
     psMatrixLUSolve(C, Aout, B, outPerm);