pub fn warp_complexity(gamma: &[f64], argvals: &[f64]) -> f64
Compute warp complexity: geodesic distance from a warp to the identity.
This is arccos(⟨ψ, ψ_id⟩) on the Hilbert sphere.
arccos(⟨ψ, ψ_id⟩)