Compare to absolute eps instead of relative eps

Jutho · Jun 14, 2024 · f3c06a8 · f3c06a8
1 parent ecbd858
commit f3c06a8
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/OptimKit.jl b/src/OptimKit.jl
@@ -17,8 +17,9 @@ function safe_sqrt(x::Real)
     return if x >= 0
         sqrt(x)
     else
-        -x < eps(x)^(3 / 4) || @error "sqrt of negative number: $x"
-        eps(x)
+        ϵ = eps(typeof(x))
+        -x < ϵ^(3 / 4) || @error "sqrt of negative number: $x"
+        ϵ
     end
 end