摘要 Abstract
本文提出了一种有理滤波器,用于计算实轴上位于某一区间内的对称正定特征值问题的所有特征值。嵌入子空间迭代框架中的滤波器所引起的线性系统通过预处理Krylov方法求解。滤波器极点的选择基于两个标准:一方面,滤波器应增强感兴趣的区间内的特征值,这表明极点应选择在区间附近或内部;另一方面,极点的选择对迭代方法的收敛速度有重要影响。对于来自振动问题的问题,这两个标准相互矛盾,因为特征值求解器的快速收敛需要极点位于或靠近区间,而当极点远离特征值时,迭代线性系统求解器变得更便宜。在本文中,我们提出了受Helmholtz方程移位Laplace预处理器启发的极点选择方案,并展示了来自有限元振动模型的数值实验。我们比较了移位Laplace有理滤波器与基于围线积分的高斯求积规则的有理滤波器。
We present a rational filter for computing all eigenvalues of a symmetric definite eigenvalue problem lying in an interval on the real axis. The linear systems arising from the filter embedded in the subspace iteration framework, are solved via a preconditioned Krylov method. The choice of the poles of the filter is based on two criteria. On the one hand, the filter should enhance the eigenvalues in the interval of interest, which suggests that the poles should be chosen close to or in the interval. On the other hand, the choice of poles has an important impact on the convergence speed of the iterative method. For the solution of problems arising from vibrations, the two criteria contradict each other, since fast convergence of the eigensolver requires poles to be in or close to the interval, whereas the iterative linear system solver becomes cheaper when the poles lie further away from the eigenvalues. In the paper, we propose a selection of poles inspired by the shifted Laplace preconditioner for the Helmholtz equation. We show numerical experiments from finite element models of vibrations. We compare the shifted Laplace rational filter with rational filters based on quadrature rules for contour integration.