Static J1_COEFFS_TAYLOR 

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

J1 taylor expansion at zeros and extremums.

Generated by SageMath:

def print_taylor_coeffs(poly, x0):
    print("[")
    for i in range(0, 24):
        coeff = poly[i]
        print_double_double("", coeff)
        # print(f"{double_to_hex(coeff)},")
    print("],")

print(f"pub(crate) static J1_COEFFS_TAYLOR: [(u64, u64); {len(j1_zeros)}] = [")

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

    x0 = mp.mpf(k_range)
    from mpmath import mp, j1, taylor
    poly = taylor(lambda val: j1(val), x0, 25)
    print_taylor_coeffs(poly, CDR(k_range))
    prev_zero = j1_zeros[i]

print("];")