num::unsafe::eig_sym

Defined in header <linear/eigen/jacobi_eig.hpp>

eigen_result eig_sym(const mat<real> &A, real tol = 1e-12, idx max_sweeps = 100);

Symmetric eigendecomposition without requiring or checking the symmetry invariant.

Reads only the lower triangle, so an asymmetric matrix yields the spectrum of its symmetric part rather than an error.