首页 | 本学科首页   官方微博 | 高级检索  
     


Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential
Authors:Jacques CR Bloch  Tobias Breu  Simon Heybrock  Tilo Wettig
Affiliation:a Institute for Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
b Department of Mathematics, University of Wuppertal, 42097 Wuppertal, Germany
Abstract:The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace approximation, which uses long recurrences. Even after the deflation of critical eigenvalues, the low efficiency of the method restricts its application to small lattices. Here we propose new short-recurrence methods which strongly enhance the efficiency of the computational method. Using rational approximations to the sign function we introduce two variants, based on the restarted Arnoldi process and on the two-sided Lanczos method, respectively, which become very efficient when combined with multishift solvers. Alternatively, in the variant based on the two-sided Lanczos method the sign function can be evaluated directly. We present numerical results which compare the efficiencies of a restarted Arnoldi-based method and the direct two-sided Lanczos approximation for various lattice sizes. We also show that our new methods gain substantially when combined with deflation.
Keywords:Overlap Dirac operator  Quark chemical potential  Sign function  Non-Hermitian matrix  Iterative methods
本文献已被 ScienceDirect 等数据库收录!
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号