LA_STEV and LA_STEVD compute all eigenvalues and,
optionally, all eigenvectors of a real symmetric tridiagonal matrix 
.
LA_STEVD uses a divide and conquer algorithm.
If eigenvectors are desired, it can be much faster than
LA_STEV for large matrices but uses more workspace.