fn correct_regions(phase: &mut [f64], visited: &[u8], num_regions: u8)
Correct each region’s median to nearest 0 by subtracting n·2π.
Matches ROMEO.jl correct_regions!.
correct_regions!