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.