Contour-integral-based rational filter leads to interior eigensolvers for non-Hermitian generalized eigenvalue problems. Based on the Zolotarev's problems, this paper proves the asymptotic optimality of the trapezoidal quadrature of the contour integral in terms of the rational function separation. A composite rule of the trapezoidal quadrature is derived. Two interior eigensolvers are proposed based on the composite rule. Both eigensolvers adopt direct factorization and multi-shift generalized minimal residual method for the inner and outer rational functions, respectively. The first eigensolver fixes the order of the outer rational function and applies the subspace iteration to achieve convergence, whereas the second eigensolver doubles the order of the outer rational function every iteration to achieve convergence without subspace iteration. The efficiency and stability of proposed eigensolvers are demonstrated on synthetic and practical sparse matrix pencils.
翻译:基于围线积分的有理滤波方法可求解非厄米广义特征值问题的内蕴特征值。本文以佐洛塔廖夫问题为基础,从有理函数分离角度证明了围线积分梯形求积的渐近最优性,并推导出梯形求积的复合规则。基于该复合规则提出了两种内蕴特征求解器:两种求解器分别对内层和外层有理函数采用直接分解法与多移位广义最小残差法。第一种求解器固定外层有理函数阶数,通过子空间迭代实现收敛;第二种求解器在每次迭代中将外层有理函数阶数加倍,无需子空间迭代即可收敛。通过合成矩阵和实际稀疏矩阵束验证了所提特征求解器的效率与稳定性。