We are interested in numerically approximating the solution ${\bf U}(t)$ of the large dimensional semilinear matrix differential equation $\dot{\bf U}(t) = { \bf A}{\bf U}(t) + {\bf U}(t){ \bf B} + {\cal F}({\bf U},t)$, with appropriate starting and boundary conditions, and $ t \in [0, T_f]$. In the framework of the Proper Orthogonal Decomposition (POD) methodology and the Discrete Empirical Interpolation Method (DEIM), we derive a novel matrix-oriented reduction process leading to an effective, structure aware low order approximation of the original problem. The reduction of the nonlinear term is also performed by means of a fully matricial interpolation using left and right projections onto two distinct reduction spaces, giving rise to a new two-sided version of DEIM. By maintaining a matrix-oriented reduction, we are able to employ first order exponential integrators at negligible costs. Numerical experiments on benchmark problems illustrate the effectiveness of the new setting.
翻译:我们有兴趣从数字上接近大型半线性矩阵差异方程式的解决方案${bf U}(t) = {bf Aunbf U}(t) + bf U}(t) {bf B}+ cals F}(bf U}),t) $(t),加上适当的起始和边界条件和[0,T_f]美元。在正正正正正方形分解法(POD)和分解经验式内插方法(DEIM)的框架内,我们产生了一个新的面向矩阵的削减进程,导致一个有效的结构,意识到最初问题的低顺序近似值。非线性术语的削减还采用完全母体间插图的方式,利用左侧和右侧的预测,在两个不同的递减空间进行,从而产生一个新的双面版本的DEIM。通过保持一个面向矩阵的缩减,我们能够利用第一个顺序的指数化实验,在可计量的成本上说明基准问题。