Function par_update 

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

f(i, &mut out[i]) for every i (over threads with the parallel feature)