add better error checking for dual numbers

JuliaMath · Jul 9, 2020 · 720e95d · 720e95d · MikaelSlevinsky · Jul 9, 2020
1 parent 123e404
commit 720e95d
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "HypergeometricFunctions"
 uuid = "34004b35-14d8-5ef3-9330-4cdb6864b03a"
-version = "0.3"
+version = "0.3.1"
 
 [deps]
 DualNumbers = "fa6b7ba4-c1ee-5f82-b5fc-ecf0adba8f74"

diff --git a/src/specialfunctions.jl b/src/specialfunctions.jl
@@ -1,4 +1,7 @@
-@inline errcheck(x, y, tol) = isfinite(x) && isfinite(y) && (norm(x-y) > max(norm(x), norm(y))*tol)
+@inline norm2(x) = norm(x)
+@inline norm2(x::Dual) = hypot(x.value, x.epsilon)
+
+@inline errcheck(x, y, tol) = isfinite(x) && isfinite(y) && (norm2(x-y) > max(norm2(x), norm2(y))*tol)
 
 # Same as in FastTransforms.jl
 """