comparison DEPENDENCIES/generic/include/boost/math/distributions/negative_binomial.hpp @ 101:c530137014c0

Update Boost headers (1.58.0)

       }
       if (P <= pow(dist.success_fraction(), dist.successes()))
       { // p <= pdf(dist, 0) == cdf(dist, 0)
         return 0;
       }
+      if(p == 0)
+      {  // Would need +infinity failures for total confidence.
+         result = policies::raise_overflow_error<RealType>(
+            function,
+            "Success fraction is 0, which implies infinite failures !", Policy());
+         return result;
+         // usually means return +std::numeric_limits<RealType>::infinity();
+         // unless #define BOOST_MATH_THROW_ON_OVERFLOW_ERROR
+      }
       /*
       // Calculate quantile of negative_binomial using the inverse incomplete beta function.
       using boost::math::ibeta_invb;
       return ibeta_invb(r, p, P, Policy()) - 1; //
       */

        if(Q == 1)
        {  // There may actually be no answer to this question,
           // since the probability of zero failures may be non-zero,
           return 0; // but zero is the best we can do:
        }
-       if (-Q <= boost::math::powm1(dist.success_fraction(), dist.successes(), Policy()))
-       {  // q <= cdf(complement(dist, 0)) == pdf(dist, 0)
-          return 0; //
-       }
        if(Q == 0)
        {  // Probability 1 - Q  == 1 so infinite failures to achieve certainty.
           // Would need +infinity failures for total confidence.
           result = policies::raise_overflow_error<RealType>(
              function,
              "Probability argument complement is 0, which implies infinite failures !", Policy());
+          return result;
+          // usually means return +std::numeric_limits<RealType>::infinity();
+          // unless #define BOOST_MATH_THROW_ON_OVERFLOW_ERROR
+       }
+       if (-Q <= boost::math::powm1(dist.success_fraction(), dist.successes(), Policy()))
+       {  // q <= cdf(complement(dist, 0)) == pdf(dist, 0)
+          return 0; //
+       }
+       if(p == 0)
+       {  // Success fraction is 0 so infinite failures to achieve certainty.
+          // Would need +infinity failures for total confidence.
+          result = policies::raise_overflow_error<RealType>(
+             function,
+             "Success fraction is 0, which implies infinite failures !", Policy());
           return result;
           // usually means return +std::numeric_limits<RealType>::infinity();
           // unless #define BOOST_MATH_THROW_ON_OVERFLOW_ERROR
        }
        //return ibetac_invb(r, p, Q, Policy()) -1;