The global extended-rational Arnoldi method for matrix function approximation

03/31/2020
by   A. H. Bentbib, et al.
0

The numerical computation of matrix functions such as f(A)V, where A is an n× n large and sparse square matrix, V is an n × p block with p≪ n and f is a nonlinear matrix function, arises in various applications such as network analysis (f(t)=exp(t) or f(t)=t^3), machine learning (f(t)=log(t)), theory of quantum chromodynamics (f(t)=t^1/2), electronic structure computation, and others. In this work, we propose the use of global extended-rational Arnoldi method for computing approximations of such expressions. The derived method projects the initial problem onto an global extended-rational Krylov subspace RK^e_m(A,V)=span({∏_i=1^m(A-s_iI_n)^-1V,...,(A-s_1I_n)^-1V,V,AV, ...,A^m-1V}) of a low dimension. An adaptive procedure for the selection of shift parameters {s_1,...,s_m} is given. The proposed method is also applied to solve parameter dependent systems. Numerical examples are presented to show the performance of the global extended-rational Arnoldi for these problems.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
08/28/2023

A Rational Krylov Subspace Method for the Computation of the Matrix Exponential Operator

The computation of approximating e^tA B, where A is a large sparse matri...
research
01/31/2017

Interior Eigensolver for Sparse Hermitian Definite Matrices Based on Zolotarev's Functions

This paper proposes an efficient method for computing selected generaliz...
research
12/24/2022

Speeding up Krylov subspace methods for computing f(A)b via randomization

This work is concerned with the computation of the action of a matrix fu...
research
06/01/2022

A barycentric trigonometric Hermite interpolant via an iterative approach

In this paper an interative approach for constructing the Hermite interp...
research
11/17/2020

Solving Symmetric and Positive Definite Second-Order Cone Linear Complementarity Problem by A Rational Krylov Subspace Method

The second-order cone linear complementarity problem (SOCLCP) is a gener...
research
02/28/2020

Generalized Rational Variable Projection With Application in ECG Compression

In this paper we develop an adaptive transform-domain technique based on...
research
08/06/2019

Rational Krylov for Stieltjes matrix functions: convergence and pole selection

Evaluating the action of a matrix function on a vector, that is x=f( M)v...

Please sign up or login with your details

Forgot password? Click here to reset