Index: /trunk/psModules/src/detrend/pmSubtractBias.c
===================================================================
--- /trunk/psModules/src/detrend/pmSubtractBias.c	(revision 7444)
+++ /trunk/psModules/src/detrend/pmSubtractBias.c	(revision 7445)
@@ -11,6 +11,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-06-02 00:55:23 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-06-08 22:14:28 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -213,4 +213,6 @@
     psVector *ordinate = psVectorAlloc(pixels->n, PS_TYPE_F32); // Ordinate
     psVector *mask = psVectorAlloc(pixels->n, PS_TYPE_U8); // Mask for fitting
+    reduced->n = ordinate->n = mask->n = pixels->n;
+
     for (int i = 0; i < pixels->n; i++) {
         psVector *values = pixels->data[i]; // Vector with overscan values
@@ -328,5 +330,4 @@
         if (overscanOpts->single) {
             psVector *pixels = psVectorAlloc(0, PS_TYPE_F32);
-            pixels->n = 0;
             psListIterator *iter = psListIteratorAlloc(overscans, PS_LIST_HEAD, false); // Iterator
             psImage *overscan = NULL;   // Overscan image from iterator
@@ -360,7 +361,7 @@
                 // The read direction is rows
                 psArray *pixels = psArrayAlloc(image->numRows); // Array of vectors containing pixels
+                pixels->n = image->numRows;
                 for (int i = 0; i < pixels->n; i++) {
                     psVector *values = psVectorAlloc(0, PS_TYPE_F32);
-                    values->n = 0;
                     pixels->data[i] = values;
                 }
@@ -404,7 +405,7 @@
                 // The read direction is columns
                 psArray *pixels = psArrayAlloc(image->numCols); // Array of vectors containing pixels
+                pixels->n = image->numCols;
                 for (int i = 0; i < pixels->n; i++) {
                     psVector *values = psVectorAlloc(0, PS_TYPE_F32);
-                    values->n = 0;
                     pixels->data[i] = values;
                 }