add better error checking for dual numbers

MikaelSlevinsky · MikaelSlevinsky · commit 720e95dd2e99 · 2020-07-09T10:40:54.000-05:00
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
 """
