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.
翻译:我们描述了一种求解通用参数(非线性)特征值问题的新算法。该方法包含两个步骤:首先,在参数某些取值处收集非参数版本问题的高精度解;随后将这些解组合以获得参数特征值的全局近似。为收集非参数数据,我们采用非侵入式基于轮廓积分的方法,但该方法无法追踪随参数变化而移入/移出轮廓的特征值。本文针对此类迁移特征值仅具有部分信息的情况,提出了执行参数组合步骤的特殊策略。此外,我们特别关注了发生分岔的特征值近似问题。最后,我们提出一种自适应策略,使得即使缺乏目标特征值行为的先验信息,也能有效应用该方法。数值实验表明,我们的算法能够实现极高的近似精度。