Premium
Domain decomposition approaches for accelerating contour integration eigenvalue solvers for symmetric eigenvalue problems
Author(s) -
Kalantzis Vassilis,
Kestyn James,
Polizzi Eric,
Saad Yousef
Publication year - 2018
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.2154
Subject(s) - eigenvalues and eigenvectors , resolvent , domain decomposition methods , mathematics , methods of contour integration , domain (mathematical analysis) , divide and conquer eigenvalue algorithm , operator (biology) , scheme (mathematics) , algorithm , eigendecomposition of a matrix , matrix (chemical analysis) , mathematical optimization , computer science , mathematical analysis , finite element method , gene , transcription factor , thermodynamics , biochemistry , physics , chemistry , repressor , quantum mechanics , materials science , composite material
Summary This paper discusses techniques for computing a few selected eigenvalue–eigenvector pairs of large and sparse symmetric matrices. A recently developed class of techniques to solve this type of problems is based on integrating the matrix resolvent operator along a complex contour that encloses the interval containing the eigenvalues of interest. This paper considers such contour integration techniques from a domain decomposition viewpoint and proposes two schemes. The first scheme can be seen as an extension of domain decomposition linear system solvers in the framework of contour integration methods for eigenvalue problems, such as FEAST. The second scheme focuses on integrating the resolvent operator primarily along the interface region defined by adjacent subdomains. A parallel implementation of the proposed schemes is described, and results on distributed computing environments are reported. These results show that domain decomposition approaches can lead to reduced run times and improved scalability.