DeepAI AI Chat
Log In Sign Up

Scheduled Relaxation Jacobi schemes for non-elliptic partial differential equations

by   Mohammad Shafaet Islam, et al.

The Scheduled Relaxation Jacobi (SRJ) method is a linear solver algorithm which greatly improves the convergence of the Jacobi iteration through the use of judiciously chosen relaxation factors (an SRJ scheme) which attenuate the solution error. Until now, the method has primarily been used to accelerate the solution of elliptic PDEs (e.g. Laplace, Poisson's equation) as the currently available schemes are restricted to solving this class of problems. The goal of this paper is to present a methodology for constructing SRJ schemes which are suitable for solving non-elliptic PDEs (or equivalent, nonsymmetric linear systems arising from the discretization of these PDEs), thereby extending the applicability of this method to a broader class of problems. These schemes are obtained by numerically solving a constrained minimization problem which guarantees the solution error will not grow as long as the linear system has eigenvalues which lie in certain regions of the complex plane. We demonstrate that these schemes are able to accelerate the convergence of standard Jacobi iteration for the nonsymmetric linear systems arising from discretization of the 1D and 2D steady advection-diffusion equations.


page 7

page 13


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

The Scheduled Relaxation Jacobi (SRJ) method is a viable candidate as a ...

Improving Solve Time of aggregation-based adaptive AMG

This paper proposes improving the solve time of a bootstrap AMG designed...

Domain Decomposition for the Closest Point Method

The discretization of elliptic PDEs leads to large coupled systems of eq...

A deep first-order system least squares method for solving elliptic PDEs

We propose a First-Order System Least Squares (FOSLS) method based on de...

Minimal positive stencils in meshfree finite difference methods for linear elliptic equations in non-divergence form

We design a monotone meshfree finite difference method for linear ellipt...

Convergence and Semi-convergence of a class of constrained block iterative methods

In this paper, we analyze the convergence projected non-stationary bloc...