Index: /trunk/psModules/src/imcombine/pmSubtraction.c
===================================================================
--- /trunk/psModules/src/imcombine/pmSubtraction.c	(revision 14715)
+++ /trunk/psModules/src/imcombine/pmSubtraction.c	(revision 14716)
@@ -4,6 +4,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.53 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-08-30 22:06:04 $
+ *  @version $Revision: 1.54 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-08-30 23:09:52 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -641,6 +641,6 @@
                     for (int jxOrder = 0; jxOrder <= spatialOrder - jyOrder;
                          jxOrder++, jIndex += numKernels) {
-                        for (int kyOrder = jyOrder, kIndex = k; kyOrder <= spatialOrder; kyOrder++) {
-                            for (int kxOrder = jxOrder; kxOrder <= spatialOrder - kyOrder;
+                        for (int kyOrder = 0, kIndex = k; kyOrder <= spatialOrder; kyOrder++) {
+                            for (int kxOrder = 0; kxOrder <= spatialOrder - kyOrder;
                                  kxOrder++, kIndex += numKernels) {
                                 double convPoly = sumCC * polyValues->data.F64[jyOrder][jxOrder] *
@@ -713,4 +713,5 @@
             psFitsClose(matrixFile);
         }
+
     }