Computing a<sup>α</sup>, log(A), and related matrix functions by contour integrals

Author: 

Hale, N
Higham, NJ
Trefethen, LN

Journal: 

SIAM Journal on Numerical Analysis

Publication Date: 

10 November 2008

Last Updated: 

2018-12-06T19:06:58.533+00:00

Issue: 

5

DOI: 

10.1137/070700607

Volume: 

46

page: 

2505-2523

abstract: 

New methods are proposed for the numerical evaluation of f(A) or f(A)b, where f(A) is a function such as A1/2or log(A) with singularities in (-∞, 0] and A is a matrix with eigenvalues on or near (0, ∞). The methods are based on combining contour integrals evaluated by the periodic trapezoid rule with conformal maps involving Jacobi elliptic functions. The convergence is geometric, so that the computation of f(A)fe is typically reduced to one or two dozen linear system solves, which can be carried out in parallel. © 2008 Society for Industrial and Applied Mathematics.

Symplectic id: 

188993

Submitted to ORA: 

Not Submitted

Publication Type: 

Journal Article