18-24 June 2017
Palacio de Congresos
Europe/Madrid timezone
Home > Timetable > Session details > Contribution details

Contribution Parallel

Seminarios 8
Algorithms and Machines

A staggered eigensolver based on sparse matrix bidiagonalization


  • James OSBORN

Primary authors



I will present a method for calculating eigenvectors of the staggered Dirac operator based on the Golub-Kahan-Lanczos bidiagonalization algorithm. Instead of using orthogonalization during the bidiagonalization procedure to increase stability, we choose to stabilize the method by combining it with an outer iteration that refines the approximate eigenvectors obtained from the inner bidiagonalization procedure. I will discuss the performance of the current implementation using QEX and compare with other methods.