@@ -34,15 +34,15 @@ namespace math {
3434 * @param alpha prior success parameter
3535 * @param beta prior failure parameter
3636 * @param precision precision for `grad_F32`, default \f$10^{-8}\f$
37- * @param max_steps max iteration allowed for `grad_F32`, default \f$10^{- 8}\f$
37+ * @param max_steps max iteration allowed for `grad_F32`, default \f$10^{8}\f$
3838 * @return log probability or log sum of probabilities
3939 * @throw std::domain_error if r, alpha, or beta fails to be positive
4040 * @throw std::invalid_argument if container sizes mismatch
4141 */
4242template <typename T_n, typename T_r, typename T_alpha, typename T_beta>
4343inline return_type_t <T_r, T_alpha, T_beta> beta_neg_binomial_lccdf (
4444 const T_n& n, const T_r& r, const T_alpha& alpha, const T_beta& beta,
45- const double precision = 1e-8 , const int max_steps = 1e6 ) {
45+ const double precision = 1e-8 , const int max_steps = 1e8 ) {
4646 static constexpr const char * function = " beta_neg_binomial_lccdf" ;
4747 check_consistent_sizes (
4848 function, " Failures variable" , n, " Number of successes parameter" , r,
@@ -100,7 +100,7 @@ inline return_type_t<T_r, T_alpha, T_beta> beta_neg_binomial_lccdf(
100100 std::initializer_list<b_t >{n_dbl + 2.0 , a_plus_r + b_plus_n + 1.0 },
101101 1.0 );
102102 auto C = lgamma (r_plus_n + 1.0 ) + lbeta (a_plus_r, b_plus_n + 1.0 )
103- - lgamma (r_dbl) - lbeta (alpha_dbl, beta_dbl) - lgamma (n_dbl + 2 );
103+ - lgamma (r_dbl) - lbeta (alpha_dbl, beta_dbl) - lgamma (n_dbl + 2.0 );
104104 log_ccdf += C + stan::math::log (F);
105105
106106 if constexpr (!is_constant_all<T_r, T_alpha, T_beta>::value) {
0 commit comments