Index: /trunk/ppStack/src/ppStackPrepare.c
===================================================================
--- /trunk/ppStack/src/ppStackPrepare.c	(revision 35399)
+++ /trunk/ppStack/src/ppStackPrepare.c	(revision 35400)
@@ -1,3 +1,9 @@
 # include "ppStack.h"
+
+
+// Option to use the old FWHM rejection code.
+# define NEW_REJECTION
+// #define OLD_REJECTION
+
 
 # define RE_PHOTOMETER 0
@@ -293,4 +299,5 @@
             psLogMsg("ppStack", PS_LOG_INFO, "PSF FWHM for image %d is too large (%f vs %f maxFWHM) --- rejected.", i, options->inputSeeing->data.F32[i], maxFWHM);
 	}
+	// End old rejection fixed limit.
 #endif
         psStringAppend(&log, "Input %d: %f\n", i, options->inputSeeing->data.F32[i]);
@@ -301,4 +308,5 @@
     psFree(log);
 
+#ifdef NEW_REJECTION
     // Do GMM test with two modes to decide where to put the break.
     double Punimodal;
@@ -324,8 +332,8 @@
       return(false);
     }
-    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]);
+    //    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,
@@ -364,4 +372,6 @@
     psFree(pi);
     psFree(P);
+    // End new rejection code
+#endif
     
 #ifdef OLD_REJECTION
@@ -389,4 +399,5 @@
         }
     }
+    // End old rejection sigma clip
 #endif
     // Generate target PSF
Index: /trunk/psLib/src/math/psMixtureModels.c
===================================================================
--- /trunk/psLib/src/math/psMixtureModels.c	(revision 35399)
+++ /trunk/psLib/src/math/psMixtureModels.c	(revision 35400)
@@ -42,5 +42,5 @@
 
 
-#define TESTING
+// #define TESTING
 /////////////////////////////////////////////////////////////////////////////////
 // Internal functions                                                          //