@@ -156,9 +156,9 @@ inline auto wiener4_distribution(const T_y& y, const T_a& a, const T_v& v,
156156 const auto arg = fmax (
157157 0.0 , fmin (1.0 , exp (one_m_w_a_neg_v + square (neg_v) * y / 2.0 + log_err)
158158 / 2.0 ));
159- const auto K2 = (arg == 0 ) ? INFTY
160- : (arg == 1 ) ? NEGATIVE_INFTY
161- : -sqrt (y) / 2.0 / a * inv_Phi (arg);
159+ const auto K2 = (arg == 0 ) ? INFTY
160+ : (arg == 1 ) ? NEGATIVE_INFTY
161+ : -sqrt (y) / 2.0 / a * inv_Phi (arg);
162162 const auto K_small_value = ceil (fmax (K1, K1 + K2));
163163
164164 const auto api = a / pi ();
@@ -326,9 +326,9 @@ inline auto wiener4_cdf_grad_a(const T_y& y, const T_a& a, const T_v& v,
326326 = make_finite (exp (log_probability_distribution (a, neg_v, one_m_w)));
327327 const auto dav = log_probability_GradAV (a, neg_v, one_m_w);
328328 auto dav_neg_v = dav * neg_v;
329- auto prob_deriv = fabs (neg_v) == 0 ? ret_t ( 0.0 )
330- : is_inf (dav_neg_v) ? NEGATIVE_INFTY
331- : dav_neg_v * prob;
329+ auto prob_deriv = fabs (neg_v) == 0
330+ ? ret_t ( 0.0 )
331+ : is_inf (dav_neg_v) ? NEGATIVE_INFTY : dav_neg_v * prob;
332332 ans = (-2.0 / a - one_m_w_neg_v) * (cdf - prob)
333333 + ans * (2.0 * pi () / square (a))
334334 * exp (-one_m_w_a_neg_v - 0.5 * square (neg_v) * y);
@@ -371,9 +371,10 @@ inline auto wiener4_cdf_grad_v(const T_y& y, const T_a& a, const T_v& v,
371371 - 2.5 * LOG_TWO - 3 * log_a - log_v),
372372 1.0 );
373373 K_large_value
374- = fmax (ceil ((alphK_large == 0 ) ? ret_t (INFTY)
375- : (alphK_large == 1 ) ? ret_t (NEGATIVE_INFTY)
376- : temp * inv_Phi (alphK_large)),
374+ = fmax (ceil ((alphK_large == 0 )
375+ ? ret_t (INFTY)
376+ : (alphK_large == 1 ) ? ret_t (NEGATIVE_INFTY)
377+ : temp * inv_Phi (alphK_large)),
377378 ret_t (1.0 ));
378379 }
379380
@@ -468,9 +469,10 @@ inline auto wiener4_cdf_grad_w(const T_y& y, const T_a& a, const T_v& v,
468469 = fmin (exp (factor + 0.5 * (LOG_PI + log_y) - 1.5 * LOG_TWO - log_a), 1.0 );
469470 alphK_large = fmax (0.0 , alphK_large);
470471 const auto K_large_value
471- = fmax (ceil ((alphK_large == 0 ) ? ret_t (INFTY)
472- : (alphK_large == 1 ) ? ret_t (NEGATIVE_INFTY)
473- : temp * inv_Phi (alphK_large)),
472+ = fmax (ceil ((alphK_large == 0 )
473+ ? ret_t (INFTY)
474+ : (alphK_large == 1 ) ? ret_t (NEGATIVE_INFTY)
475+ : temp * inv_Phi (alphK_large)),
474476 ret_t (1.0 ));
475477
476478 const auto sqrt_y = sqrt (y);
@@ -479,9 +481,10 @@ inline auto wiener4_cdf_grad_w(const T_y& y, const T_a& a, const T_v& v,
479481 const auto lv = log1p (square (neg_v) * y);
480482 const auto alphK_small = factor - LOG_TWO - lv;
481483 const auto arg = fmin (exp (alphK_small), 1.0 );
482- const auto K_small = (arg == 0 ) ? INFTY
483- : (arg == 1 ) ? NEGATIVE_INFTY
484- : -sqrt_y / a * inv_Phi (arg) - wdash;
484+ const auto K_small
485+ = (arg == 0 )
486+ ? INFTY
487+ : (arg == 1 ) ? NEGATIVE_INFTY : -sqrt_y / a * inv_Phi (arg) - wdash;
485488 const auto K_small_value = ceil (fmax (fmax (K_small, K_large), ret_t (1.0 )));
486489
487490 if (K_large_value > 4 * K_small_value) {
0 commit comments