Index: trunk/psModules/src/objects/models/pmModel_GAUSS.c
===================================================================
--- trunk/psModules/src/objects/models/pmModel_GAUSS.c	(revision 6329)
+++ trunk/psModules/src/objects/models/pmModel_GAUSS.c	(revision 6511)
@@ -71,4 +71,7 @@
     *params_min = psVectorAlloc (7, PS_TYPE_F32);
     *params_max = psVectorAlloc (7, PS_TYPE_F32);
+    (*beta_lim)->n = (*beta_lim)->nalloc;
+    (*params_min)->n = (*params_min)->nalloc;
+    (*params_max)->n = (*params_max)->nalloc;
 
     beta_lim[0][0].data.F32[0] = 1000;
Index: trunk/psModules/src/objects/models/pmModel_PGAUSS.c
===================================================================
--- trunk/psModules/src/objects/models/pmModel_PGAUSS.c	(revision 6329)
+++ trunk/psModules/src/objects/models/pmModel_PGAUSS.c	(revision 6511)
@@ -45,4 +45,7 @@
     *params_min = psVectorAlloc (7, PS_TYPE_F32);
     *params_max = psVectorAlloc (7, PS_TYPE_F32);
+    (*beta_lim)->n = (*beta_lim)->nalloc;
+    (*params_min)->n = (*params_min)->nalloc;
+    (*params_max)->n = (*params_max)->nalloc;
 
     beta_lim[0][0].data.F32[0] = 1000;
Index: trunk/psModules/src/objects/models/pmModel_QGAUSS.c
===================================================================
--- trunk/psModules/src/objects/models/pmModel_QGAUSS.c	(revision 6329)
+++ trunk/psModules/src/objects/models/pmModel_QGAUSS.c	(revision 6511)
@@ -11,6 +11,6 @@
     params->data.F32[5] = sqrt(2.0) / SigmaY;
     params->data.F32[6] = Sxy;
-    params->data.F32[7] = 
-    params->data.F32[8] = 
+    params->data.F32[7] =
+    params->data.F32[8] =
 *****************************************************************************/
 
@@ -53,4 +53,7 @@
     *params_min = psVectorAlloc (8, PS_TYPE_F32);
     *params_max = psVectorAlloc (8, PS_TYPE_F32);
+    (*beta_lim)->n = (*beta_lim)->nalloc;
+    (*params_min)->n = (*params_min)->nalloc;
+    (*params_max)->n = (*params_max)->nalloc;
 
     beta_lim[0][0].data.F32[0] = 1000;
Index: trunk/psModules/src/objects/models/pmModel_SGAUSS.c
===================================================================
--- trunk/psModules/src/objects/models/pmModel_SGAUSS.c	(revision 6329)
+++ trunk/psModules/src/objects/models/pmModel_SGAUSS.c	(revision 6511)
@@ -64,4 +64,7 @@
     *params_min = psVectorAlloc (9, PS_TYPE_F32);
     *params_max = psVectorAlloc (9, PS_TYPE_F32);
+    (*beta_lim)->n = (*beta_lim)->nalloc;
+    (*params_min)->n = (*params_min)->nalloc;
+    (*params_max)->n = (*params_max)->nalloc;
 
     beta_lim[0][0].data.F32[0] = 1000;
@@ -116,4 +119,5 @@
 
     psVector *contour = psVectorAlloc (Nt, PS_TYPE_F32);
