Abstract:The computation of off-diagonal blocks of matrix functions $f(T)$, where $T$ is block triangular, poses a challenging problem in scientific computing. We present a novel algorithm that exploits the structure of block triangular matrices, generalizing the algorithm of Al-Mohy and Higham for computing the Fréchet derivative of the matrix exponential. This work has significant applications in fields such as exponential integrators for solving systems of first-order differential equations, Hamiltonian linear systems in control theory, and option pricing in finance. Our approach introduces a linear operator that maps off-diagonal blocks of $T$ into their counterparts in $f(T)$. By studying the algebraic properties of the operator, we establish a comprehensive computational framework, paving the way to extend existing Fréchet derivative algorithms of matrix functions to more general settings. For the matrix exponential, in particular, the algorithm employs the scaling and squaring method with diagonal Padé approximants to $\exp(x)$, with parameters chosen based on a rigorous backward error analysis, which notably does not depend on the norm of the off-diagonal blocks. The numerical experiment demonstrates that our algorithm surpasses existing algorithms in terms of accuracy and efficiency, making it highly valuable for a wide range of applications.

Fast Hankel tensor–vector product and its application to exponential data fitting

Fast Tensor-Vector Products

A Nonuniform Fast Hankel Transform

Fast Hankel Transform And Its Application To Beam

A Lanczos Bidiagonalization Algorithm for Hankel Matrices

A Split Fast Fourier Transform Algorithm for Block Toeplitz Matrix-Vector Multiplication

Efficient algorithms for Tucker decomposition via approximate matrix multiplication

On the finite section method for computing exponentials of doubly-infinite skew-Hermitian matrices

A fast algorithm for solving circulant tensor systems

Generalizing the Fréchet Derivative Algorithm for the Matrix Exponential

A Fast Algorithm for Cosine Transform Based Tensor Singular Value Decomposition

Fast Fourier Solvers for the Tensor Product High-Order FEM for a Poisson Type Equation

High Performance Hardware Architecture for Singular Spectrum Analysis of Hankel Tensors.

A fast algorithm to factorize high-dimensional tensor product matrices used in genetic models

Numerical Methods for Fast Nonlinear Fourier Transformation, Part I: Exponential Runge-Kutta and Linear Multistep Methods

Sum-of-Squares Decomposition and Augmented Vandermonde Decomposition of Strong Hankel Tensors: Theory and Algorithms

Tucker Tensor Approach for Accelerating Fock Exchange Computations in a Real-Space Finite-Element Discretization of Generalized Kohn-Sham Density Functional Theory

Fast Method to Compute Tensor Product 2-D Wavelet Transforms

Efficient scaling and squaring method for the matrix exponential

Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction Via Low-Rank Hankel Matrix Completion

A Fast Matrices Decomposition Algorithm for Two-Dimensional Elliptic Problem Using Legendre Spectral Collocation Approximation