Index: trunk/psLib/src/math/psMinimizePowell.c
===================================================================
--- trunk/psLib/src/math/psMinimizePowell.c	(revision 6449)
+++ trunk/psLib/src/math/psMinimizePowell.c	(revision 6484)
@@ -11,6 +11,6 @@
  *  NOTE: XXX: The SDR is silent about data types.  F32 is implemented here.
  *
- *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-02-18 00:56:44 $
+ *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-02-24 23:43:15 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -525,6 +525,7 @@
         p_psMemSetPersistent(myParamMask, true);
         p_psMemSetPersistent(myParamMask->data.U8, true);
-        for (i=0;i<myParamMask->n;i++) {
+        for (i=0;i<myParamMask->nalloc;i++) {
             myParamMask->data.U8[i] = 0;
+            myParamMask->n++;
         }
     } else {
@@ -543,4 +544,5 @@
                 ((psVector *) (v->data[i]))->data.F32[j] = 0.0;
             }
+            ((psVector *)(v->data[i]))->n++;
         }
         v->n++;
@@ -550,4 +552,5 @@
     for (i=0;i<numDims;i++) {
         Q->data.F32[i] = params->data.F32[i];
+        Q->n++;
     }
 
@@ -595,4 +598,5 @@
             if (myParamMask->data.U8[i] == 0) {
                 u->data.F32[i] = Q->data.F32[i] - params->data.F32[i];
+                u->n++;
 
                 psTrace(__func__, 6, "u[i]=Q[i]-P[i] (%f = %f - %f)\n", u->data.F32[i],
@@ -602,4 +606,5 @@
             } else {
                 u->data.F32[i] = 0.0;
+                u->n++;
             }
         }