Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides

Details Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides

2009-11-12

arxiv.org/abs/0911.2285v1

PoS LAT2009:036,2009

Physics

High Energy Physics

High Energy Physics - Lattice

8 pages, 3 figures, presented at the XXVII International Symposium on Lattice Field Theory, July 26-31, 2009, Peking Universit

Scientific paper

For Hermitian positive definite linear systems and eigenvalue problems, the eigCG algorithm is a memory efficient algorithm that solves the linear system and simultaneously computes some of its eigenvalues. The algorithm is based on the Conjugate-Gradient (CG) algorithm, however, it uses only a window of the vectors generated by the CG algorithm to compute approximate eigenvalues. The number and accuracy of the eigenvectors can be increased by solving more right-hand sides. For Hermitian systems with multiple right-hand sides, the computed eigenvectors can be used to speed up the solution of subsequent systems. The algorithm was tested on Lattice QCD problems by solving the normal equations and was shown to give large speed up factors and to remove the critical slowing down as we approach light quark masses. Here, an extension to the non-symmetric case based on the two-sided Lanczos algorithm is given. The new algorithm is tested on Lattice QCD problems and is shown to give promising results. We also study the removal of the critical slowing down and compare results with those of the eigCG algorithm. We also discuss the case when the system is gamma5-Hermitian.

Affiliated with

Also associated with

No associations

Rating

Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-150071