pub fn warp_smoothness(gamma: &[f64], argvals: &[f64]) -> f64
Compute warp smoothness (bending energy): ∫(γ’’)² dt.