Very last change: log1m(1-w) to log(w)

Franzi2114 · Franzi2114 · commit e190fe032ef8 · 2025-12-08T09:16:36.000+01:00
diff --git a/stan/math/prim/prob/wiener4_lccdf.hpp b/stan/math/prim/prob/wiener4_lccdf.hpp
@@ -26,7 +26,7 @@ inline auto log_wiener_prob_hit_upper(const T_a& a, const T_v& v,
   const auto neg_v = -v;
   const auto one_m_w = 1.0 - w;
   if (fabs(v) == 0.0) {
-    return ret_t(log1m(one_m_w));
+    return ret_t(log(w));
   }
   const auto exponent = 2.0 * v * a * w;
   // This branch is for numeric stability

Original file line number	Diff line number	Diff line change
`@@ -26,7 +26,7 @@ inline auto log_wiener_prob_hit_upper(const T_a& a, const T_v& v,`
`26`	`26`	`const auto neg_v = -v;`
`27`	`27`	`const auto one_m_w = 1.0 - w;`
`28`	`28`	`if (fabs(v) == 0.0) {`
`29`		`- return ret_t(log1m(one_m_w));`
	`29`	`+ return ret_t(log(w));`
`30`	`30`	`}`
`31`	`31`	`const auto exponent = 2.0 * v * a * w;`
`32`	`32`	`// This branch is for numeric stability`