Abstract:This paper presents a new algorithm KIOPS for computing linear combinations of $\varphi$-functions that appear in exponential integrators. This algorithm is suitable for large-scale problems in computational physics where little or no information about the spectrum or norm of the Jacobian matrix is known a priori. We first show that such problems can be solved efficiently by computing a single exponential of a modified matrix. Then our approach is to compute an appropriate basis for the Krylov subspace using the incomplete orthogonalization procedure and project the matrix exponential on this subspace. We also present a novel adaptive procedure that significantly reduces the computational complexity of exponential integrators. Our numerical experiments demonstrate that KIOPS outperforms the current state-of-the-art adaptive Krylov algorithm phipm.

What problem does this paper attempt to address?

### Problems the paper attempts to solve The paper aims to solve the problems of efficiency and accuracy in the large - scale rigid matrix exponential integration calculation. Specifically, the paper proposes a new algorithm KIOPS (Krylov with Incomplete Orthogonalization Procedure Solver) for efficiently calculating the linear combination of ϕ - functions that appear in exponential integrators. #### Main problem background When solving large - scale ordinary differential equations (ODEs), exponential integration methods have attracted attention due to their high accuracy and good stability. However, the main computational cost of these methods lies in calculating the product of the exponential function of the Jacobian matrix or its approximation and a vector. For large - scale problems, especially when little is known about the spectral or norm information of the Jacobian matrix, existing methods often face high computational costs or numerical stability problems. #### Core improvements of the KIOPS algorithm 1. **Single - exponential calculation**: The paper first shows how to efficiently solve these problems by calculating the exponential of a modified matrix. 2. **Incomplete orthogonalization**: Use the Incomplete Orthogonalization Procedure (IOP) to construct the Krylov subspace basis and project the matrix exponential onto this subspace. 3. **Adaptive strategy**: Propose a new adaptive procedure, which significantly reduces the computational complexity of the exponential integrator. #### Advantages and innovation points - **Wide range of applicability**: KIOPS is suitable for large - scale physical computing problems, especially when the spectral or norm information of the Jacobian matrix is not known in advance. - **High efficiency**: By reducing the number of explicit matrix - vector multiplications, KIOPS avoids rounding errors that become more sensitive as the order increases. - **Adaptive optimization**: KIOPS further improves the computational efficiency by dynamically adjusting the sub - step size and the size of the Krylov subspace. #### Experimental verification The paper verifies the performance of KIOPS through numerical experiments. The results show that when dealing with large - scale rigid matrix problems, KIOPS has higher efficiency and accuracy compared to the current state - of - the - art adaptive Krylov algorithm phipm. ### Summary By introducing the KIOPS algorithm, this paper solves the key challenges in the large - scale rigid matrix exponential integration calculation and provides a more efficient and stable solution, especially suitable for large - scale simulation problems in fields such as atmospheric fluid mechanics.

KIOPS: A fast adaptive Krylov subspace solver for exponential integrators

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Low-synchronization Arnoldi Methods for the Matrix Exponential with Application to Exponential Integrators

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

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

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

A Hybrid Direct-Iterative Method for Solving KKT Linear Systems

EPIRK-W and EPIRK-K time discretization methods

Parallel-in-time high-order multiderivative IMEX solvers

A Flexible Krylov Solver for Shifted Systems with Application to Oscillatory Hydraulic Tomography

GPU-based Parallel Solver Via the Kantorovich Theorem for the Nonlinear Bernstein Polynomial Systems.

Parallel-in-time Multilevel Krylov Methods: A Prototype

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

A robust and time-parallel preconditioner for parabolic reconstruction problems using Isogeometric Analysis

Quantum Krylov-Subspace Method Based Linear Solver

A preconditioned Krylov subspace method for linear inverse problems with general-form Tikhonov regularization

A splitting-based KPIK method for eddy current optimal control problems in an all-at-once approach

Libkrylov: A modular open‐source software library for extremely large on‐the‐fly matrix computations

A new approach to constructing efficient stiffly accurate exponential propagation iterative methods of Runge-Kutta type (EPIRK)

Krylov Solvers for Interior Point Methods with Applications in Radiation Therapy and Support Vector Machines