### Abstract

We describe a matrix diagonalization algorithm for complex symmetric (not Hermitian) matrices, A̲=A̲^{T}, which is based on a two-step algorithm involving generalized Householder reflections based on the indefinite inner product 〈u̲,v̲〉_{∗}=∑_{i}u_{i}v_{i}. This inner product is linear in both arguments and avoids complex conjugation. The complex symmetric input matrix is transformed to tridiagonal form using generalized Householder transformations (first step). An iterative, generalized QL decomposition of the tridiagonal matrix employing an implicit shift converges toward diagonal form (second step). The QL algorithm employs iterative deflation techniques when a machine-precision zero is encountered “prematurely” on the super-/sub-diagonal. The algorithm allows for a reliable and computationally efficient computation of resonance and antiresonance energies which emerge from complex-scaled Hamiltonians, and for the numerical determination of the real energy eigenvalues of pseudo-Hermitian and PT-symmetric Hamilton matrices. Numerical reference values are provided. Program summary Program Title: HTDQLS Program Files doi: http://dx.doi.org/10.17632/x24wjxtrsg.1 Licensing provisions: GPLv3 Programming language: Fortran 90 using fixed form notation Nature of problem: Calculating the eigenvalues and optionally the eigenvectors of complex symmetric (non-Hermitian), densely populated matrices. Solution method: The complex symmetric (not Hermitian) input matrix is diagonalized in two steps. First step: The matrix is tridiagonalized via a series of (n−2) generalized Householder reflections, where n is the rank of the input matrix. Second step: The tridiagonal matrix is diagonalized via a generalization of the “chasing the bulge” technique, which is an iterative process utilizing an implicitly shifted initial rotation followed by (n−2) Givens rotations. This technique is an implementation of QL factorization, and converges roughly as [(λ_{i}−σ_{i})∕(λ_{i+1}−σ_{i})]^{j} where λ_{i} is the eigenvalue located in the (i,i) position of the final diagonal matrix and the eigenvalues are ordered (|λ_{1}|<|λ_{2}|<…<|λ_{n}|), and j is the iteration. The “educated guess” σ_{i} for the eigenvalue λ_{i} is obtained from the analytic determination of the eigenvalues of (k×k)-submatrices of A̲, in the vicinity of the ith element, where k=0,1,2,3 (here, k=0 means that the implicit shift vanishes, σ_{i}=0). The routine optionally calculates the rotation matrix Z̲, such that Z̲^{−1}A̲Z̲=D̲ where A̲ is the input matrix and D̲ is the diagonal matrix containing the eigenvalues. The ith column of Z̲ then is the eigenvector of A̲ corresponding to the eigenvalue found at the element D̲(i,i), in the ith position on the diagonal of the matrix D̲. Unusual features: For simplicity, the “wrapper” program which contains an example application and the HTDQLS routine are distributed in the same file.

Original language | English |
---|---|

Pages (from-to) | 304-316 |

Number of pages | 13 |

Journal | Computer Physics Communications |

Volume | 221 |

DOIs | |

Publication status | Published - Dec 1 2017 |

### Fingerprint

### Keywords

- Complex symmetric matrix diagonalization
- Deflation techniques
- Implicit shift
- Indefinite inner product

### ASJC Scopus subject areas

- Hardware and Architecture
- Physics and Astronomy(all)

### Cite this

*Computer Physics Communications*,

*221*, 304-316. https://doi.org/10.1016/j.cpc.2017.06.014

**Diagonalization of complex symmetric matrices : Generalized Householder reflections, iterative deflation and implicit shifts.** / Noble, J. H.; Lubasch, M.; Stevens, J.; Jentschura, U.

Research output: Contribution to journal › Article

*Computer Physics Communications*, vol. 221, pp. 304-316. https://doi.org/10.1016/j.cpc.2017.06.014

}

TY - JOUR

T1 - Diagonalization of complex symmetric matrices

T2 - Generalized Householder reflections, iterative deflation and implicit shifts

AU - Noble, J. H.

AU - Lubasch, M.

AU - Stevens, J.

AU - Jentschura, U.

PY - 2017/12/1

Y1 - 2017/12/1

