accelerating the reduction to upper hessenberg , tridiagonal , and bidiagonal forms through hybrid gpu - based computing .
matrix computations ( 3 . ed . ) .