Preconditioners for model order reduction by interpolation and random sketching of operators

by   Oleg Balabanov, et al.

The performance of projection-based model order reduction methods for solving parameter-dependent systems of equations highly depends on the properties of the operator, which can be improved by preconditioning. In this paper we present strategies to construct a parameter-dependent preconditioner by an interpolation of operator's inverse. The interpolation is obtained by minimizing a discrepancy between the (preconditioned) operator and the matrix defining the metric of interest. The discrepancy measure is chosen such that its minimization can be efficiently performed online for each parameter value by the solution of a small least-squares problem. Furthermore, we show how to tune the discrepancy measure for improving the quality of Petrov-Galerkin projection or residual-based error estimation. This paper also addresses preconditioning for the randomized model order reduction methods from [Balabanov and Nouy 2019, Part I]. Our methodology can be readily used for efficient and stable solution of ill-conditioned parametric systems and an effective error estimation/certification without the need to estimate expensive stability constants. The proposed approach involves heavy computations in both offline and online stages that are circumvented by random sketching. The norms of high-dimensional matrices and vectors are estimated by l2-norms of their low-dimensional images, called sketches, through random embeddings. For this we extend the framework from [Balabanov and Nouy 2019, Part I] to random embeddings of operators.



page 1

page 2

page 3

page 4


Randomized linear algebra for model reduction. Part II: minimal residual methods and dictionary-based approximation

A methodology for using random sketching in the context of model order r...

Basic Ideas and Tools for Projection-Based Model Reduction of Parametric Partial Differential Equations

We provide first the functional analysis background required for reduced...

Interpolation Operator on negative Sobolev Spaces

We introduce a Scott–Zhang type projection operator mapping to Lagrange ...

Parametric model reduction via rational interpolation along parameters

We present a novel projection-based model reduction framework for parame...

Tractability properties of the discrepancy in Orlicz norms

We show that the minimal discrepancy of a point set in the d-dimensional...

A Differentiable Solver Approach to Operator Inference

Model Order Reduction is a key technology for industrial applications in...

Randomized Gram-Schmidt process with application to GMRES

A randomized Gram-Schmidt algorithm is developed for orthonormalization ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.