pub fn compute_rmsd(coords: &[f64], reference: &[f64]) -> f64
Compute RMSD between two coordinate sets after Kabsch alignment.