Update cbrt to use ldexpk

musm · Jun 3, 2017 · 50badb5 · 50badb5
1 parent e1c240b
commit 50badb5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/misc.jl b/src/misc.jl
@@ -96,7 +96,7 @@ function cbrt(d::T) where {T<:IEEEFloat}
     v  = dadd(dsqu(z), y)
     v  = dmul(v, d)
     v  = dmul(v, q2)
-    z  = ldexp(T(v), (e + 6144) ÷ 3 - 2048)
+    z  = ldexpk(T(v), (e + 6144) ÷ 3 - 2048)
     isinf(d) && (z = flipsign(T(Inf), q2.hi))
     d == 0   && (z = flipsign(T(0), q2.hi))
     return z