DeepAI AI Chat
Log In Sign Up

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

by   Mohammad Shafaet Islam, et al.

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.


page 1

page 2

page 3

page 4


Scheduled Relaxation Jacobi schemes for non-elliptic partial differential equations

The Scheduled Relaxation Jacobi (SRJ) method is a linear solver algorith...

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

In this work, we propose three Braess-Sarazin-type multigrid relaxation ...

PCPATCH: software for the topological construction of multigrid relaxation methods

Effective relaxation methods are necessary for good multigrid convergenc...

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

Fast and accurate solution of time-dependent partial differential equati...

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

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