Changeset 41520 for trunk/psLib

Timestamp:

Apr 2, 2021, 2:12:52 PM (5 years ago)

Author:

eugene

Message:

fix index for UQ and LQ (was ignored, but still)

File:

• trunk/psLib/src/math/psStats.c (modified) (3 diffs)

trunk/psLib/src/math/psStats.c

@@ -377 +377 @@
     }
 
-    int Qmin = (int)0.25*count;
-    int Qmax = (int)0.75*count;
+    int Qmin = (int)(0.25*count);
+    int Qmax = (int)(0.75*count);
     psAssert (Qmin >=           0, "oops");
     psAssert (Qmin < outVector->n, "oops");
@@ -385 +385 @@
 
     // Calculate the quartile points exactly.
-    stats->sampleUQ = output[(int)(0.75*count)];
-    stats->sampleLQ = output[(int)(0.25*count)];
-
+    stats->sampleUQ = output[Qmax];
+    stats->sampleLQ = output[Qmin];
+     
     stats->results |= PS_STAT_SAMPLE_MEDIAN;
     stats->results |= PS_STAT_SAMPLE_QUARTILE;
@@ -1915 +1915 @@
         // NOTE: vectorSampleMedian only returns 'false' for very bad cases:
         // NULL input vector, psVectorCopy failure, invalid vector type (likely programming errors)
-        if (!vectorSampleMedian(inF32, maskVector, maskVal, stats)) {
-            psError(PS_ERR_UNKNOWN, false, "Failed to calculate sample median");
-            status &= false;
-        }
+        if (!vectorSampleMedian(inF32, maskVector, maskVal, stats)) {
+            psError(PS_ERR_UNKNOWN, false, "Failed to calculate sample median");
+            status &= false;
+        }
     }