Merge pull request #27 from brgirgis/fix-fdist

Fix F dist at X=Inf
kthohr · Mar 13, 2021 · fc1f7d8 · fc1f7d8
2 parents 590bfc1 + 794f802
commit fc1f7d8
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/include/stats_incl/dens/df.ipp b/include/stats_incl/dens/df.ipp
@@ -100,6 +100,8 @@ noexcept
                 log_zero_if<T>(log_form) :
             x == T(0) ? \
                 log_if(df_limit_vals_x(df1_par), log_form) :
+            GCINT::is_posinf(x) ? \
+                T(0) :
             //
             GCINT::any_posinf(df1_par,df2_par) ? \
                 df_limit_vals_dof(x,df1_par,df2_par,log_form) :

diff --git a/include/stats_incl/prob/pf.ipp b/include/stats_incl/prob/pf.ipp
@@ -66,6 +66,8 @@ pf_vals_check(const T x, const T df1_par, const T df2_par, const bool log_form)
             //
             STLIM<T>::epsilon() > x ? \
                 log_zero_if<T>(log_form) :
+            GCINT::is_posinf(x) ? \
+                T(1) :
             //
             GCINT::any_posinf(df1_par,df2_par) ? \
                 pf_limit_vals_dof(x,df1_par,df2_par,log_form) :