Parallel matrix function evaluation via initial value ODE modeling

DOI10.1016/J.CAMWA.2016.04.036MaRDI QIDQ2007305zbMATH OpenOpenAlexFDO

Authors Jean-Paul Chehab, Mădălina Petcu

Publication date 12 October 2020

Published in Computers & Mathematics with Applications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1505.07959

zbMATH Keywords

ODEsfunctions of matricesparareal methods

Mathematics Subject Classification ID

Parallel numerical computation (65Y05) Stability of solutions to ordinary differential equations (34D20)

Abstract: The purpose of this article is to propose ODE based approaches for the numerical evaluation of matrix functions $f(A)$, a question of major interest in the numerical linear algebra. To this end, we model $f(A)$ as the solution at a finite time $T$ of a time dependent equation. We use parallel algorithms, such as the parareal method, on the time interval $[0,T]$ in order to solve the evolution equation obtained. When $f(A)$ is reached as a stable steady state, it can be computed by combining parareal algorithms and optimal control techniques. Numerical illustrations are given.

Recommendations

• Parallel Factorizations in Numerical Analysis
• Parallel solution in time of ODEs: Some achievements and perspectives
• Computing A^\alpha, \log(A), and Related Matrix Functions by Contour Integrals
• Computing matrix functions
• Sparse approximations of matrix functions via numerical integration of ODEs

Cites work

• A ``parareal in time discretization of PDE's
• A new inversion-free method for a rational matrix equation
• A parallel block cyclic reduction algorithm for the fast solution of elliptic equations
• Analysis of a Krylov subspace enhanced parareal algorithm for linear problems
• Analysis of a Modified Parareal Algorithm for Second-Order Ordinary Differential Equations
• Analysis of the Parareal Time‐Parallel Time‐Integration Method
• Chebyshev rational approximations to \(e^{-x}\) in \([0,+\infty)\) and applications to heat-conduction problems
• Control Perspectives on Numerical Algorithms and Matrix Problems
• Convergence properties of some block Krylov subspace methods for multiple linear systems
• Decay bounds and \(O(n)\) algorithms for approximating functions of sparse matrices
• Differential equations and solution of linear systems
• Efficient algorithms for the matrix cosine and sine
• Efficient and stable Arnoldi restarts for matrix functions based on quadrature
• Functions of Matrices
• Inexact and truncated parareal-in-time Krylov subspace methods for parabolic optimal control problems
• Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later
• PARAEXP: a parallel integrator for linear initial-value problems
• Rational Krylov approximation of matrix functions: numerical methods and optimal pole selection
• Residual methods for the large-scale matrix \(p\)th root and some related problems
• Time-parallel implicit integrators for the near-real-time prediction of linear structural dynamic responses
• Unconditionally stable scheme for Riccati equation

Cited in

(4)

• Sparse approximations of matrix functions via numerical integration of ODEs
• A parallel spectral deferred correction method for first-order evolution problems
• A parareal algorithm for a highly oscillating Vlasov-Poisson system with reduced models for the coarse solving
• Parallel computation of functions of matrices and their action on vectors for exponential integrators

Describes a project that uses

Uses Software

• mftoolbox
• PARAEXP

This page was built for publication: Parallel matrix function evaluation via initial value ODE modeling

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2007305)