Score: 0

Interpolatory Dynamical Low-Rank Approximation: Theoretical Foundations and Algorithms

Published: October 22, 2025 | arXiv ID: 2510.19518v1

By: Benjamin Carrel , Daniel Kressner , Hei Yin Lam and more

Potential Business Impact:

Makes computer models run much faster.

Business Areas:
Big Data Data and Analytics

Dynamical low-rank approximation (DLRA) is a widely used paradigm for solving large-scale matrix differential equations, as they arise, for example, from the discretization of time-dependent partial differential equations on tensorized domains. Through orthogonally projecting the dynamics onto the tangent space of a low-dimensional manifold, DLRA achieves a significant reduction of the storage required to represent the solution. However, the need for evaluating the velocity field can make it challenging to attain a corresponding reduction of computational cost in the presence of nonlinearities. In this work, we address this challenge by replacing orthogonal tangent space projections with oblique, data-sparse projections selected by a discrete empirical interpolation method (DEIM). At the continuous-time level, this leads to DLRA-DEIM, a well-posed differential inclusion (in the Filippov sense) that captures the discontinuities induced by changes in the indices selected by DEIM. We establish an existence result, exactness property and error bound for DLRA-DEIM that match existing results for DLRA. For the particular case of QDEIM, a popular variant of DEIM, we provide an explicit convex-polytope characterization of the differential inclusion. Building on DLRA-DEIM, we propose a new class of projected integrators, called PRK-DEIM, that combines explicit Runge--Kutta methods with DEIM-based projections. We analyze the convergence order of PRK-DEIM and show that it matches the accuracy of previously proposed projected Runge-Kutta methods, while being significantly cheaper. Extensions to exponential Runge--Kutta methods and low-order tensor differential equations demonstrate the versatility of our framework.

Country of Origin
🇨🇭 Switzerland

Page Count
33 pages

Category
Mathematics:
Numerical Analysis (Math)