Adding possibility of not computing eigenvectors: this should speed things up by a factor 3 in the QR factorization, if I'm right.