[kepler_solver] Improve type-stability

JuliaAstro · Apr 1, 2021 · 5493395 · 5493395
1 parent 7fe0411
commit 5493395
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/kepler_solver.jl b/src/kepler_solver.jl
@@ -11,7 +11,7 @@ function kepler_solver(_M::Real, e::Real)
     end
     pi2 = abs2(T(pi))
     # equation (20)
-    α = (3 * pi2 + 1.6 * (pi2 - pi * abs(M))/(1 + e))/(pi2 - 6)
+    α = (3 * pi2 + 8 * (pi2 - pi * abs(M)) / (5 * (1 + e)))/(pi2 - 6)
     # equation (5)
     d = 3 * (1 - e) + α * e
     # equation (9)