Practical shift choice in the shift-and-invert Krylov subspace evaluations of the matrix exponential

09/28/2019
by   Alexandr Katrutsa, et al.
0

We propose two methods to find a proper shift parameter in the shift-and-invert method for computing matrix exponential matrix-vector products. These methods are useful in the case of matrix exponential action has to be computed for a number of vectors. The first method is based on the zero-order optimization of the mean residual norm for a given number of initial vectors. The second method processes the vectors one-by-one and estimates, for each vector, the derivative of the residual norm as a function of the shift parameter. The estimated derivative value is then used to update the shift value for the next vector. To demonstrate the performance of the proposed methods we perform numerical experiments for two-dimensional non-stationary convection-diffusion equation with discontinuous coefficients and two-dimensional anisotropic diffusion equation.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/05/2019

An accurate restarting for shift-and-invert Krylov subspaces computing matrix exponential actions of nonsymmetric matrices

An accurate residual–time (AccuRT) restarting for computing matrix expon...
research
12/25/2018

ART: adaptive residual--time restarting for Krylov subspace matrix exponential evaluations

In this paper a new restarting method for Krylov subspace matrix exponen...
research
07/15/2021

Coarse grid corrections in Krylov subspace evaluations of the matrix exponential

A coarse grid correction (CGC) approach is proposed to enhance the effic...
research
12/22/2017

Monotone Difference Schemes for Convection-Dominated Diffusion-Reaction Equations Based on Quadratic Spline

A three-point monotone difference scheme is proposed for solving a one-d...
research
09/14/2022

Unprojected recycled block Krylov subspace methods for shifted systems

The use of block Krylov subspace methods for computing the solution to a...
research
03/09/2020

Vector logic and counterfactuals

In this work we investigate the representation of counterfactual conditi...
research
10/25/2021

Operator Augmentation for Model-based Policy Evaluation

In model-based reinforcement learning, the transition matrix and reward ...

Please sign up or login with your details

Forgot password? Click here to reset