Function par_fill 

pub(crate) fn par_fill<T: Send>(
    out: &mut [T],
    f: impl Fn(usize) -> T + Sync + Send,
)

out[i] = f(i) for every i (over threads with the parallel feature)