Function sici_with 

pub fn sici_with<P, E, V, const NSI: usize, const NCI: usize, const NP: usize, const NQ: usize>(
    x: V,
    x0: E,
    map: E,
    cutoff: E,
    cheb_si: &[E; NSI],
    cheb_cin: &[E; NCI],
    aux_p: &[E; NP],
    aux_q: &[E; NQ],
) -> (V, V)
where
    E: FloatElement,
    V: FloatVector<Element = E> + SpecializedTranscendentalMath<E>,
    P: Policy,

(Si(x), Ci(x)).

cheb_si and cheb_cin are Chebyshev coefficients for Si(x)/x and Cin(x)/x^2 in v = x^2 mapped onto [-1, 1] by v*map - 1. aux_p and aux_q are ascending monomial coefficients in v = 1/x^2.