An adaptive homotopy approach for non-selfadjoint eigenvalue problems

This paper presents adaptive algorithms for eigenvalue problems associated with non-selfadjoint partial differential operators. The basis for the developed algorithms is a homotopy method which departs from a well-understood selfadjoint problem. Apart from the adaptive grid refinement, the progress of the homotopy as well as the solution of the iterative method are adapted to balance the contributions of the different error sources. The first algorithm balances the homotopy, discretization and approximation errors with respect to a fixed stepsize τ in the homotopy. The second algorithm combines the adaptive stepsize control for the homotopy with an adaptation in space that ensures an error below a fixed tolerance ε. The outcome of the analysis leads to the third algorithm which allows the complete adaptivity in space, homotopy stepsize as well as the iterative algebraic eigenvalue solver. All three algorithms are compared in numerical examples.

[1]  Merico E. Argentati,et al.  Bounds on Changes in Ritz Values for a Perturbed Invariant Subspace of a Hermitian Matrix , 2008, SIAM J. Matrix Anal. Appl..

[2]  D. Boffi,et al.  Computational Models of Electromagnetic Resonators: Analysis of Edge Element Approximation , 1999 .

[3]  Lloyd N. Trefethen,et al.  Computed eigenmodes of planar regions , 2005 .

[4]  Alexandre Ern,et al.  A Posteriori Control of Modeling Errors and Discretization Errors , 2003, Multiscale Model. Simul..

[5]  Andrew V. Knyazev,et al.  New A Priori FEM Error Estimates for Eigenvalues , 2006, SIAM J. Numer. Anal..

[6]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[7]  Richard B. Lehoucq,et al.  Uniform accuracy of eigenpairs from a shift-invert Lanczos method , 2005, SIAM J. Matrix Anal. Appl..

[8]  Tosio Kato.  Perturbation theory for linear operators , 1966 .

[9]  I. Babuska,et al.  Finite element-galerkin approximation of the eigenvalues and Eigenvectors of selfadjoint problems , 1989 .

[10]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[11]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis: Oden/A Posteriori , 2000 .

[12]  Stefan A. Sauter,et al.  hp-Finite Elements for Elliptic Eigenvalue Problems: Error Estimates Which Are Explicit with Respect to Lambda, h, and p , 2010, SIAM J. Numer. Anal..

[13]  Carsten Carstensen,et al.  An oscillation-free adaptive FEM for symmetric eigenvalue problems , 2011, Numerische Mathematik.

[14]  G. Stewart,et al.  Matrix Perturbation Theory , 1990 .

[15]  T. Y. Li,et al.  Solving eigenvalue problems of real nonsymmetric matrices with real homotopies , 1992 .

[16]  Mats G. Larson,et al.  A Posteriori and a Priori Error Analysis for Finite Element Approximations of Self-Adjoint Elliptic Eigenvalue Problems , 2000, SIAM J. Numer. Anal..

[17]  Endre Süli,et al.  Adaptive finite element methods for differential equations , 2003, Lectures in mathematics.

[18]  H. Keller,et al.  Homotopy Method for the Large, Sparse, Real Nonsymmetric Eigenvalue Problem , 1997, SIAM J. Matrix Anal. Appl..

[19]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[20]  J. M. Thomas,et al.  Introduction à l'analyse numérique des équations aux dérivées partielles , 1983 .

[21]  Eduardo M. Garau,et al.  Convergence of adaptive finite element methods for eigenvalue problems , 2008, 0803.0365.

[22]  Jack Dongarra,et al.  Templates for the Solution of Algebraic Eigenvalue Problems , 2000, Software, environments, tools.

[23]  Luka Grubisic,et al.  On estimators for eigenvalue/eigenvector approximations , 2009, Math. Comput..

[24]  Gene H. Golub,et al.  Matrix computations , 1983 .

[25]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[26]  Stefano Giani,et al.  A Convergent Adaptive Method for Elliptic Eigenvalue Problems , 2009, SIAM J. Numer. Anal..

[27]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[28]  I. Babuška,et al.  Corrigendum: “Finite element-Galerkin approximation of the eigenvalues and eigenvectors of selfadjoint problems” [Math. Comp. 52 (1989), no. 186, 275–297; MR0962210 (89k:65132)] , 1994 .

[29]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[30]  Gene H. Golub,et al.  Homotopy method for the numerical solution of the eigenvalue problem of self-adjoint partial differential operators , 2005, Numerical Algorithms.

[31]  W. Dörfler.  A convergent adaptive algorithm for Poisson's equation , 1996 .

[32]  Klaus Neymeyr,et al.  A posteriori error estimation for elliptic eigenproblems , 2002, Numer. Linear Algebra Appl..

[33]  Volker Mehrmann,et al.  Adaptive computation of smallest eigenvalues of self-adjoint elliptic partial differential equations , 2011, Numer. Linear Algebra Appl..

[34]  Stig Larsson,et al.  Partial differential equations with numerical methods , 2003, Texts in applied mathematics.

[35]  Jinchao Xu,et al.  A two-grid discretization scheme for eigenvalue problems , 2001, Math. Comput..

[36]  Rodolfo Rodríguez,et al.  A Posteriori Error Estimates and a Local Refinement Strategy for a Finite Element Method to Solve Structural-Acoustic Vibration Problems , 2001, Adv. Comput. Math..

[37]  Aihui Zhou,et al.  Adaptive finite element algorithms for eigenvalue problems based on local averaging type a posteriori error estimates , 2006, Adv. Comput. Math..

[38]  Zhonggang Zeng,et al.  The Homotopy Continuation Algorithm for the Real Nonsymmetric Eigenproblem: Further Development and Implementation , 1999, SIAM J. Sci. Comput..

[39]  Merico E. Argentati,et al.  Rayleigh-Ritz Majorization Error Bounds with Applications to FEM , 2009, SIAM J. Matrix Anal. Appl..

[40]  F. Chatelin.  Spectral approximation of linear operators , 2011 .

[41]  Ricardo G. Durán,et al.  A Posteriori Error Estimates for the Finite Element Approximation of Eigenvalue Problems , 2003 .

[42]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis , 2000 .

[43]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[44]  Carsten Carstensen,et al.  Remarks around 50 lines of Matlab: short finite element implementation , 1999, Numerical Algorithms.

[45]  Andrew Knyazev,et al.  New estimates for Ritz vectors , 1997, Math. Comput..

[46]  William H. Press,et al.  Numerical recipes in C , 2002 .

[47]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[48]  Rolf Rannacher,et al.  A posteriori error control for finite element approximations of elliptic eigenvalue problems , 2001, Adv. Comput. Math..

[49]  Rüdiger Verfürth,et al.  A posteriori error estimation and adaptive mesh-refinement techniques , 1994 .

[50]  A. C. Aitken.  XXV.—On Bernoulli's Numerical Solution of Algebraic Equations , 1927 .

[51]  Tien-Yien Li,et al.  Homotopy-determinant algorithm for solving nonsymmetric eigenvalue problems , 1992 .

[52]  Carsten Carstensen,et al.  A Posteriori Error estimators for Non-Symmetric Eigenvalue Problems , 2009 .