JuliaMath / JuliaMath/SpecialFunctions.jl

Overflow for besselh(85.5,0.001)

Open
#381 1 comment 0 reactions 0 assignees View on GitHub
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.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.