pub fn grid_normal(grid: &VoxelGrid, x: i32, y: i32, z: i32) -> Vec3
Estimate normal at a grid point using central differences.