Natural Factor Based Solvers

06/09/2021
by   O. Andrés Cuervo, et al.
0

We consider parametric families of partial differential equations–PDEs where the parameter κ modifies only the (1,1) block of a saddle point matrix product of a discretization below. The main goal is to develop an algorithm that removes, as much as possible, the dependence of iterative solvers on the parameter κ. The algorithm we propose requires only one matrix factorization which does not depend on κ, therefore, allows to reuse it for solving very fast a large number of discrete PDEs for different κ and forcing terms. The design of the proposed algorithm is motivated by previous works on natural factor of formulation of the stiffness matrices and their stable numerical solvers. As an application, in two dimensions, we consider an iterative preconditioned solver based on the null space of Crouzeix-Raviart discrete gradient represented as the discrete curl of P_1 conforming finite element functions. For the numerical examples, we consider the case of random coefficient pressure equation where the permeability is modeled by an stochastic process. We note that contrarily from recycling Krylov subspace techniques, the proposed algorithm does not require fixed forcing terms.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

12/22/2020

Numerical approximation of singular-degenerate parabolic stochastic PDEs

We study a general class of singular degenerate parabolic stochastic par...
08/30/2019

Matrix equation techniques for certain evolutionary partial differential equations

We show that the discrete operator stemming from the time and space disc...
11/09/2020

Numerically Solving Parametric Families of High-Dimensional Kolmogorov Partial Differential Equations via Deep Learning

We present a deep learning algorithm for the numerical solution of param...
01/14/2022

Evaluating Accuracy and Efficiency of HPC Solvers for Sparse Linear Systems with Applications to PDEs

Partial Differential Equations (PDEs) describe several problems relevant...
05/23/2019

Recursive blocked algorithms for linear systems with Kronecker product structure

Recursive blocked algorithms have proven to be highly efficient at the n...
10/25/2019

A computational study of preconditioning techniques for the stochastic diffusion equation with lognormal coefficient

We present a computational study of several preconditioning techniques f...
07/13/2020

Discretization-error-accurate mixed-precision multigrid solvers

This paper builds on the algebraic theory in the companion paper [Algebr...
This week in AI

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