Function angular_dist_pi 

pub fn angular_dist_pi(a: f32, b: f32) -> f32

Smallest angular distance on the circle with period π. Result in [0, π/2].