A Hamiltonian Krylov-Schur-type method based on the symplectic Lanczos process

被引：19

作者：

Benner, Peter ^{[1
]}

Fassbender, Heike ^{[2
]}

Stoll, Martin ^{[3
]}

机构：

[1] TU Chemnitz, Fak Math Math Ind & Tech, D-09107 Chemnitz, Germany

[2] TU Braunschweig, AG Numer, Inst Computat Math, D-38092 Braunschweig, Germany

[3] Math Inst, Oxford Ctr Collaborat Appl Mat, Oxford OX1 3LB, England

来源：

LINEAR ALGEBRA AND ITS APPLICATIONS | 2011年 / 435卷 / 03期

关键词：

Hamiltonian eigenproblem; Symplectic Lanczos method; Krylov-Schur method; Implicit restarting; SR algorithm; MODEL-REDUCTION; ALGORITHM; MATRIX;

D O I：

10.1016/j.laa.2010.04.048

中图分类号：

O29 [应用数学];

学科分类号：

070104 ;

摘要：

We discuss a Krylov-Schur-like restarting technique applied within the symplectic Lanczos algorithm for the Hamiltonian eigenvalue problem. This allows us to easily implement a purging and locking strategy in order to improve the convergence properties of the symplectic Lanczos algorithm. The Krylov-Schur-like restarting is based on the SR algorithm. Some ingredients of the latter need to be adapted to the structure of the symplectic Lanczos recursion. We demonstrate the efficiency of the new method for several Hamiltonian eigenproblems. (C) 2010 Elsevier Inc. All rights reserved.

引用

页码：578 / 600

页数：23

共 58 条

[21] Algorithm 854: Fortran 77 subroutines for computing the eigenvalues of Hamiltonian matrices II [J].

Benner, Peter ;

Kressner, Daniel .

ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2006, 32 (02) :352-373

[22]

Benner P, 2007, ELECTRON T NUMER ANA, V29, P212

[23]

Boyd S., 1989, Mathematics of Control, Signals, and Systems, V2, P207, DOI 10.1007/BF02551385

[24] A SYMPLECTIC QR LIKE ALGORITHM FOR THE SOLUTION OF THE REAL ALGEBRAIC RICCATI EQUATION [J].

BUNSEGERSTNER, A ;

MEHRMANN, V .

IEEE TRANSACTIONS ON AUTOMATIC CONTROL, 1986, 31 (12) :1104-1113

[25] MATRIX FACTORIZATIONS FOR SYMPLECTIC QR-LIKE METHODS [J].

BUNSEGERSTNER, A .

LINEAR ALGEBRA AND ITS APPLICATIONS, 1986, 83 :49-77

[26] Robust stability and a criss-cross algorithm for pseudospectra [J].

Burke, JV ;

Lewis, AS ;

Overton, ML .

IMA JOURNAL OF NUMERICAL ANALYSIS, 2003, 23 (03) :359-375

[27] A BISECTION METHOD FOR MEASURING THE DISTANCE OF A STABLE MATRIX TO THE UNSTABLE MATRICES [J].

BYERS, R .

SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1988, 9 (05) :875-881

[28]

Fassbender H, 2007, ELECTRON T NUMER ANA, V26, P121

[29] The shift-inverted J-Lanczos algorithm for the numerical solutions of large sparse algebraic Riccati equations [J].

Ferng, WR ;

Lin, WW ;

Wang, CS .

COMPUTERS & MATHEMATICS WITH APPLICATIONS, 1997, 33 (10) :23-40

[30] QR TRANSFORMATION - A UNITARY ANALOG TO LR TRANSFORMATION .1. [J].

FRANCIS, J .

COMPUTER JOURNAL, 1961, 4 :265-&

← 1 2 3 4 5 6 →