Natural Factor Based Solvers

by   O. Andrés Cuervo, et al.

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.



There are no comments yet.


page 1

page 2

page 3

page 4


Numerical approximation of singular-degenerate parabolic stochastic PDEs

We study a general class of singular degenerate parabolic stochastic par...

Matrix equation techniques for certain evolutionary partial differential equations

We show that the discrete operator stemming from the time and space disc...

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...

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

Partial Differential Equations (PDEs) describe several problems relevant...

Recursive blocked algorithms for linear systems with Kronecker product structure

Recursive blocked algorithms have proven to be highly efficient at the n...

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

We present a computational study of several preconditioning techniques f...

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.