Solving large-scale nonlinear eigenvalue problems by rational interpolation and resolvent sampling based Rayleigh–Ritz method

Jinyou Xiao, Chuanzeng Zhang, Tsung-Min Hwang, Tetsuya Sakurai

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by the matrix structures, properties of the eigen-solutions, sizes of the problems, etc. This paper aims to remove those limitations and develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach, the RIA provides the possibility to select sampling points in more general regions and has advantages in improving the accuracy and reducing the computational cost. Second, a resolvent sampling scheme using the RIA is proposed to construct reliable search spaces for the Rayleigh–Ritz procedure, based on which a robust eigen-solver, called resolvent sampling based Rayleigh–Ritz method (RSRR), is developed for solving general NEPs. The RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of the RSRR are demonstrated by a variety of benchmark and application examples.

Original languageEnglish
Pages (from-to)776-800
Number of pages25
JournalInternational Journal for Numerical Methods in Engineering
Volume110
Issue number8
DOIs
Publication statusPublished - 2017 May 25

Fingerprint

Rayleigh-Ritz Method
Rational Interpolation
Nonlinear Eigenvalue Problem
Resolvent
Interpolation
Sampling
Contour integral
Computational Science
Matrix Function
Engineering Application
Search Space
Computational Cost
Analytic function
Numerical Solution
Benchmark
Engineering
Theorem
Costs

Keywords

  • Rayleigh–Ritz procedure
  • boundary element methods
  • eigenvalue problems
  • finite element methods
  • nonlinear solvers

ASJC Scopus subject areas

  • Numerical Analysis
  • Engineering(all)
  • Applied Mathematics

Cite this

Solving large-scale nonlinear eigenvalue problems by rational interpolation and resolvent sampling based Rayleigh–Ritz method. / Xiao, Jinyou; Zhang, Chuanzeng; Hwang, Tsung-Min; Sakurai, Tetsuya.

In: International Journal for Numerical Methods in Engineering, Vol. 110, No. 8, 25.05.2017, p. 776-800.

Research output: Contribution to journalArticle

@article{a4dec2638ff947f298cb2cca2ebbd6b5,
title = "Solving large-scale nonlinear eigenvalue problems by rational interpolation and resolvent sampling based Rayleigh–Ritz method",
abstract = "Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by the matrix structures, properties of the eigen-solutions, sizes of the problems, etc. This paper aims to remove those limitations and develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach, the RIA provides the possibility to select sampling points in more general regions and has advantages in improving the accuracy and reducing the computational cost. Second, a resolvent sampling scheme using the RIA is proposed to construct reliable search spaces for the Rayleigh–Ritz procedure, based on which a robust eigen-solver, called resolvent sampling based Rayleigh–Ritz method (RSRR), is developed for solving general NEPs. The RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of the RSRR are demonstrated by a variety of benchmark and application examples.",
keywords = "Rayleigh–Ritz procedure, boundary element methods, eigenvalue problems, finite element methods, nonlinear solvers",
author = "Jinyou Xiao and Chuanzeng Zhang and Tsung-Min Hwang and Tetsuya Sakurai",
year = "2017",
month = "5",
day = "25",
doi = "10.1002/nme.5441",
language = "English",
volume = "110",
pages = "776--800",
journal = "International Journal for Numerical Methods in Engineering",
issn = "0029-5981",
publisher = "John Wiley and Sons Ltd",
number = "8",

}

TY - JOUR

T1 - Solving large-scale nonlinear eigenvalue problems by rational interpolation and resolvent sampling based Rayleigh–Ritz method

AU - Xiao, Jinyou

AU - Zhang, Chuanzeng

AU - Hwang, Tsung-Min

AU - Sakurai, Tetsuya

PY - 2017/5/25

Y1 - 2017/5/25

N2 - Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by the matrix structures, properties of the eigen-solutions, sizes of the problems, etc. This paper aims to remove those limitations and develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach, the RIA provides the possibility to select sampling points in more general regions and has advantages in improving the accuracy and reducing the computational cost. Second, a resolvent sampling scheme using the RIA is proposed to construct reliable search spaces for the Rayleigh–Ritz procedure, based on which a robust eigen-solver, called resolvent sampling based Rayleigh–Ritz method (RSRR), is developed for solving general NEPs. The RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of the RSRR are demonstrated by a variety of benchmark and application examples.

AB - Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by the matrix structures, properties of the eigen-solutions, sizes of the problems, etc. This paper aims to remove those limitations and develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach, the RIA provides the possibility to select sampling points in more general regions and has advantages in improving the accuracy and reducing the computational cost. Second, a resolvent sampling scheme using the RIA is proposed to construct reliable search spaces for the Rayleigh–Ritz procedure, based on which a robust eigen-solver, called resolvent sampling based Rayleigh–Ritz method (RSRR), is developed for solving general NEPs. The RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of the RSRR are demonstrated by a variety of benchmark and application examples.

KW - Rayleigh–Ritz procedure

KW - boundary element methods

KW - eigenvalue problems

KW - finite element methods

KW - nonlinear solvers

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

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

U2 - 10.1002/nme.5441

DO - 10.1002/nme.5441

M3 - Article

AN - SCOPUS:85006086299

VL - 110

SP - 776

EP - 800

JO - International Journal for Numerical Methods in Engineering

JF - International Journal for Numerical Methods in Engineering

SN - 0029-5981

IS - 8

ER -