Towards Neural Sparse Linear Solvers

03/14/2022
by   Luca Grementieri, et al.
0

Large sparse symmetric linear systems appear in several branches of science and engineering thanks to the widespread use of the finite element method (FEM). The fastest sparse linear solvers available implement hybrid iterative methods. These methods are based on heuristic algorithms to permute rows and columns or find a preconditioner matrix. In addition, they are inherently sequential, making them unable to leverage the GPU processing power entirely. We propose neural sparse linear solvers, a deep learning framework to learn approximate solvers for sparse symmetric linear systems. Our method relies on representing a sparse symmetric linear system as an undirected weighted graph. Such graph representation is inherently permutation-equivariant and scale-invariant, and it can become the input to a graph neural network trained to regress the solution. We test neural sparse linear solvers on static linear analysis problems from structural engineering. Our method is less accurate than classic algorithms, but it is hardware-independent, fast on GPUs, and applicable to generic sparse symmetric systems without any additional hypothesis. Although many limitations remain, this study shows a general approach to tackle problems involving sparse symmetric matrices using graph neural networks.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/12/2020

Learning Algebraic Multigrid Using Graph Neural Networks

Efficient numerical solvers for sparse linear systems are crucial in sci...
research
01/14/2022

StAnD: A Dataset of Linear Static Analysis Problems

Static analysis of structures is a fundamental step for determining the ...
research
10/13/2017

On Parallel Solution of Sparse Triangular Linear Systems in CUDA

The acceleration of sparse matrix computations on modern many-core proce...
research
06/25/2021

Linear solvers for power grid optimization problems: a review of GPU-accelerated linear solvers

The linear equations that arise in interior methods for constrained opti...
research
11/10/2019

MATLAB Implementation of Element-based Solvers

Rahman and Valdman (2013) introduced a vectorized way to assemble finite...
research
09/15/2020

A two-level iterative scheme for general sparse linear systems based on approximate skew-symmetrizers

We propose a two-level iterative scheme for solving general sparse linea...
research
09/04/2020

Hybrid DCOP Solvers: Boosting Performance of Local Search Algorithms

We propose a novel method for expediting both symmetric and asymmetric D...

Please sign up or login with your details

Forgot password? Click here to reset