Function unwrap_against 

fn unwrap_against(phase: f32, reference: f32) -> f32

phase - round((phase - reference)/(2*PI))*2*PI with mexUnwrap.c’s types (module docs, The unwrap step).