diff --git a/apps/probability/law/exponential_law.cpp b/apps/probability/law/exponential_law.cpp index f6a93b3a3..dd197ff23 100644 --- a/apps/probability/law/exponential_law.cpp +++ b/apps/probability/law/exponential_law.cpp @@ -69,7 +69,7 @@ float ExponentialLaw::cumulativeDistributiveFunctionAtAbscissa(float x) const { float ExponentialLaw::cumulativeDistributiveInverseForProbability(float * probability) { if (*probability >= 1.0f) { - return FLT_MAX+1.0f; + return INFINITY; } if (*probability <= 0.0f) { return 0.0f; diff --git a/apps/probability/law/law.cpp b/apps/probability/law/law.cpp index acdc90041..ae75b25fe 100644 --- a/apps/probability/law/law.cpp +++ b/apps/probability/law/law.cpp @@ -58,7 +58,7 @@ float Law::finiteIntegralBetweenAbscissas(float a, float b) const { float Law::cumulativeDistributiveInverseForProbability(float * probability) { if (*probability >= 1.0f) { - return FLT_MAX+1.0f; + return INFINITY; } if (!isContinuous()) { if (*probability <= 0.0f) { diff --git a/apps/probability/law/normal_law.cpp b/apps/probability/law/normal_law.cpp index 676439561..20eba841e 100644 --- a/apps/probability/law/normal_law.cpp +++ b/apps/probability/law/normal_law.cpp @@ -106,10 +106,10 @@ float NormalLaw::standardNormalCumulativeDistributiveFunctionAtAbscissa(float ab float NormalLaw::standardNormalCumulativeDistributiveInverseForProbability(float probability) { if (probability >= 1.0f) { - return FLT_MAX+1.0f; + return INFINITY; } if (probability <= 0.0f) { - return -FLT_MAX-1.0f; + return -INFINITY; } if (probability < 0.5f) { return -standardNormalCumulativeDistributiveInverseForProbability(1-probability); diff --git a/liba/include/math.h b/liba/include/math.h index 5495a952e..293a8b5b4 100644 --- a/liba/include/math.h +++ b/liba/include/math.h @@ -6,6 +6,8 @@ LIBA_BEGIN_DECLS #define NAN (0.0f/0.0f) +#define HUGE_VALF 1e50f +#define INFINITY HUGE_VALF #define M_E 2.71828182845904523536028747135266250 #define M_PI 3.14159265358979323846264338327950288