pub fn split_with_halo( global: &Array2<f64>, slices: &[DomainSlice], ) -> FusionResult<Vec<Array2<f64>>>