Abstract:We present a method for computing actions of the exponential-like $\varphi$-functions for a Kronecker sum $K$ of $d$ arbitrary matrices $A_\mu$. It is based on the approximation of the integral representation of the $\varphi$-functions by Gaussian quadrature formulas combined with a scaling and squaring technique. The resulting algorithm, which we call PHIKS, evaluates the required actions by means of $\mu$-mode products involving exponentials of the small sized matrices $A_\mu$, without forming the large sized matrix $K$ itself. PHIKS, which profits from the highly efficient level 3 BLAS, is designed to compute different $\varphi$-functions applied on the same vector or a linear combination of actions of $\varphi$-functions applied on different vectors. In addition, thanks to the underlying scaling and squaring techniques, the desired quantities are available simultaneously at suitable time scales. All these features allow the effective usage of PHIKS in the exponential integration context. In fact, our newly designed method has been tested on popular exponential Runge--Kutta integrators of stiff order from one to four, in comparison with state-of-the-art algorithms for computing actions of $\varphi$-functions. The numerical experiments with discretized semilinear evolutionary 2D or 3D advection--diffusion--reaction, Allen--Cahn, and Brusselator equations show the superiority of the proposed $\mu$-mode approach.

LeXInt: Package for Exponential Integrators employing Leja interpolation

LeXInt: GPU-accelerated Exponential Integrators package

Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE

Exponential Integrators for Phase-Field Equations using Pseudo-spectral Methods: A Python Implementation

EPIRK-W and EPIRK-K time discretization methods

Full- and low-rank exponential Euler integrators for the Lindblad equation

A $μ$-mode approach for exponential integrators: actions of $φ$-functions of Kronecker sums

Exact Integration for singular Zienkiewicz and Guzman-Neilan Finite Elements with Implementation

differint: A Python Package for Numerical Fractional Calculus

A discrete framework for the interpolation of Banach spaces

ZERNIPAX: A Fast and Accurate Zernike Polynomial Calculator in Python

A -mode approach for exponential integrators: actions of -functions of Kronecker sums

Parallel-in-time high-order multiderivative IMEX solvers

High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation

KIOPS: A fast adaptive Krylov subspace solver for exponential integrators

Geometric Exponential Integrators

A novel high-order explicit exponential integrator scheme for the space–time fractional Bloch–Torrey equation

Numerical methods for computing the discrete and continuous Laplace transforms

Trigonometric time integrators for the Zakharov system

Symmetric and Arbitrarily High-Order Birkhoff-Hermite Time Integrators and Their Long-Time Behaviour for Solving Nonlinear Klein-Gordon Equations

Third order, uniform in low to high oscillatory coefficients, exponential integrators for Klein-Gordon equations