[apps/proba] Use DBL_EPSILON in cnditions

This commit is contained in:
Léa Saviot
2019-09-06 10:23:03 +02:00
parent d9e49cb369
commit 80f3c1fd0a

View File

@@ -105,10 +105,10 @@ double Distribution::evaluateAtDiscreteAbscissa(int k) const {
double Distribution::cumulativeDistributiveInverseForProbabilityUsingIncreasingFunctionRoot(double * probability, double ax, double bx) {
assert(ax < bx);
if (*probability >= 1.0) {
if (*probability > 1.0 - DBL_EPSILON) {
return INFINITY;
}
if (*probability <= 0.0) {
if (*probability < DBL_EPSILON) {
return -INFINITY;
}
Poincare::Coordinate2D<double> result = Poincare::Solver::IncreasingFunctionRoot(