Problem with airyaiprime
macd opened this issue · 3 comments
macd commented
Maybe I'm doing it wrong, but it looks like there is a problem with airyaiprime and BigFloats:
_ _ _(_)_ | Documentation: https://docs.julialang.org
(_) | (_) (_) |
_ _ _| |_ __ _ | Type "?" for help, "]?" for Pkg help.
| | | | | | |/ _` | |
| | |_| | | | (_| | | Version 1.9.0 (2023-05-07)
_/ |\__'_|_|_|\__'_| | Official https://julialang.org/ release
|__/ |
julia @v1.9> using ArbNumerics
[ Info: Precompiling ArbNumerics [7e558dbc-694d-5a72-987c-6f4ebed21442]
julia @v1.9> ArbNumerics.airyai(BigFloat(0))
0.3550280538878172392600631860041831763979791741991772405833265103008100424501266
julia @v1.9> ArbNumerics.airyaiprime(BigFloat(0))
ERROR: MethodError: no method matching airyaiprime(::Complex{BigFloat})
Closest candidates are:
airyaiprime(::Complex{<:AbstractFloat})
@ SpecialFunctions ~/julia-versions/dot_julia/packages/SpecialFunctions/sKqs4/src/bessel.jl:168
airyaiprime(::Complex)
@ SpecialFunctions ~/julia-versions/dot_julia/packages/SpecialFunctions/sKqs4/src/bessel.jl:165
airyaiprime(::ArbComplex{P}) where P
@ ArbNumerics ~/julia-versions/dot_julia/packages/ArbNumerics/OJfxd/src/float/airy.jl:53
...
Stacktrace:
[1] airyaiprime(z::Complex{BigFloat})
@ SpecialFunctions ~/julia-versions/dot_julia/packages/SpecialFunctions/sKqs4/src/bessel.jl:168
[2] airyaiprime(x::BigFloat)
@ SpecialFunctions ~/julia-versions/dot_julia/packages/SpecialFunctions/sKqs4/src/bessel.jl:173
[3] top-level scope
@ REPL[3]:1
JeffreySarnoff commented
what does this display?
using Pkg; Pkg.status("ArbNumerics")
JeffreySarnoff commented
Actually -- something is wrong.
(while it does not help for Airy, you should prefer to use)
using ArbNumerics
setprecision(ArbFloat, precision(BigFloat))
a = ArbFloat(0.0)
result = function(a)
I need time to look into this.
macd commented
No worries and no rush. Just thought you'd like to know. For the record:
julia @v1.11> using Pkg
julia @v1.11> Pkg.status("ArbNumerics")
Status `~/julia-versions/dot_julia/environments/v1.11/Project.toml`
[7e558dbc] ArbNumerics v1.3.3
julia @v1.11>