Function warp_statistics 

pub fn warp_statistics(
    gammas: &FdMatrix,
    argvals: &[f64],
    confidence_level: f64,
) -> Result<WarpStatistics, FdarError>

Compute summary statistics for a set of warping functions.

Given an n x m matrix of warping functions (one per row) and the common evaluation grid, computes pointwise statistics, confidence bands, the Karcher mean warp on the Hilbert sphere, and per-warp geodesic distances.

Arguments

• gammas — Warping functions matrix (n x m), one warp per row.
• argvals — Common evaluation points (length m).
• confidence_level — Confidence level for the bands (e.g. 0.95).

Errors

Returns FdarError::InvalidDimension if n < 2 or m does not match. Returns FdarError::InvalidParameter if confidence_level is not in (0, 1).