Function pairwise_consistency 

pub fn pairwise_consistency(
    data: &FdMatrix,
    argvals: &[f64],
    lambda: f64,
    max_triplets: usize,
) -> f64

Measure pairwise alignment consistency via triplet checks.

For triplets (i,j,k), checks γ_ij ∘ γ_jk ≈ γ_ik by measuring the L2 deviation of the composed warp from the direct warp.

Arguments

• data — Functional data (n × m)
• argvals — Evaluation points (length m)
• lambda — Penalty weight
• max_triplets — Maximum number of triplets to check (0 = all)