+    contour->n = contour->nalloc;
     for (t = 0, N = 0; (t < 2*M_PI) && (N < Nt); t += dt) {
         x = ct*axes.major*cos(t) + st*axes.minor*sin(t) + xo;
Index: trunk/psModules/src/objects/pmObjects.c
===================================================================
--- trunk/psModules/src/objects/pmObjects.c	(revision 6329)
+++ trunk/psModules/src/objects/pmObjects.c	(revision 6511)
@@ -6,6 +6,6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.9 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-02-06 22:17:54 $
+ *  @version $Revision: 1.10 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-04 01:01:33 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -96,4 +96,5 @@
 
     psVector *tmpVector = psVectorAlloc(image->numCols, PS_TYPE_F32);
+    tmpVector->n = tmpVector->nalloc;
     for (psU32 col = 0; col < image->numCols ; col++) {
         tmpVector->data.F32[col] = image->data.F32[row][col];
@@ -282,4 +283,6 @@
     tmp->params  = psVectorAlloc(Nparams, PS_TYPE_F32);
     tmp->dparams = psVectorAlloc(Nparams, PS_TYPE_F32);
+    tmp->params->n = tmp->params->nalloc;
+    tmp->dparams->n = tmp->dparams->nalloc;
 
     for (psS32 i = 0; i < tmp->params->n; i++) {
@@ -348,4 +351,5 @@
         if (vector->data.F32[0] > threshold) {
             tmpVector = psVectorAlloc(1, PS_TYPE_U32);
+            tmpVector->n = 1;
             tmpVector->data.U32[0] = 0;
         } else {
@@ -388,4 +392,5 @@
     //
     psVector *tmpVector = psVectorAlloc(count, PS_TYPE_U32);
+    tmpVector->n = tmpVector->nalloc;
     count = 0;
 
@@ -1428,4 +1433,5 @@
     //
     psVector *x = psVectorAlloc(2, PS_TYPE_F32);
+    x->n = x->nalloc;
     x->data.F32[0] = (psF32) (col + image->col0);
     x->data.F32[1] = (psF32) (row + image->row0);
@@ -1472,5 +1478,6 @@
     psVector *xVec = psVectorAlloc(2 * source->pixels->numRows, PS_TYPE_F32);
     psVector *yVec = psVectorAlloc(2 * source->pixels->numRows, PS_TYPE_F32);
-
+    xVec->n = xVec->nalloc;
+    yVec->n = yVec->nalloc;
     //
     // Start at the row with peak pixel, then decrement.
@@ -1552,4 +1559,5 @@
     //
     psArray *tmpArray = psArrayAlloc(2);
+    tmpArray->n = 2;
     tmpArray->data[0] = (psPtr *) yVec;
     tmpArray->data[1] = (psPtr *) xVec;
@@ -1616,4 +1624,7 @@
     psVector *y = psVectorAlloc(count, PS_TYPE_F32);
     psVector *yErr = psVectorAlloc(count, PS_TYPE_F32);
+    x->n = x->nalloc;
+    y->n = y->nalloc;
+    yErr->n = yErr->nalloc;
     psS32 tmpCnt = 0;
     for (psS32 i = 0; i < source->pixels->numRows; i++) {
@@ -1621,4 +1632,5 @@
             if (source->mask->data.U8[i][j] == 0) {
                 psVector *coord = psVectorAlloc(2, PS_TYPE_F32);
+                coord->n = 2;
                 // XXX: Convert i/j to image space:
                 // XXX EAM: coord order is (x,y) == (col,row)
@@ -1641,4 +1653,5 @@
     if (PSF) {
         paramMask = psVectorAlloc (params->n, PS_TYPE_U8);
+        paramMask->n = paramMask->nalloc;
         for (int i = 0; i < 4; i++) {
             paramMask->data.U8[i] = 0;
@@ -1673,4 +1686,5 @@
     if (paramMask != NULL) {
         psVector *delta = psVectorAlloc (params->n, PS_TYPE_F64);
+        delta->n = delta->nalloc;
         psMinimizeGaussNewtonDelta (delta, params, NULL, x, y, yErr, modelFunc);
         for (int i = 0; i < dparams->n; i++) {
@@ -1763,4 +1777,7 @@
     psVector *y = psVectorAlloc(count, PS_TYPE_F32);
     psVector *yErr = psVectorAlloc(count, PS_TYPE_F32);
+    x->n = x->nalloc;
+    y->n = y->nalloc;
+    yErr->n = yErr->nalloc;
     psS32 tmpCnt = 0;
     for (psS32 i = 0; i < source->pixels->numRows; i++) {
@@ -1768,4 +1785,5 @@
             if (source->mask->data.U8[i][j] == 0) {
                 psVector *coord = psVectorAlloc(2, PS_TYPE_F32);
+                coord->n = 2;
                 // XXX: Convert i/j to image space:
                 // XXX EAM: coord order is (x,y) == (col,row)
@@ -1779,5 +1797,6 @@
 
                 // this test enhances the weight based on deviation from the model flux
-                Ro = 1.0 + fabs (y->data.F32[tmpCnt] - ymodel) / sqrt(PS_SQR(ymodel - So) + PS_SQR(So));
+                Ro = 1.0 + fabs (y->data.F32[tmpCnt] - ymodel) / sqrt(PS_SQR(ymodel - So)
+                        + PS_SQR(So));
 
                 // psMinimizeLMChi2_EAM takes wt = 1/dY^2
@@ -1798,4 +1817,5 @@
     if (PSF) {
         paramMask = psVectorAlloc (params->n, PS_TYPE_U8);
+        paramMask->n = paramMask->nalloc;
         for (int i = 0; i < 4; i++) {
             paramMask->data.U8[i] = 0;
@@ -1846,4 +1866,5 @@
     if (paramMask != NULL) {
         psVector *delta = psVectorAlloc (params->n, PS_TYPE_F64);
+        delta->n = delta->nalloc;
         psMinimizeGaussNewtonDelta(delta, params, NULL, x, y, yErr, modelFunc);
         for (int i = 0; i < dparams->n; i++) {
@@ -1897,4 +1918,5 @@
 
     psVector *x = psVectorAlloc(2, PS_TYPE_F32);
+    x->n = 2;
     psVector *params = model->params;
     pmModelFunc modelFunc = pmModelFunc_GetFunction (model->type);
Index: trunk/psModules/src/objects/pmPSF.c
===================================================================
--- trunk/psModules/src/objects/pmPSF.c	(revision 6329)
+++ trunk/psModules/src/objects/pmPSF.c	(revision 6511)
@@ -6,6 +6,6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-12-24 01:24:32 $
+ *  @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
@@ -91,4 +91,5 @@
 
     psf->params = psArrayAlloc(Nparams - 4);
+    psf->params->n = psf->params->nalloc;
     for (int i = 0; i < psf->params->n; i++) {
         // XXX EAM : make this a user-defined value?
@@ -119,4 +120,8 @@
     psVector *z  = psVectorAlloc (models->n, PS_TYPE_F64);
     psVector *dz = psVectorAlloc (models->n, PS_TYPE_F64);
+    x->n = x->nalloc;
+    y->n = y->nalloc;
+    z->n = z->nalloc;
+    dz->n = dz->nalloc;
 
     for (int i = 0; i < models->n; i++) {
Index: trunk/psModules/src/objects/pmPSFtry.c
===================================================================
--- trunk/psModules/src/objects/pmPSFtry.c	(revision 6329)
+++ trunk/psModules/src/objects/pmPSFtry.c	(revision 6511)
@@ -5,6 +5,6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-12-24 01:24:32 $
+ *  @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
@@ -59,4 +59,10 @@
     test->mask     = psVectorAlloc (sources->n, PS_TYPE_U8);
 
+    test->modelFLT->n = test->modelFLT->nalloc;
+    test->modelPSF->n = test->modelPSF->nalloc;
+    test->metric->n   = test->metric->nalloc;
+    test->fitMag->n   = test->fitMag->nalloc;
+    test->mask->n     = test->mask->nalloc;
+
     for (int i = 0; i < test->modelFLT->n; i++) {
         test->mask->data.U8[i]  = 0;
@@ -204,4 +210,5 @@
     // rflux = ten(0.4*fitMag);
     psVector *rflux = psVectorAlloc (psfTry->sources->n, PS_TYPE_F64);
+    rflux->n = rflux->nalloc;
     for (int i = 0; i < psfTry->sources->n; i++) {
         if (psfTry->mask->data.U8[i] & PSFTRY_MASK_ALL)
@@ -219,4 +226,6 @@
     psVector *daBin = psVectorAlloc (rfBin->n, PS_TYPE_F64);
     psVector *maskB = psVectorAlloc (rfBin->n, PS_TYPE_U8);
+    daBin->n = daBin->nalloc;
+    maskB->n = maskB->nalloc;
     psFree (stats);
 
@@ -272,4 +281,5 @@
 
         psVector *tmp2 = psVectorAlloc (nKeep, PS_TYPE_F64);
+        tmp2->n = tmp2->nalloc;
         for (int j = 0; j < tmp2->n; j++) {
             tmp2->data.F64[j] = tmp->data.F64[j + nSkip];
@@ -337,4 +347,5 @@
     psVector *rflux = psVectorAlloc (try
                                      ->sources->n, PS_TYPE_F64);
+    rflux->n = rflux->nalloc;
     for (int i = 0; i < try
                 ->sources->n; i++) {