Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Problem with airyaiprime #61

Open
macd opened this issue Aug 5, 2023 · 3 comments
Open

Problem with airyaiprime #61

macd opened this issue Aug 5, 2023 · 3 comments

Comments

@macd
Copy link

macd commented Aug 5, 2023

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
Copy link
Owner

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

@JeffreySarnoff
Copy link
Owner

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
Copy link
Author

macd commented Aug 7, 2023

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> 

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants