JeffreySarnoff/ArbNumerics.jl

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

what does this display?
using Pkg; Pkg.status("ArbNumerics")

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>