Function euclidean_cluster_roots_grid 

pub fn euclidean_cluster_roots_grid(
    x: &[f32],
    y: &[f32],
    z: &[f32],
    cluster_tolerance: f32,
) -> Result<Vec<u32>, SpatialError>

Connected-component roots via uniform-grid union-find (minimum index per component).