Index: trunk/ppStack/src/ppStackConvolve.c
===================================================================
--- trunk/ppStack/src/ppStackConvolve.c	(revision 35400)
+++ trunk/ppStack/src/ppStackConvolve.c	(revision 35455)
@@ -403,5 +403,6 @@
                     continue;
                 }
-                if ((options->matchChi2->data.F32[i] > thresh) || ! isfinite(options->matchChi2->data.F32[i])) {
+                if ((options->matchChi2->data.F32[i] > thresh) ||
+		    ! isfinite(options->matchChi2->data.F32[i])) {
                     numRej++;
                     options->inputMask->data.PS_TYPE_VECTOR_MASK_DATA[i] |= PPSTACK_MASK_CHI2;
Index: trunk/ppStack/src/ppStackMatch.c
===================================================================
--- trunk/ppStack/src/ppStackMatch.c	(revision 35400)
+++ trunk/ppStack/src/ppStackMatch.c	(revision 35455)
@@ -467,7 +467,10 @@
 		sum += kernels->rms;
                 num++;
+		//		psLogMsg("ppStack", PS_LOG_INFO, "KERNMATCHDUMP %d %d %g %g %g\n",
+		//			 index,num,kernels->mean,kernels->rms,sum);
             }
             psFree(iter);
             options->matchChi2->data.F32[index] = sum / (psImageCovarianceFactor(readout->covariance) * num);
+
         }
 
Index: trunk/ppStack/src/ppStackPrepare.c
===================================================================
--- trunk/ppStack/src/ppStackPrepare.c	(revision 35400)
+++ trunk/ppStack/src/ppStackPrepare.c	(revision 35455)
@@ -334,11 +334,5 @@
     //    fprintf(stderr,"means: %g %g\n",means->data.F32[0],means->data.F32[1]);
     //    fprintf(stderr,"sigma: %g %g \n",S->data.F32[0],S->data.F32[1]);
-
     //    fprintf(stderr,"pi:    %g %g\n",pi->data.F32[0],pi->data.F32[1]);
-    
-    //    unwrittenGMMfunction(options->inputSeeing, options->inputMask, 0xff,
-    //		 &Punimodal,
-    //		 &m1,&s1,&pi1,
-    //		 &m2,&s2,&pi2);
 
     // Use Gaussian mixture model analysis of the FWHM distribution to decide an optional FWHM limit
@@ -358,5 +352,9 @@
     if (limit > maxFWHM)    { limit = maxFWHM; }    // We should not be larger than our max
     if (limit < threshFWHM) { limit = threshFWHM; } // Nor smaller than our min
-
+    psLogMsg("ppStack",PS_LOG_INFO,
+	     "PSF FWHM distribution: limit: %f (%f %f %f) (%f %f %f) %f",
+	     limit,means->data.F32[0],S->data.F32[0],pi->data.F32[0],
+	     means->data.F32[1],S->data.F32[1],pi->data.F32[1],
+	     Punimodal);
     for (int i = 0; i < num; i++) {
       if (options->inputSeeing->data.F32[i] > limit) {
@@ -404,5 +402,5 @@
     if (options->convolve) {
         options->psf = ppStackPSF(config, numCols, numRows, psfs, options);
-        psFree(psfs);
+	psFree(psfs);
         if (!options->psf) {
 #if 1