Index: /branches/eam_branches/ipp-20211108/psModules/src/imcombine/pmStack.c
===================================================================
--- /branches/eam_branches/ipp-20211108/psModules/src/imcombine/pmStack.c	(revision 41948)
+++ /branches/eam_branches/ipp-20211108/psModules/src/imcombine/pmStack.c	(revision 41949)
@@ -1657,4 +1657,8 @@
 	    int Nhi = nGood; 
 
+	    if ((x == 3161) && (y == 450)) {
+	      fprintf (stderr, "test point\n");
+	    }
+
 	    if (nGood >= 5) {
 	      int midPoint = nGood / 2;
@@ -1665,6 +1669,6 @@
 	      int P75 = 0.75*nGood;
 	      float rawSigma = 0.74*(pixelData->data.F32[P75] - pixelData->data.F32[P25]);
-	      float minThresh = rawMedian - 7.0*rawSigma;
-	      float maxThresh = rawMedian + 7.0*rawSigma;
+	      float minThresh = rawMedian - 5.0*rawSigma;
+	      float maxThresh = rawMedian + 5.0*rawSigma;
 
 	      // find the entries which are in the range
@@ -1691,9 +1695,9 @@
 # if (AT_LEAST)
 	    int Ns = MIN(MAX(1, 0.5*rejectFraction * nGoodClip), nGoodClip) + Nlo;
-	    int Ne = nGood - Ns + Nlo;
+	    int Ne = nGoodClip - Ns + Nlo;
 	    int Npt = Ne - Ns;
 # else
-	    int Ns = MIN(MAX(0, 0.5*rejectFraction * nGood), nGood) + Nlo;
-	    int Ne = nGood - Ns + Nlo;
+	    int Ns = MIN(MAX(0, 0.5*rejectFraction * nGoodClip), nGoodClip) + Nlo;
+	    int Ne = nGoodClip - Ns + Nlo;
 	    int Npt = Ne - Ns;
 # endif