fn symmetric_eigen(a: &[f64], n: usize) -> (Vec<f64>, Vec<f64>)
Symmetric eigendecomposition via Jacobi rotations. a is row-major n x n (symmetric). Returns (eigenvalues, eigenvectors_col_major) sorted descending.
a
n x n
(eigenvalues, eigenvectors_col_major)