Function phase_distance_pair 

pub fn phase_distance_pair(
    f1: &[f64],
    f2: &[f64],
    argvals: &[f64],
    lambda: f64,
) -> f64

Compute the phase distance between two curves (geodesic distance of optimal warp from identity).