One more edge case-test+fix

stan-dev · Mar 12, 2020 · 6aed316 · 6aed316
1 parent c84a259
commit 6aed316
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/stan/math/prim/fun/binomial_coefficient_log.hpp b/stan/math/prim/fun/binomial_coefficient_log.hpp
@@ -84,7 +84,7 @@ inline return_type_t<T_n, T_k> binomial_coefficient_log(const T_n n,
   }
 
   // Choosing the more stable of the symmetric branches
-  if (n > 0 && k > value_of_rec(n) / 2.0 + 1e-8) {
+  if (n > -1 && k > value_of_rec(n) / 2.0 + 1e-8) {
     return binomial_coefficient_log(n, n - k);
   }
 

diff --git a/test/unit/math/rev/fun/binomial_coefficient_log_test.cpp b/test/unit/math/rev/fun/binomial_coefficient_log_test.cpp
@@ -106,6 +106,8 @@ std::vector<TestValue> testValuesEdge = {
     {3, -1, stan::math::NEGATIVE_INFTY, -0.25, stan::math::INFTY},
     {-1, -0.2, stan::math::INFTY, stan::math::NEGATIVE_INFTY,
      -4.324031329886049836},
+    {-0.5, 0.5, stan::math::NEGATIVE_INFTY, stan::math::INFTY,
+     stan::math::NEGATIVE_INFTY},
     {4.0, 5.0, stan::math::NEGATIVE_INFTY, stan::math::INFTY,
      stan::math::NEGATIVE_INFTY},
     {1, 0, 0, 0, 1},