Is your feature request related to a problem? Please describe. Eigendecomposition of symmetric 2x2 and 3x3 matrices (scipy.linalg.eigh) is currently performed via LAPACK's general algorithms (various ...