34
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Computing and deflating eigenvalues while solving multiple right hand side linear systems in Quantum Chromodynamics

      Preprint
      ,

      Read this article at

      ScienceOpenPublisherArXiv
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          We present a new algorithm that computes eigenvalues and eigenvectors of a Hermitian positive definite matrix while solving a linear system of equations with Conjugate Gradient (CG). Traditionally, all the CG iteration vectors could be saved and recombined through the eigenvectors of the tridiagonal projection matrix, which is equivalent theoretically to unrestarted Lanczos. Our algorithm capitalizes on the iteration vectors produced by CG to update only a small window of vectors that approximate the eigenvectors. While this window is restarted in a locally optimal way, the CG algorithm for the linear system is unaffected. Yet, in all our experiments, this small window converges to the required eigenvectors at a rate identical to unrestarted Lanczos. After the solution of the linear system, eigenvectors that have not accurately converged can be improved in an incremental fashion by solving additional linear systems. In this case, eigenvectors identified in earlier systems can be used to deflate, and thus accelerate, the convergence of subsequent systems. We have used this algorithm with excellent results in lattice QCD applications, where hundreds of right hand sides may be needed. Specifically, about 70 eigenvectors are obtained to full accuracy after solving 24 right hand sides. Deflating these from the large number of subsequent right hand sides removes the dreaded critical slowdown, where the conditioning of the matrix increases as the quark mass reaches a critical value. Our experiments show almost a constant number of iterations for our method, regardless of quark mass, and speedups of 8 over original CG for light quark masses.

          Related collections

          Most cited references28

          • Record: found
          • Abstract: not found
          • Article: not found

          Confinement of quarks

            • Record: found
            • Abstract: not found
            • Article: not found

            Methods of conjugate gradients for solving linear systems

              • Record: found
              • Abstract: not found
              • Article: not found

              The Geometry of Algorithms with Orthogonality Constraints

                Author and article information

                Journal
                01 July 2007
                2008-06-11
                Article
                10.1137/080725532
                0707.0131
                51f1bc0c-c4e6-44e3-8b33-926e9d171fe8

                http://arxiv.org/licenses/nonexclusive-distrib/1.0/

                History
                Custom metadata
                SIAM J. Sci. Comput. Vol. 32, No. 1, 439--462, 2010
                22 pages, 26 eps figures
                hep-lat

                Comments

                Comment on this article

                Related Documents Log