Function wasserstein_drift

pub fn wasserstein_drift(
    dist_t1: &[f32],
    dist_t2: &[f32],
    centroids: &[&[f32]],
    n_projections: usize,
) -> f64

Compute Wasserstein drift between two region distributions.

Convenience function that wraps sliced_wasserstein for the common use case of comparing region distributions at two time points.

Arguments

• dist_t1 - Region distribution at time T₁ (K floats summing to ~1.0).
• dist_t2 - Region distribution at time T₂.
• centroids - Region centroid vectors from index.regions(level).
• n_projections - Number of projections (default: 50).