Abstract:The numerical integration of stiff equations is a challenging problem that needs to be approached by specialized numerical methods. Exponential integrators form a popular class of such methods since they are provably robust to stiffness and have been successfully applied to a variety of problems. The dynamical low- \rank approximation is a recent technique for solving high-dimensional differential equations by means of low-rank approximations. However, the domain is lacking numerical methods for stiff equations since existing methods are either not robust-to-stiffness or have unreasonably large hidden constants. In this paper, we focus on solving large-scale stiff matrix differential equations with a Sylvester-like structure, that admit good low-rank approximations. We propose two new methods that have good convergence properties, small memory footprint and that are fast to compute. The theoretical analysis shows that the new methods have order one and two, respectively. We also propose a practical implementation based on Krylov techniques. The approximation error is analyzed, leading to a priori error bounds and, therefore, a mean for choosing the size of the Krylov space. Numerical experiments are performed on several examples, confirming the theory and showing good speedup in comparison to existing techniques.

Large-Scale Stein and Lyapunov Equations 3 2 Large-Scale Stein Equations

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

Scalable Distributed Least Square Algorithms for Large-Scale Linear Equations Via an Optimization Approach

Low-Rank Methods for Solving Discrete-Time Projected Lyapunov Equations

A new projection method for solving large Sylvester equations

Linear Scaling Algorithms for Solving High-Dimensional Nonlinear Parabolic Differential Equations

An Algebraic Substructuring Method for Large-Scale Eigenvalue Calculation

Decoupled Structure-Preserving Doubling Algorithm with Truncation for Large-Scale Algebraic Riccati Equations

A low-rank algorithm for solving Lyapunov operator $φ$-functions within the matrix-valued exponential integrators

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

Stein's Method for Law of Large Numbers under Sublinear Expectations

Sharpness on the Lower Bound of the Lifespan of Solutions to Nonlinear Wave Equations

Symmetric and arbitrarily high-order time-stepping methods for solving nonlinear Klein–Gordon equations

The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

The Nonlinear Steepest Descent Method to Long-Time Asymptotics of the Coupled Nonlinear Schrödinger Equation

Soliton Solutions For Some Equations In The (1+2)-Dimensional Hyperbolic Su(N) Akns System

Large Time Behavior of Solutions to SemiLinear Equations with Quadratic Growth in the Gradient

Projected exponential methods for stiff dynamical low-rank approximation problems

On Nonlinear Matrix Equations from the First Standard Form

A Priori Estimates Versus Arbitrarily Large Solutions for Fractional Semi-Linear Elliptic Equations with Critical Sobolev Exponent

Inexact linear solves in the low-rank ADI iteration for large Sylvester equations