@@ -48,8 +48,7 @@ struct log_gamma_q_result {
4848 */
4949template <typename T_a, typename T_z>
5050inline log_gamma_q_result<return_type_t <T_a, T_z>> log_gamma_q_dgamma (
51- const T_a& a, const T_z& z, int max_steps = 250 ,
52- double precision = 1e-16 ) {
51+ const T_a& a, const T_z& z, int max_steps = 250 , double precision = 1e-16 ) {
5352 using std::exp;
5453 using std::fabs;
5554 using std::log;
@@ -99,7 +98,8 @@ inline log_gamma_q_result<return_type_t<T_a, T_z>> log_gamma_q_dgamma(
9998 // For gradient, use: d/da log(Q) = (1/Q) * dQ/da
10099 // grad_reg_inc_gamma computes dQ/da
101100 const double Q_val = exp (result.log_q );
102- const double dQ_da = grad_reg_inc_gamma (a_dbl, z_dbl, tgamma (a_dbl), digamma (a_dbl));
101+ const double dQ_da
102+ = grad_reg_inc_gamma (a_dbl, z_dbl, tgamma (a_dbl), digamma (a_dbl));
103103 result.dlog_q_da = dQ_da / Q_val;
104104
105105 } else {
@@ -111,7 +111,8 @@ inline log_gamma_q_result<return_type_t<T_a, T_z>> log_gamma_q_dgamma(
111111 // grad_reg_inc_gamma computes dQ/da
112112 const double Q_val = exp (result.log_q );
113113 if (Q_val > 0 ) {
114- const double dQ_da = grad_reg_inc_gamma (a_dbl, z_dbl, tgamma (a_dbl), digamma (a_dbl));
114+ const double dQ_da
115+ = grad_reg_inc_gamma (a_dbl, z_dbl, tgamma (a_dbl), digamma (a_dbl));
115116 result.dlog_q_da = dQ_da / Q_val;
116117 } else {
117118 // Fallback if Q rounds to zero - use asymptotic approximation
0 commit comments