JuliaMath / JuliaMath/SpecialFunctions.jl
Overflow for besselh(85.5,0.001)
- Dominant language
- Julia
- Stars
- 381
- Forks
- 113
- Avg merge
- 1d 23h
- Merged PRs (30d)
- 1
Description
The following code gives an overflow. Is there a way to support BigFloat or could it return infinity instead of an exception?
```
julia> besselh(85.5,0.001)
ERROR: AmosException with id 2: overflow.
Stacktrace:
[1] _bessely(nu::Float64, z::ComplexF64, kode::Int32)
@ SpecialFunctions ~/.julia/packages/SpecialFunctions/9pXme/src/bessel.jl:292
[2] bessely(nu::Float64, z::ComplexF64)
@ SpecialFunctions ~/.julia/packages/SpecialFunctions/9pXme/src/bessel.jl:429
[3] bessely(nu::Float64, x::Float64)
@ SpecialFunctions ~/.julia/packages/SpecialFunctions/9pXme/src/bessel.jl:560
[4] besselh(nu::Float64, k::Int64, x::Float64)
@ SpecialFunctions ~/.julia/packages/SpecialFunctions/9pXme/src/bessel.jl:320
[5] besselh(nu::Float64, z::Float64)
@ SpecialFunctions ~/.julia/packages/SpecialFunctions/9pXme/src/bessel.jl:597
[6] top-level scope
@ REPL[12]:1
julia>
```
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.