×

Solving large-scale nonlinear eigenvalue problems by rational interpolation and resolvent sampling based Rayleigh-Ritz method. (English) Zbl 1367.65162

Summary: Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by the matrix structures, properties of the eigen-solutions, sizes of the problems, etc. This paper aims to remove those limitations and develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach, the RIA provides the possibility to select sampling points in more general regions and has advantages in improving the accuracy and reducing the computational cost. Second, a resolvent sampling scheme using the RIA is proposed to construct reliable search spaces for the Rayleigh-Ritz procedure, based on which a robust eigen-solver, called resolvent sampling based Rayleigh-Ritz method (RSRR), is developed for solving general NEPs. The RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of the RSRR are demonstrated by a variety of benchmark and application examples.

MSC:

65F15 Numerical computation of eigenvalues and eigenvectors of matrices
65H17 Numerical solution of nonlinear eigenvalue and eigenvector problems

Software:

NLEVP; FEAST; NLEIGS
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Tisseur, The quadratic eigenvalue problem, SIAM Review 43 (2) pp 235– (2001) · Zbl 0985.65028
[2] Mehrmann, Nonlinear eigenvalue and frequency response problems in industrial practice, Journal of Mathematics in Industry 1 (1) pp 1– (2011) · Zbl 1273.65058
[3] Effenberger C Robust solution methods for nonlinear eigenvalue problems Ph.D. Thesis 2013
[4] Van Beeumen R Rational Krylov methods for nonlinear eigenvalue problems Ph.D. Thesis 2015
[5] Opstal, A finite-element/boundary-element method for three-dimensional, large-displacement fluid-structure-interaction, Computer Methods in Applied Mechanics and Engineering 284 pp 637– (2015) · Zbl 1423.74287
[6] Kimeswenger, Coupled finite and boundary element methods for fluid-solid interaction eigenvalue problems, SIAM Journal on Numerical Analysis 52 (5) pp 2400– (2014) · Zbl 1307.74032
[7] Abdel-Aziz, Safeguarded use of the implicit restarted Lanczos technique for solving non-linear structural eigensystems, International Journal for Numerical Methods in Engineering 37 (18) pp 3117– (1994) · Zbl 0812.73073
[8] Dumont, On the solution of generalized non-linear complex-symmetric eigenvalue problems, International Journal for Numerical Methods in Engineering 71 (13) pp 1534– (2007) · Zbl 1194.74133
[9] Effenberger, Robust successive computation of eigenpairs for nonlinear eigenvalue problems, SIAM Journal on Matrix Analysis and Applications 34 (3) pp 1231– (2013) · Zbl 1297.47067
[10] Salas, Spectral recycling strategies for the solution of nonlinear eigenproblems in thermoacoustics, Numerical Linear Algebra with Applications 22 (6) pp 1039– (2015) · Zbl 1374.65090
[11] Lu, A Padé approximate linearization algorithm for solving the quadratic eigenvalue problem with low-rank damping, International Journal for Numerical Methods in Engineering 103 (11) pp 840– (2015) · Zbl 1352.65120
[12] Sakurai, A projection method for generalized eigenvalue problems using numerical integration, Journal of Computational and Applied Mathematics 159 (1) pp 119– (2003) · Zbl 1037.65040
[13] Polizzi, Density-matrix-based algorithm for solving eigenvalue problems, Physical Review B 79 (11) pp 1– (2009)
[14] Asakura, A numerical method for nonlinear eigenvalue problems using contour integrals, JSIAM Letters 1 (0) pp 52– (2009) · Zbl 1278.65072
[15] Beyn, An integral method for solving nonlinear eigenvalue problems, Linear Algebra and Its Applications 436 (10) pp 3839– (2012) · Zbl 1237.65035
[16] Yokota, A projection method for nonlinear eigenvalue problems using contour integrals, JSIAM Letters 5 (0) pp 41– (2013)
[17] Gavin, Non-linear eigensolver-based alternative to traditional SCF methods, The Journal of Chemical Physics 138 (19) pp 1– (2013)
[18] Austin, Computing eigenvalues of real symmetric matrices with rational filters in real arithmetic, SIAM Journal on Scientific Computing 37 (3) pp A1365– (2015) · Zbl 1328.15016
[19] Güttel, Zolotarev quadrature rules and load balancing for the FEAST eigensolver, SIAM Journal on Scientific Computing 37 (4) pp A2100– (2015) · Zbl 1321.65055
[20] Sakurai, Error analysis for a matrix pencil of Hankel matrices with perturbed complex moments, JSIAM Letters 1 (0) pp 76– (2009) · Zbl 1278.65041
[21] Van Beeumen, A rational Krylov method based on Hermite interpolation for nonlinear eigenvalue problems, SIAM Journal on Scientific Computing 35 (1) pp A327– (2013) · Zbl 1284.47041
[22] Neumaier, Residual inverse iteration for the nonlinear eigenvalue problem, SIAM Journal on Numerical Analysis 22 (5) pp 914– (1985) · Zbl 0594.65026
[23] Voss, A Jacobi-Davidson method for nonlinear and nonsymmetric eigenproblems, Computers & Structures 85 (17) pp 1284– (2007)
[24] Kressner, A block Newton method for nonlinear eigenvalue problems, Numerische Mathematik 114 (2) pp 355– (2009) · Zbl 1191.65054
[25] Hochstenbach, Controlling inner iterations in the Jacobi-Davidson method, SIAM Journal on Matrix Analysis and Applications 31 (2) pp 460– (2009) · Zbl 1191.65032
[26] Effenberger, Linearization techniques for band structure calculations in absorbing photonic crystals, International Journal for Numerical Methods in Engineering 89 (2) pp 180– (2012) · Zbl 1242.74110
[27] Kirkup, Solution of the Helmholtz eigenvalue problem via the boundary element method, International Journal for Numerical Methods in Engineering 36 (2) pp 321– (1993) · Zbl 0825.76470
[28] Ali, Advances in acoustic eigenvalue analysis using boundary element method, Computers & Structures 56 (5) pp 837– (1995) · Zbl 1002.76539
[29] Effenberger, Chebyshev interpolation for nonlinear eigenvalue problems, BIT Numerical Mathematics 52 (4) pp 933– (2012) · Zbl 1263.65048
[30] Van Beeumen, Compact rational Krylov methods for nonlinear eigenvalue problems, SIAM Journal on Matrix Analysis and Applications 36 (2) pp 820– (2015) · Zbl 1319.65042
[31] Effenberger, Interpolation-based solution of a nonlinear eigenvalue problem in fluid-structure interaction, PAMM 12 (1) pp 633– (2012)
[32] SolovŠëv, Preconditioned iterative methods for a class of nonlinear eigenvalue problems, Linear Algebra and Its Applications 415 (1) pp 210– (2006) · Zbl 1095.65033
[33] Liao, Nonlinear Rayleigh-Ritz iterative method for solving large scale nonlinear eigenvalue problems, Taiwanese Journal of Mathematics 14 (3A) pp 869– (2010)
[34] Cao, A fast directional BEM for large-scale acoustic problems based on the Burton-Miller formulation, Engineering Analysis with Boundary Elements 50 pp 47– (2015) · Zbl 1403.76062
[35] Nicoud, Acoustic modes in combustors with complex impedances and multidimensional active flames, AIAA journal 45 (2) pp 426– (2007)
[36] Van Barel, Designing rational filter functions for solving eigenvalue problems by contour integration, Linear Algebra and its Applications 502 pp 346– (2016) · Zbl 1386.65115
[37] Austin, Numerical algorithms based on analytic function values at roots of unity, SIAM Journal on Numerical Analysis 52 (4) pp 1795– (2014) · Zbl 1305.41006
[38] Saff, An extension of Montessus de Ballore’s theorem on the convergence of interpolating rational functions, Journal of Approximation Theory 6 (1) pp 63– (1972) · Zbl 0241.30013
[39] Eğecioğlu, A fast algorithm for rational interpolation via orthogonal polynomials, Mathematics of Computation 53 (187) pp 249– (1989)
[40] Berrut, Barycentric Lagrange interpolation, SIAM Review 46 (3) pp 501– (2004) · Zbl 1061.65006
[41] Sakurai, Efficient parameter estimation and implementation of a contour integral-based eigensolver, Journal of Algorithms & Computational Technology 7 (3) pp 249– (2013)
[42] Wang, Explicit barycentric weights for polynomial interpolation in the roots or extrema of classical orthogonal polynomials, Mathematics of Computation 83 (290) pp 2893– (2014) · Zbl 1297.41001
[43] Parks, Recycling Krylov subspaces for sequences of linear systems, SIAM Journal on Scientific Computing 28 (5) pp 1651– (2006) · Zbl 1123.65022
[44] Hasegawa, Recovering from accuracy deterioration in the contour integral-based eigensolver, JSIAM Letters 8 (0) pp 1– (2015)
[45] Betcke, NLEVP: a collection of nonlinear eigenvalue problems, ACM Transactions on Mathematical Software (TOMS) 39 (2) pp 1– (2013) · Zbl 1295.65140
[46] Güttel, NLEIGS: a class of fully rational Krylov methods for nonlinear eigenvalue problems, SIAM Journal on Scientific Computing 36 (6) pp A2842– (2014) · Zbl 1321.47128
[47] Xiao, Resolvent sampling based Rayleigh-Ritz method for large-scale nonlinear eigenvalue problems, Computer Methods in Applied Mechanics and Engineering 310 pp 33– (2016)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.