4 ms·
https://github.com/JuliaMath/Bessels.jl/blob/master/src/besselj.jl https://github.com/JuliaMath/Bessels.jl/blob/master/src/bess... Thanks! I love it, so easy t
by nhatcher 4y ago
https://github.com/JuliaMath/Bessels.jl/blob/master/src/besselj.jl https://github.com/JuliaMath/Bessels.jl/blob/master/src/bess...
Thanks! I love it, so easy to understand and follow.
My favourite work on the subject is Fredrik Johansson's:
https://github.com/fredrik-johansson/arb https://github.com/fredrik-johansson/arb
Whenever I feel down and without energy I just read something in there
- adgjlsfhk1 4y agoyeah, Fredrik's papers are great. A lot of the reason Bessels.jl is so good is because ARB makes it easy for us to have test against a known source of truth (which make things like minimax polynomials and picking optimal tradeoffs a ton easier).
- nhatcher 4y agoYou work with Julia? Nice! I actually use Julia to generate my test cases: using Nemo CC = AcbField(100) values = [1, 2, 3, -2, 5, 30, 2e-8] for value in values y_acb = besseli(CC(1), CC(value)) real64 = convert(Float64, real(y_acb)) im64 = convert(Float64, real(y_acb)) println("(", value, ", ", real64, "),") end