Bessels.jl

Bessels.jl copied to clipboard

Consider calling `besselj(n, x)` where `x` is negative, while `n` is a large negative integer: ```julia julia> Bessels.besselj(-172, -2.0) NaN julia> Bessels.besselj(-37, -2.0f0) NaN32 ``` When at least one argument...

yakovbraver

Bumps [julia-actions/setup-julia](https://github.com/julia-actions/setup-julia) from 1 to 2. Release notes Sourced from julia-actions/setup-julia's releases. v2.0.0 - Update to Node20 What's Changed update to node20 by @​ranocha in julia-actions/setup-julia#209 Note the rationale for...

dependabot[bot]

Use better sin_sum for F32

6

This fixes #90 where performance was fixed in #92. ```julia # before julia> Bessels.besselj0(328049.34f0) -0.0013240778f0 # after julia> Bessels.besselj0(328049.34f0) -0.0013258625f0 # Float64 number julia> Bessels.besselj0(Float64(328049.34f0)) -0.001325862383187567 ``` This significantly improves...

heltonmc

Hankel functions sometimes throw NaN values

1

I recently bumped into hankelh1() and hankelh2() returning NaN real parts for nu = 7, x= \sqrt(15^2-1). Values close to this point both in nu and x behave nicely. Is...

olaca09

Support for BigFloat arguments

4

Why cant we get arbitrary precision on the besselj functions?

fda-tome

A different approach to compute besselk in intermediate zones

33

@cgeoga I was thinking some more about how we compute the region between asymptotic expansions (power series, uniform expansions, and large argument). Which corresponds to roughly `0 @benchmark test(2.2, x)...

heltonmc

```julia julia> Bessels.besselj0(328049.34f0) -0.0013240778f0 julia> Bessels.besselj0(Float64(328049.34f0)) -0.0013258623831875669 ``` I'm imagining the better sin sum here could help. _Originally posted by @heltonmc in https://github.com/JuliaMath/Bessels.jl/pull/88#discussion_r1162007124_ The fix should be with the better...

heltonmc

The test structure is fairly comprehensive but we are starting to run into a few things as the library grows. The test files are becoming a conglomeration of new functions...

heltonmc

Variable order in large-order expansions (at least in `besselk`)

11

I finally got around to cleaning up the uniform asymptotic expansion code in `BesselK.jl` and I ended up with a solution that is within a few `ns` of `Bessels.besselk_large_orders` while...

cgeoga

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...

heltonmc