We describe a novel algorithm for solving general parametric (nonlinear) eigenvalue problems. Our method has two steps: first, high-accuracy solutions of non-parametric versions of the problem are gathered at some values of the parameters; these are then combined to obtain global approximations of the parametric eigenvalues. To gather the non-parametric data, we use non-intrusive contour-integration-based methods, which, however, cannot track eigenvalues that migrate into/out of the contour as the parameter changes. Special strategies are described for performing the combination-over-parameter step despite having only partial information on such migrating eigenvalues. Moreover, we dedicate a special focus to the approximation of eigenvalues that undergo bifurcations. Finally, we propose an adaptive strategy that allows one to effectively apply our method even without any a priori information on the behavior of the sought-after eigenvalues. Numerical tests are performed, showing that our algorithm can achieve remarkably high approximation accuracy.
翻译:我们描述了一种求解一般参数(非线性)特征值问题的新型算法。该方法包含两个步骤:首先,在参数某些取值处收集参数无关版本问题的高精度解;随后将这些解组合以获得参数特征值的全局近似。为收集参数无关数据,我们采用非侵入式基于围道积分的方法,但当参数变化导致特征值进出围道时,该方法无法追踪此类迁移特征值。针对仅有此类迁移特征值的部分信息的情况,我们提出了特殊策略执行参数组合步骤。此外,我们特别关注了经历分岔的特征值的近似问题。最后,我们提出一种自适应策略,使得即使没有关于目标特征值行为的先验信息,也能有效应用该方法。数值实验表明,该算法能够达到显著高的近似精度。