Abstract:We develop two new sets of stable, rank-adaptive Dynamically Orthogonal Runge-Kutta (DORK) schemes that capture the high-order curvature of the nonlinear low-rank manifold. The DORK schemes asymptotically approximate the truncated singular value decomposition at a greatly reduced cost while preserving mode continuity using newly derived retractions. We show that arbitrarily high-order optimal perturbative retractions can be obtained, and we prove that these new retractions are stable. In addition, we demonstrate that repeatedly applying retractions yields a gradient-descent algorithm on the low-rank manifold that converges superlinearly when approximating a low-rank matrix. When approximating a higher-rank matrix, iterations converge linearly to the best low-rank approximation. We then develop a rank-adaptive retraction that is robust to overapproximation. Building off of these retractions, we derive two rank-adaptive integration schemes that dynamically update the subspace upon which the system dynamics are projected within each time step: the stable, optimal Dynamically Orthogonal Runge-Kutta (so-DORK) and gradient-descent Dynamically Orthogonal Runge-Kutta (gd-DORK) schemes. These integration schemes are numerically evaluated and compared on an ill-conditioned matrix differential equation, an advection-diffusion partial differential equation, and a nonlinear, stochastic reaction-diffusion partial differential equation. Results show a reduced error accumulation rate with the new stable, optimal and gradient-descent integrators. In addition, we find that rank adaptation allows for highly accurate solutions while preserving computational efficiency.

Stabilization of a matrix via a low-rank-adaptive ODE

Stabilization of a matrix via a low rank-adaptive ODE

Robust Decentralised Stabilisation of Interval Discrete-Time Singular Large-Scale Systems

Robust D Stabilization of Singular Systems with Polytopic Uncertainties

Robust Stability and Stabilisation of Discrete-Time Descriptor Systems with Uncertainties in the Difference Matrix

Static Output Feedback Stabilization Based on Structured Lyapunov Matrix

A low rank ODE for spectral clustering stability

Implementation improvements and extensions of an ODE-based algorithm for structured low-rank approximation

Preconditioned Low-Rank Riemannian Optimization for Symmetric Positive Definite Linear Matrix Equations

Low-Rank Matrix Approximation with Stability.

Adaptive Low-Rank Methods: Problems on Sobolev Spaces

An Analysis of a Multi-Level Projected Steepest Descent Iteration for Nonlinear Inverse Problems in Banach Spaces Subject to Stability Constraints

A rank-adaptive robust integrator for dynamical low-rank approximation

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

A Relaxed Algorithm for Online Matrix Inversion

Projected exponential methods for stiff dynamical low-rank approximation problems

A Quadratically Convergent Algorithm for Structured Low-Rank Approximation

Optimization schemes on manifolds for structured matrices with fixed eigenvalues

Stable rank-adaptive Dynamically Orthogonal Runge-Kutta schemes

Novel Alternating Update Method for Low Rank Approximation of Structured Matrices

A non-damped stabilization algorithm for multibody dynamics