My problem is to compute eigenvectors and eigenvalues of a lot of small (n < 30) symetric, positive definite matrices.
So far I am using LAPACK's DSYEV.
The priority is speed more than accuracy. Is there a faster algorithm (maybe not asymptotically, but in practice for small input)?