Add rtol/atol options for computing roots to avoid numerical errors

serenity4 · Aug 17, 2024 · 3cf433e · 3cf433e
1 parent 0c42cad
commit 3cf433e
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/frontend/MathFunctions.jl b/src/frontend/MathFunctions.jl
@@ -36,8 +36,8 @@ function rotate_2d(v, angle)
   typeof(v)(real(cv′), imag(cv′))
 end
 
-function compute_roots(a, b, c)
-  if isapprox(a, zero(a), atol = 1e-7)
+function compute_roots(a, b, c; rtol = 1e-6F, atol = 0F)
+  if isapprox(a, zero(a); rtol, atol)
     t₁ = t₂ = c / 2b
     return (t₁, t₂)
   end