Coarsening in Algebraic Multigrid using Gaussian Processes

04/23/2020
by   Hanno Gottschalk, et al.
0

Multigrid methods have proven to be an invaluable tool to efficiently solve large sparse linear systems arising in the discretization of partial differential equations (PDEs). Algebraic multigrid methods and in particular adaptive algebraic multigrid approaches have shown that multigrid efficiency can be obtained without having to resort to properties of the PDE. Yet the required setup of these methods poses a not negligible overhead cost. Methods from machine learning have attracted attention to streamline processes based on statistical models being trained on the available data. Interpreting algebraically smooth error as an instance of a Gaussian process, we develop a new, data driven approach to construct adaptive algebraic multigrid methods. Based on Gaussian a priori distributions, Kriging interpolation minimizes the mean squared error of the a posteriori distribution, given the data on the coarse grid. Going one step further, we exploit the quantification of uncertainty in the Gaussian process model in order to construct efficient variable splittings. Using a semivariogram fit of a suitable covariance model we demonstrate that our approach yields efficient methods using a single algebraically smooth vector.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/28/2018

Algorithmic Linearly Constrained Gaussian Processes

We algorithmically construct multi-output Gaussian process priors which ...
research
05/29/2014

Functional Gaussian processes for regression with linear PDE models

In this paper, we present a new statistical approach to the problem of i...
research
01/07/2020

On efficient numerical solution of linear algebraic systems arising in goal-oriented error estimates

We deal with the numerical solution of linear partial differential equat...
research
12/29/2022

Gaussian Process Priors for Systems of Linear Partial Differential Equations with Constant Coefficients

Partial differential equations (PDEs) are important tools to model physi...
research
03/25/2021

On the quantification of discretization uncertainty: comparison of two paradigms

Numerical models based on partial differential equations (PDE), or integ...
research
08/30/2022

A greedy MOR method for the tracking of eigensolutions to parametric elliptic PDEs

In this paper, we introduce an algorithm based on a sparse grid adaptive...
research
05/23/2019

A Bulirsch-Stoer algorithm using Gaussian processes

In this paper, we treat the problem of evaluating the asymptotic error i...

Please sign up or login with your details

Forgot password? Click here to reset