pub fn get_default_iterations(res: (f32, f32, f32), step_size: f32) -> usize
Get default number of iterations based on voxel size and step size.
Matches the Julia reference: max(1000, 3200 / prod(res)^0.42) / step_size^0.6
max(1000, 3200 / prod(res)^0.42) / step_size^0.6