@@ -79,46 +79,45 @@ inline return_type_t<T_r, T_alpha, T_beta> beta_neg_binomial_lpmf(
7979 if constexpr (include_summand<propto>::value) {
8080 logp -= lgamma (n_vec[i] + 1 );
8181 }
82- T_partials_return lbeta_denominator
83- = lbeta (r_vec.val (i), alpha_vec.val (i));
84- T_partials_return lgamma_numerator
85- = lgamma (n_vec[i] + beta_vec.val (i));
82+ T_partials_return lbeta_denominator = lbeta (r_vec.val (i), alpha_vec.val (i));
83+ T_partials_return lgamma_numerator = lgamma (n_vec[i] + beta_vec.val (i));
8684 T_partials_return lgamma_denominator = lgamma (beta_vec.val (i));
8785 T_partials_return lbeta_numerator
8886 = lbeta (n_vec[i] + r_vec.val (i), alpha_vec.val (i) + beta_vec.val (i));
8987 logp += lbeta_numerator + lgamma_numerator - lbeta_denominator
9088 - lgamma_denominator;
9189 if (!is_constant_all<T_r, T_alpha, T_beta>::value) {
92- T_partials_return digamma_n_r_alpha_beta
93- = digamma (n_vec[i] + r_vec.val (i) + alpha_vec.val (i)
94- + beta_vec.val (i));
90+ T_partials_return digamma_n_r_alpha_beta = digamma (
91+ n_vec[i] + r_vec.val (i) + alpha_vec.val (i) + beta_vec.val (i));
9592
9693 if constexpr (!is_constant<T_r>::value || !is_constant<T_alpha>::value) {
97- T_partials_return digamma_r_alpha = digamma (r_vec.val (i) + alpha_vec.val (i));
94+ T_partials_return digamma_r_alpha
95+ = digamma (r_vec.val (i) + alpha_vec.val (i));
9896 if constexpr (!is_constant_all<T_r>::value) {
9997 partials<0 >(ops_partials)[i]
10098 += digamma (n_vec[i] + r_vec.val (i)) - digamma_n_r_alpha_beta
101- - (digamma (r_vec.val (i)) - digamma_r_alpha);
99+ - (digamma (r_vec.val (i)) - digamma_r_alpha);
102100 }
103101 if constexpr (!is_constant_all<T_alpha>::value) {
104102 partials<1 >(ops_partials)[i]
105103 += -digamma_n_r_alpha_beta
106- - (digamma (alpha_vec.val (i)) - digamma_r_alpha);
104+ - (digamma (alpha_vec.val (i)) - digamma_r_alpha);
107105 }
108106 }
109- if constexpr (!is_constant<T_beta>::value || !is_constant<T_alpha>::value) {
107+ if constexpr (!is_constant<T_beta>::value
108+ || !is_constant<T_alpha>::value) {
110109 T_partials_return digamma_alpha_beta
111110 = digamma (alpha_vec.val (i) + beta_vec.val (i));
112111 if constexpr (!is_constant_all<T_beta>::value) {
113- partials<2 >(ops_partials)[i]
114- += digamma_alpha_beta - digamma_n_r_alpha_beta
115- + digamma (n_vec[i] + beta_vec.val (i)) - digamma (beta_vec.val (i));
112+ partials<2 >(ops_partials)[i] += digamma_alpha_beta
113+ - digamma_n_r_alpha_beta
114+ + digamma (n_vec[i] + beta_vec.val (i))
115+ - digamma (beta_vec.val (i));
116116 }
117117 if constexpr (!is_constant_all<T_alpha>::value) {
118118 partials<1 >(ops_partials)[i] += digamma_alpha_beta;
119119 }
120120 }
121-
122121 }
123122 }
124123 return ops_partials.build (logp);
0 commit comments