Static J0_COEFFS_TAYLOR 

pub(crate) static J0_COEFFS_TAYLOR: [[(u64, u64); 24]; 47]

Taylor expansion at zeros and extremums for J0.

Generated by SageMath:

mp.prec = 180
def print_taylor_coeffs(poly):
    print("[")
    for i in range(0, 24):
        coeff = poly[i]
        print_double_double("", coeff)
    print("],")

def print_taylor_coeffs_dyad(poly):
    print("[")
    for i in range(0, 24):
        coeff = poly[i]
        print_dyadic(coeff)
    print("],")

print(f"pub(crate) static J0_COEFFS_TAYLOR: [[(u64, u64); 24]; {len(j0_zeros)}] = [")

prev_zero = 0

for i in range(0, len(j0_zeros)):
    k_range = j0_zeros[i]
    range_diff = k_range - prev_zero

    mp.prec = 180
    x0 = mp.mpf(k_range)
    from mpmath import mp, j0, taylor
    poly = taylor(lambda val: j0(val), x0, 24)
    print_taylor_coeffs(poly)
    # print(poly)
    prev_zero = j0_zeros[i]

print("];")