Function find_close_pairs 

pub fn find_close_pairs(
    positions: &[Vec3],
    max_dist: f32,
) -> Vec<(usize, usize, f32)>

Find all pairs of points closer than max_dist.

Returns (i, j, distance) for each pair where i < j. Uses a spatial grid for O(n log n) performance.