use log1p

Expander · Oct 4, 2023 · e1bb472 · e1bb472
1 parent afe6979
commit e1bb472
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Li2.jl b/src/Li2.jl
@@ -47,7 +47,7 @@ end
 # approximation of Re[Li2(x)] for x in [0, 1/2]
 # todo(alex): benchmark this routine
 function reli2_approx(x::BigFloat)::BigFloat
-    u = -log(one(x) - x)
+    u = -log1p(-x)
     u2 = u*u
 
     p = u2 # powers of u2