Fix detection of non strictly increasing CDF

josephmure · Apr 12, 2024 · 93f1404 · 93f1404
1 parent caa14a0
commit 93f1404
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/lib/src/Uncertainty/Model/DistributionImplementation.cxx b/lib/src/Uncertainty/Model/DistributionImplementation.cxx
@@ -2371,7 +2371,7 @@ Scalar DistributionImplementation::computeScalarQuantile(const Scalar prob,
   LOGDEBUG(OSS() << "root=" << root);
 
   // special case non strictly increasing CDF: retain the inf of the interval veryfing F(x)=p
-  if (computeCDF(root - quantileEpsilon_) == prob)
+  if (root > quantileEpsilon_ && computeCDF(root - quantileEpsilon_) == prob)
   {
     solver.setResidualError(0.5 * cdfEpsilon_);
     root = solver.solve(f, p - cdfEpsilon_, leftTau, root, leftCDF, prob);