fix AD tests

Author: heltonmc

src/BesselFunctions/besselk.jl

@@ -244,7 +244,7 @@ function _besselk(v::T, x::T) where T <: Union{Float32, Float64}
 
     else
         v_floor, v_int = modf(v)
-        if x >= 1.5 # determine cutoff as function for differnet types
+        if x > 1.5 # determine cutoff as function for differnet types
             kv, kvp1 = besselkx_levin(v_floor, x, Val(16)), besselkx_levin(v_floor + 1, x, Val(16))
             return besselk_up_recurrence(x, kvp1, kv, v_floor + 1, v)[1] * exp(-x)
         else

test/besselk_enzyme_test.jl

@@ -8,10 +8,10 @@ dbesselkx_dv(v, x) = autodiff(Forward, _v->besselkx_levin(_v, x, Val(30)),
 dbesselkx_dx(v, x) = autodiff(Forward, _x->besselkx_levin(v, _x, Val(30)), 
                               Duplicated, Duplicated(x, 1.0))[2]
 
-dbesselk_ps_dv(v, x) = autodiff(Forward, _v->besselk_power_series(_v, x), 
+dbesselk_ps_dv(v, x) = autodiff(Forward, _v->besselk(_v, x), 
                                 Duplicated, Duplicated(v, 1.0))[2]
 
-dbesselk_ps_dx(v, x) = autodiff(Forward, _x->besselk_power_series(v, _x), 
+dbesselk_ps_dx(v, x) = autodiff(Forward, _x->besselk(v, _x), 
                                 Duplicated, Duplicated(x, 1.0))[2]