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 #88 (comment)

The fix should be with the better sin_sum PR. Need to test further.