Function peak_localisation_error 

pub fn peak_localisation_error(resolution: &Array2<f64>) -> Array1<usize>

Compute peak localisation error for each source.

For each source i, finds the index of the maximum in the PSF and reports the index offset |argmax(R[:, i]) − i|.

Returns an array of length n_sources.