JuliaMath / JuliaMath/Bessels.jl

Improve accuracy of asymptotic expansion for very large arguments (besselj, bessely)

Open
#37 0 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
90
Forks
11
PR merge metrics
No merged PRs in 30d

Description

Right now, for very large arguments there are some roundoff errors when the argument is very large.

```julia
julia> Bessels.besseljy_large_argument(11, 2e9)
(-6.441564914025336e-6, 1.663779215046842e-5)

julia> Bessels.besseljy_large_argument(big"11", big"2e9")
(-6.441565416278640966871190587378818114705713999893782928810537933625658595387592e-06, 1.663779195601368418718142484679994485342917679346562260323801403887512758594912e-05)

```

It would be nice to improve the accuracy of these expansions.

Edit: Though it does look like it then improves...
```
julia> SpecialFunctions.besselj(11.0, 2e9)
-6.441565416278642e-6

julia> besselj(11.0, 2e9)
-6.441564914025336e-6

julia> besselj(11.0, 2e20)
4.600631613945902e-11

julia> SpecialFunctions.besselj(11.0, 2e20)
4.600631613945894e-11
```

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.