qimpy.math.spherical_bessel.jl_by_xl

jl_by_xl(l_max, x)

Compute spherical bessel functions j_l(x)/x^l for each l <= l_max. This is optimized to calculate j_l up to l = 6 efficiently combining recursion relations and Taylor expansions to achieve both absolute and relative errors < 1e-15 for all x. (The errors will grow beyond l = 6 due to instability of the efficient recursion relation chosen, so do not use this routine for higher l without testing.)

Parameters: