An adaptive homotopy approach for non-selfadjoint eigenvalue problems

C. Carstensen, J. Gedicke, V. Mehrmann, A. Miedlar

Research output: Contribution to journalArticlepeer-review

12 Scopus citations

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 languageEnglish (US)
Pages (from-to)557-583
Number of pages27
JournalNumerische Mathematik
Volume119
Issue number3
DOIs
StatePublished - Nov 2011

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

Fingerprint

Dive into the research topics of 'An adaptive homotopy approach for non-selfadjoint eigenvalue problems'. Together they form a unique fingerprint.

Cite this