Abstract
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.
Original language | English (US) |
---|---|
Pages (from-to) | 557-583 |
Number of pages | 27 |
Journal | Numerische Mathematik |
Volume | 119 |
Issue number | 3 |
DOIs | |
State | Published - Nov 2011 |
Externally published | Yes |
Bibliographical note
Funding Information:V. Mehrmann: Supported by the DFG Research Center Matheon “Mathematics for Key Technologies” in Berlin. J. Gedicke and A. Miedlar: Supported by the DFG graduate school BMS “Berlin Mathematical School” in Berlin. C. Carstensen: World Class University (WCU) program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology R31-2008-000-10049-0.
Keywords
- 65F15
- 65H20
- 65M60
- 65N15
- 65N25
- 65N50