Michael Helton
Michael Helton
I was hoping to write some specific functions for DoubleFloats using `evalpoly` routines but it appears that using BigFloats is significantly faster... ```julia using DoubleFloats, BenchmarkTools function d64_test(x::Double64) c =...
It would be nice if the example function could also work for `NTuple` types such as ```julia @. model(x, p) = p[1]*exp(-x*p[2]) xdata = (0.5, 1.5, 2.5, 4.5, 6.5) ydata...
I noticed that some of the coefficients could be computed significantly faster using Horner's scheme and combining evaluations to reuse as many computations as possible. This significantly reduces the computation...
# Hyperbola This pull requests implements the numerical inversion of the Laplace transform by using a hyperbola for the Bromwhich contour. It's important to say that the inversion is very...
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...
@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)...
```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...
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...
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...
**EDIT: This issue started with discussion on the `gamma` implementation but has moved to discussing the scope of this package and inclusion of other special functions or splitting into smaller...