N2 - We describe a matrix diagonalization algorithm for complex symmetric (not Hermitian) matrices, A̲=A̲T, which is based on a two-step algorithm involving generalized Householder reflections based on the indefinite inner product 〈u̲,v̲〉∗=∑iuivi. This inner product is linear in both arguments and avoids complex conjugation. The complex symmetric input matrix is transformed to tridiagonal form using generalized Householder transformations (first step). An iterative, generalized QL decomposition of the tridiagonal matrix employing an implicit shift converges toward diagonal form (second step). The QL algorithm employs iterative deflation techniques when a machine-precision zero is encountered “prematurely” on the super-/sub-diagonal. The algorithm allows for a reliable and computationally efficient computation of resonance and antiresonance energies which emerge from complex-scaled Hamiltonians, and for the numerical determination of the real energy eigenvalues of pseudo-Hermitian and PT-symmetric Hamilton matrices. Numerical reference values are provided. Program summary Program Title: HTDQLS Program Files doi: http://dx.doi.org/10.17632/x24wjxtrsg.1 Licensing provisions: GPLv3 Programming language: Fortran 90 using fixed form notation Nature of problem: Calculating the eigenvalues and optionally the eigenvectors of complex symmetric (non-Hermitian), densely populated matrices. Solution method: The complex symmetric (not Hermitian) input matrix is diagonalized in two steps. First step: The matrix is tridiagonalized via a series of (n−2) generalized Householder reflections, where n is the rank of the input matrix. Second step: The tridiagonal matrix is diagonalized via a generalization of the “chasing the bulge” technique, which is an iterative process utilizing an implicitly shifted initial rotation followed by (n−2) Givens rotations. This technique is an implementation of QL factorization, and converges roughly as [(λi−σi)∕(λi+1−σi)]j where λi is the eigenvalue located in the (i,i) position of the final diagonal matrix and the eigenvalues are ordered (|λ1|<|λ2|<…<|λn|), and j is the iteration. The “educated guess” σi for the eigenvalue λi is obtained from the analytic determination of the eigenvalues of (k×k)-submatrices of A̲, in the vicinity of the ith element, where k=0,1,2,3 (here, k=0 means that the implicit shift vanishes, σi=0). The routine optionally calculates the rotation matrix Z̲, such that Z̲−1A̲Z̲=D̲ where A̲ is the input matrix and D̲ is the diagonal matrix containing the eigenvalues. The ith column of Z̲ then is the eigenvector of A̲ corresponding to the eigenvalue found at the element D̲(i,i), in the ith position on the diagonal of the matrix D̲. Unusual features: For simplicity, the “wrapper” program which contains an example application and the HTDQLS routine are distributed in the same file.

AB - We describe a matrix diagonalization algorithm for complex symmetric (not Hermitian) matrices, A̲=A̲T, which is based on a two-step algorithm involving generalized Householder reflections based on the indefinite inner product 〈u̲,v̲〉∗=∑iuivi. This inner product is linear in both arguments and avoids complex conjugation. The complex symmetric input matrix is transformed to tridiagonal form using generalized Householder transformations (first step). An iterative, generalized QL decomposition of the tridiagonal matrix employing an implicit shift converges toward diagonal form (second step). The QL algorithm employs iterative deflation techniques when a machine-precision zero is encountered “prematurely” on the super-/sub-diagonal. The algorithm allows for a reliable and computationally efficient computation of resonance and antiresonance energies which emerge from complex-scaled Hamiltonians, and for the numerical determination of the real energy eigenvalues of pseudo-Hermitian and PT-symmetric Hamilton matrices. Numerical reference values are provided. Program summary Program Title: HTDQLS Program Files doi: http://dx.doi.org/10.17632/x24wjxtrsg.1 Licensing provisions: GPLv3 Programming language: Fortran 90 using fixed form notation Nature of problem: Calculating the eigenvalues and optionally the eigenvectors of complex symmetric (non-Hermitian), densely populated matrices. Solution method: The complex symmetric (not Hermitian) input matrix is diagonalized in two steps. First step: The matrix is tridiagonalized via a series of (n−2) generalized Householder reflections, where n is the rank of the input matrix. Second step: The tridiagonal matrix is diagonalized via a generalization of the “chasing the bulge” technique, which is an iterative process utilizing an implicitly shifted initial rotation followed by (n−2) Givens rotations. This technique is an implementation of QL factorization, and converges roughly as [(λi−σi)∕(λi+1−σi)]j where λi is the eigenvalue located in the (i,i) position of the final diagonal matrix and the eigenvalues are ordered (|λ1|<|λ2|<…<|λn|), and j is the iteration. The “educated guess” σi for the eigenvalue λi is obtained from the analytic determination of the eigenvalues of (k×k)-submatrices of A̲, in the vicinity of the ith element, where k=0,1,2,3 (here, k=0 means that the implicit shift vanishes, σi=0). The routine optionally calculates the rotation matrix Z̲, such that Z̲−1A̲Z̲=D̲ where A̲ is the input matrix and D̲ is the diagonal matrix containing the eigenvalues. The ith column of Z̲ then is the eigenvector of A̲ corresponding to the eigenvalue found at the element D̲(i,i), in the ith position on the diagonal of the matrix D̲. Unusual features: For simplicity, the “wrapper” program which contains an example application and the HTDQLS routine are distributed in the same file.

KW - Complex symmetric matrix diagonalization

KW - Deflation techniques

KW - Implicit shift

KW - Indefinite inner product

UR - http://www.scopus.com/inward/record.url?scp=85028728933&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85028728933&partnerID=8YFLogxK

U2 - 10.1016/j.cpc.2017.06.014

DO - 10.1016/j.cpc.2017.06.014

M3 - Article

AN - SCOPUS:85028728933

VL - 221

SP - 304

EP - 316

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

ER -