Deep Learning of Preconditioners for Conjugate Gradient Solvers in Urban Water Related Problems

06/17/2019
by   Johannes Sappl, et al.
0

Solving systems of linear equations is a problem occuring frequently in water engineering applications. Usually the size of the problem is too large to be solved via direct factorization. One can resort to iterative approaches, in particular the conjugate gradients method if the matrix is symmetric positive definite. Preconditioners further enhance the rate of convergence but hitherto only handcrafted ones requiring expert knowledge have been used. We propose an innovative approach employing Machine Learning, in particular a Convolutional Neural Network, to unassistedly design preconditioning matrices specifically for the problem at hand. Based on an in-depth case study in fluid simulation we are able to show that our learned preconditioner is able to improve the convergence rate even beyond well established methods like incomplete Cholesky factorization or Algebraic MultiGrid.

READ FULL TEXT
research
02/18/2022

Optimization of the Sparse Multi-Threaded Cholesky Factorization for A64FX

Sparse linear algebra routines are fundamental building blocks of a larg...
research
01/28/2021

Two-level Nyström–Schur preconditioner for sparse symmetric positive definite matrices

Randomized methods are becoming increasingly popular in numerical linear...
research
11/21/2022

Hierarchical LU preconditioning for the time-harmonic Maxwell equations

The time-harmonic Maxwell equations are used to study the effect of elec...
research
10/07/2021

A Hybrid Direct-Iterative Method for Solving KKT Linear Systems

We propose a solution strategy for linear systems arising in interior me...
research
06/22/2021

Toward a new fully algebraic preconditioner for symmetric positive definite problems

A new domain decomposition preconditioner is introduced for efficiently ...
research
07/19/2021

A Robust Algebraic Domain Decomposition Preconditioner for Sparse Normal Equations

Solving the normal equations corresponding to large sparse linear least-...
research
05/08/2023

Parallel Cholesky Factorization for Banded Matrices using OpenMP Tasks

Cholesky factorization is a widely used method for solving linear system...

Please sign up or login with your details

Forgot password? Click here to reset