Scheduled Relaxation Jacobi method: improvements and applications

Abstract

Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficiency in the reduction of the residual increases with the number of levels employed in the algorithm. Applying the original methodology to compute the algorithm parameters with more than 5 levels notably hinders obtaining optimal SRJ schemes, as the mixed (non-linear) algebraic-differential equations from which they result become notably stiff. Here we present a new methodology for obtaining the parameters of SRJ schemes that overcomes the limitations of the original algorithm and provide parameters for SRJ schemes with up to 15 levels and resolutions of up to $2^{15}$ points per dimension, allowing for acceleration factors larger than several hundreds with respect to the Jacobi method for typical resolutions and, in some high resolution cases, close to 1000. Furthermore, we extend the original algorithm to apply it to certain systems of non-linear ePDEs.

37 pages, 8 figures, submitted to JCP

Cited by in corpus (13)

• Towards asteroseismology of core-collapse supernovae with gravitational-wave observations - II. Spacetime perturbations
• Topological Defects and nano-Hz Gravitational Waves in Aligned Axion Models
• Estimation of the mechanical properties of the eye through the study of its vibrational modes
• Numerically solving the relativistic Grad-Shafranov equation in Kerr spacetimes: Numerical techniques
• Hyperbolic Relaxation Method for Elliptic Equations
• On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method
• On the existence of a luminosity threshold of GRB jets in massive stars
• Modelling 3D Force-Free Neutron Star Magnetospheres
• Global optimization for data assimilation in landslide tsunamis models
• Speeding up a few orders of magnitude the Jacobi method: high order Chebyshev-Jacobi over GPUs
• Numerical analysis of shell-based geometric image inpainting algorithms and their semi-implicit extension
• Non-Hydrostatic Pressure Shallow Flows: GPU Implementation Using Finite Volume and Finite Difference Scheme
• Constraint-satisfying binary boson star initial data via XCFC