Complex moment-based methods for differential eigenvalue problems

05/02/2022
by   Akira Imakura, et al.
0

This paper considers computing partial eigenpairs of differential eigenvalue problems (DEPs) such that eigenvalues are in a certain region on the complex plane. Recently, based on a "solve-then-discretize" paradigm, an operator analogue of the FEAST method has been proposed for DEPs without discretization of the coefficient operators. Compared to conventional "discretize-then-solve" approaches that discretize the operators and solve the resulting matrix problem, the operator analogue of FEAST exhibits much higher accuracy; however, it involves solving a large number of ordinary differential equations (ODEs). In this paper, to reduce the computational costs, we propose operation analogues of Sakurai-Sugiura-type complex moment-based eigensolvers for DEPs using higher-order complex moments and analyze the error bound of the proposed methods. We show that the number of ODEs to be solved can be reduced by a factor of the degree of complex moments without degrading accuracy, which is verified by numerical results. Numerical results demonstrate that the proposed methods are over five times faster compared with the operator analogue of FEAST for several DEPs while maintaining almost the same high accuracy. This study is expected to promote the "solve-then-discretize" paradigm for solving DEPs and contribute to faster and more accurate solutions in real-world applications.

READ FULL TEXT

page 20

page 21

research
05/07/2021

A semigroup method for high dimensional elliptic PDEs and eigenvalue problems based on neural networks

In this paper, we propose a semigroup method for solving high-dimensiona...
research
11/27/2019

A Scalable Framework for Solving Fractional Diffusion Equations

The study of fractional order differential operators is receiving renewe...
research
04/17/2023

A novel higher-order numerical method for parabolic integro-fractional differential equations based on wavelets and L2-1_σ scheme

This paper aims to construct an efficient and highly accurate numerical ...
research
10/27/2019

A Jacobi spectral method for computing eigenvalue gaps and their distribution statistics of the fractional Schrödinger operator

We propose a spectral method by using the Jacobi functions for computing...
research
07/02/2020

Projection method for partial eigenproblems of linear matrix pencils

In this contribution, we extend a projection method for square matrices ...
research
09/15/2021

Hybrid quadrature moment method for accurate and stable representation of non-Gaussian processes and their dynamics

Solving the population balance equation (PBE) for the dynamics of a disp...

Please sign up or login with your details

Forgot password? Click here to reset