I am trying to link together an eigenvector solver using MAGMA functions, using this standard method:

a symmetric matrix is converted to tridiagonal form

matrix goes through QL transformation

eigenvalues are found, eigenvectors found through reverse iteration (code?), eigenvectors are converted back by the tridiagonal transform in reverse

In my understanding, this would require magma SQEHRD\DGEHRD for the hessenberg form,

followed by SQEQLF/DQEQLF for QL transformation

Do you have sources available for these routines? Am I missing a few steps? Any beta eigenvector codes would also be greatly appreciated.

Best,

Matt Slezak