pub fn cartesian_gto_cint_batch_rayon(
    a: f64,
    l: usize,
    c: &[f64; 3],
    r: &[[f64; 3]]
) -> MatrixFull<f64>