A framework of the harmonic Arnoldi method for evaluating -functions with applications to exponential integrators
From MaRDI portal
Publication:295370
DOI10.1007/S10444-015-9433-0zbMATH Open1338.65123arXiv1409.0349OpenAlexW2202953938MaRDI QIDQ295370FDOQ295370
Authors: Gang Wu, Lu Zhang, Tingting Xu 
Publication date: 13 June 2016
Published in: Advances in Computational Mathematics (Search for Journal in Brave)
Abstract: In recent years, a great deal of attention has been focused on numerically solving exponential integrators. The important ingredient to the implementation of exponential integrators is the efficient and accurate evaluation of the so called -functions on a given vector. The Krylov subspace method is an important technique for this problem. For this type of method, however, restarts become essential for the sake of storage requirements or due to the growing computational complexity of evaluating the matrix function on a Hessenberg matrix of growing size. Another problem in computing -functions is the lack of a clear residual notion. The contribution of this work is threefold. First, we introduce a framework of the harmonic Arnoldi method for -functions, which is based on the residual and the oblique projection technique. Second, we establish the relationship between the harmonic Arnoldi approximation and the classical Arnoldi approximation, and compare the harmonic Arnoldi method with the Arnoldi method from a theoretical point of view. Third, we apply the thick-restarting strategy to the harmonic Arnoldi method, and propose a thick-restated harmonic Arnoldi algorithm for evaluating -functions. An advantage of the new algorithm is that we can compute several -functions simultaneously in the same search subspace. We show the merit of augmenting approximate eigenvectors in the search subspace, and give insight into the relationship between the error and the residual of -functions. Numerical experiments show the superiority of our new algorithm over many state-of-the-art algorithms for the computation of -functions.
Full work available at URL: https://arxiv.org/abs/1409.0349
Recommendations
- Efficient and accurate computation for the \(\varphi\)-functions arising from exponential integrators
- A moment-matching Arnoldi iteration for linear combinations of \(\phi\) functions
- Algorithm 919: A Krylov subspace algorithm for evaluating the \(\phi\)-functions appearing in exponential integrators
- Residual and restarting in Krylov subspace evaluation of the \(\varphi\) function
- Block Krylov subspace methods for approximating the linear combination of \(\varphi\)-functions arising in exponential integrators
exponential integratorsmatrix exponential\(\varphi\)-functionsharmonic Arnoldi methodoblique projection methodthick-restarting strategy
Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Iterative numerical methods for linear systems (65F10) Numerical computation of matrix exponential and similar matrix functions (65F60)
Cites Work
- Expokit
- Algorithm 919: A Krylov subspace algorithm for evaluating the \(\phi\)-functions appearing in exponential integrators
- Title not available (Why is that?)
- The scaling and squaring method for the matrix exponential revisited
- Functions of Matrices
- Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later
- Matrix computations.
- Fractional differential equations. An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications
- Finite difference approximations for fractional advection-dispersion flow equations
- A direct \(O(N \log ^{2} N)\) finite difference method for fractional diffusion equations
- Matrix algorithms. Vol. 2: Eigensystems
- Multigrid method for fractional diffusion equations
- Option pricing when underlying stock returns are discontinuous
- Generalized integrating factor methods for stiff PDEs
- Exponential integrators
- Shift-invert Arnoldi approximation to the Toeplitz matrix exponential
- Computing the action of the matrix exponential, with an application to exponential integrators
- Exponential Integrators for Large Systems of Differential Equations
- GMRES with Deflated Restarting
- New adaptive exponential propagation iterative methods of Runge-Kutta type
- On Krylov Subspace Approximations to the Matrix Exponential Operator
- The scaling and modified squaring method for matrix functions related to the exponential
- Algorithmic challenges in computational molecular biophysics
- A Krylov projection method for systems of ODEs
- Resolvent Krylov subspace approximation to operator functions
- RD-rational approximations of the matrix exponential
- A residual based error estimate for Leja interpolation of matrix functions
- A harmonic restarted Arnoldi algorithm for calculating eigenvalues and determining multiplicity
- Comparison of software for computing the action of the matrix exponential
- The Gautschi time stepping scheme for edge finite element discretizations of the Maxwell equations
- Thick-restart Lanczos method for large symmetric eigenvalue problems
- Residual, restarting, and Richardson iteration for the matrix exponential
- Convergence analysis of an extended Krylov subspace method for the approximation of operator functions in exponential integrators
- Rational Krylov approximation of matrix functions: numerical methods and optimal pole selection
- Efficient and stable Arnoldi restarts for matrix functions based on quadrature
- Error Estimates and Evaluation of Matrix Functions via the Faber Transform
- A restarted Lanczos approximation to functions of a symmetric matrix
- On analysis of projection methods for rational function approximation to the matrix exponential
- A new investigation of the extended Krylov subspace method for matrix function evaluations
- Deflated restarting for matrix functions
- Using the restricted-denominator rational Arnoldi method for exponential integrators
- On Restart and Error Estimation for Krylov Approximation of $w=f(A)v$
- A black-box rational Arnoldi variant for Cauchy-Stieltjes matrix functions
- Analysis of Some Krylov Subspace Approximations to the Matrix Exponential Operator
- Efficient Solution of Parabolic Equations by Krylov Approximation Methods
- Extended Krylov Subspaces: Approximation of the Matrix Square Root and Related Functions
- Using Nonorthogonal Lanczos Vectors in the Computation of Matrix Functions
- Krylov subspace approximations for the exponential Euler method: error estimates and the harmonic Ritz approximant
- Approximate solutions and eigenvalue bounds from Krylov subspaces
- The restarted shift-and-invert Krylov method for matrix functions.
- Regularization of nonlinear ill-posed problems by exponential integrators
- Uniform approximation of \(\varphi\)-functions in exponential integrators by a rational Krylov subspace method with simple poles
- Convergence of Restarted Krylov Subspace Methods for Stieltjes Functions of Matrices
- A Restarted Krylov Subspace Method for the Evaluation of Matrix Functions
- Preconditioning Lanczos Approximations to the Matrix Exponential
- Rational Krylov methods in exponential integrators for European option pricing.
- A thick-restarted block Arnoldi algorithm with modified Ritz vectors for large eigenproblems
- A generalization of the steepest descent method for matrix functions
- Evaluating matrix functions for exponential integrators via Carathéodory-Fejér approximation and contour integrals
- Implementation of a restarted Krylov subspace method for the evaluation of matrix functions
Cited In (7)
- Computing a matrix function for exponential integrators.
- A moment-matching Arnoldi iteration for linear combinations of \(\phi\) functions
- A study of defect-based error estimates for the Krylov approximation of \(\varphi\)-functions
- Efficient and accurate computation for the \(\varphi\)-functions arising from exponential integrators
- A \(\mu\)-mode approach for exponential integrators: actions of \(\varphi\)-functions of Kronecker sums
- Using the restricted-denominator rational Arnoldi method for exponential integrators
- A shifted block FOM algorithm with deflated restarting for matrix exponential computations
Uses Software
This page was built for publication: A framework of the harmonic Arnoldi method for evaluating \(\varphi\)-functions with applications to exponential integrators
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q295370)
