Changeset 26031

Timestamp:

Nov 4, 2009, 2:29:14 PM (17 years ago)

Author:

Paul Price

Message:

Ensure threshold is positive.

File:

• trunk/psModules/src/imcombine/pmSubtractionEquation.c (modified) (1 diff)

trunk/psModules/src/imcombine/pmSubtractionEquation.c

@@ -1016 +1016 @@
             #define TOL 1.0e-5
             int normIndex = PM_SUBTRACTION_INDEX_NORM(kernels); // Index for normalisation
-            double norm = solution->data.F64[normIndex];        // Normalisation
+            double norm = fabs(solution->data.F64[normIndex]);  // Normalisation
             double thresh = norm * TOL;                         // Threshold for low parameters
             for (int i = 0; i < numKernels; i++) {