A data driven heuristic for rapid convergence of general Scheduled Relaxation Jacobi (SRJ) schemes

11/12/2020
by   Mohammad Shafaet Islam, et al.
0

The Scheduled Relaxation Jacobi (SRJ) method is a viable candidate as a high performance linear solver for elliptic PDEs. The method greatly improves the convergence of standard Jacobi iteration by applying a sequence of M overrelaxation and underrelaxation steps in each cycle of the algorithm. In previous work, the relaxation factors associated with each of the M steps (which characterize an SRJ scheme) were derived to be specific to the problem of interest and its discretization. In this work we develop a class of SRJ schemes which could be applied to solve any linear system as long as the original Jacobi iterative method would converge. Furthermore, we use data to train an algorithm to select which scheme to use at each cycle of the SRJ method for rapid convergence. Specifically, the algorithm is trained using convergence data obtained from randomly applying SRJ schemes to the 1D Poisson problem. The automatic selection heuristic that is developed based on this limited data is found to provide good convergence for a wide range of problems.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/13/2021

Scheduled Relaxation Jacobi schemes for non-elliptic partial differential equations

The Scheduled Relaxation Jacobi (SRJ) method is a linear solver algorith...
research
04/22/2022

Parameter-robust Braess-Sarazin-type smoothers for linear elasticity problems

In this work, we propose three Braess-Sarazin-type multigrid relaxation ...
research
12/18/2019

PCPATCH: software for the topological construction of multigrid relaxation methods

Effective relaxation methods are necessary for good multigrid convergenc...
research
10/08/2019

Implicit Neural Solver for Time-dependent Linear PDEs with Convergence Guarantee

Fast and accurate solution of time-dependent partial differential equati...
research
10/15/2021

Tweed and wireframe: accelerated relaxation algorithms for multigrid solution of elliptic PDEs on stretched structured grids

Two new relaxation schemes are proposed for the smoothing step in the ge...
research
11/14/2022

SlabLU: A Sparse Direct Solver for Elliptic PDEs on Rectangular Domains

The paper describes a sparse direct solver for the linear systems that a...

Please sign up or login with your details

Forgot password? Click here to reset