Speaker: Abdou Abdel-Rehim (College of William & Mary)
Title: Memory Efficient Algorithm for Non-symmetric Linear
Systems and Eigenvalue Problems with Application to
Lattice QCD
Abstract:
In Lattice QCD, the calculation of the quark propagator requires
the solution of a linear system M*x=b, where M is a large, sparse,
non-symmetric matrix that corresponds to a particular discretization
of the Dirac operator. We often also interested in computing the
quark propagator from many points on the lattice which corresponds
to solving the linear system for multiple right-hand sides b. I will
discuss a new memory efficient algorithm for solving non-symmetric
linear systems and computing eigenvalues based on the Bi-Lanczos
algorithm. The algorithm is based on a similar approach that worked
successfully in the symmetric case.