Abstract:The International Journal of High Performance Computing Applications, Ahead of Print. Recently, global and local relaxation Runge–Kutta methods have been developed for guaranteeing the conservation, dissipation, or other solution properties for general convex functionals whose dynamics are crucial for an ordinary differential equation solution. These novel time integration procedures have an application in a wide range of problems that require dynamics-consistent and stable numerical methods. The application of a relaxation scheme involves solving scalar nonlinear algebraic equations to find the relaxation parameter. Even though root-finding may seem to be a problem technically straightforward and computationally insignificant, we address the problem at scale as we solve full-scale industrial problems on a CPU-powered supercomputer and show its cost to be considerable. In particular, we apply the relaxation schemes in the context of the compressible Navier–Stokes equations and use them to enforce the correct entropy evolution. We use seven different algorithms to solve for the global and local relaxation parameters and analyze their strong scalability. As a result of this analysis, within the global relaxation scheme, we recommend using Brent's method for problems with a low polynomial degree and of small sizes for the global relaxation scheme, while secant proves to be the best choice for higher polynomial degree solutions and large problem sizes. For the local relaxation scheme, we recommend secant. Further, we compare the schemes' performance using their most efficient implementations, where we look at their effect on the timestep size, overhead, and weak scalability. We show the global relaxation scheme to be always more expensive than the local approach—typically 1.1–1.5 times the cost. At the same time, we highlight scenarios where the global relaxation scheme might underperform due to its increased communication requirements. Finally, we present an analysis that sets expectations on the computational overhead anticipated based on the system properties.

A data driven heuristic for rapid convergence of Scheduled Relaxation Jacobi schemes

Acceleration of the Scheduled Relaxation Jacobi Method: Promising Strategies for Solving Large, Sparse Linear Systems

On the Convergence of the Accelerated Riccati Iteration Method.

Performance analysis of relaxation Runge–Kutta methods

Adaptive Relaxation Strategy on Basic Iterative Methods for Solving Linear Systems with Single and Multiple Right-Hand Sides

Energy Conservative Relaxation-Free Runge-Kutta Schemes

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

The relaxing schemes for Hamilton-Jacobi equations

Splitting Schemes for Coupled Differential Equations: Block Schur-Based Approaches and Partial Jacobi Approximation

Stabilized Times Schemes for High Accurate Finite Differences Solutions of Nonlinear Parabolic Equations

A Class of High Resolution Difference Schemes for Nonlinear Hamilton--Jacobi Equations with Varying Time and Space Grids

Implicit Block Data-Parallel Relaxation Scheme of Navier–Stokes Equations Using Graphics Processing Units

A Comparative Study of Iterative Riemann Solvers for the Shallow Water and Euler Equations

Optimal a priori error estimate of relaxation-type linear finite element method for nonlinear Schrödinger equation

A PARALLEL DECOUPLING TECHNIQUE TO ACCELERATE CONVERGENCE OF RELAXATION SOLUTION OF INTEGRAL-DIFFERENTIAL-ALGEBRAIC EQUATIONS

The Parallel Waveform Relaxation Stochastic Runge–Kutta Method for Stochastic Differential Equations

Efficient Same-Dimensional Implicit Time Advancement Parallel Scheme and Optimization Methods for the Iteration Parameters Using a Graphics-Processing Unit

An Element-wise RSAV Algorithm for Unconstrained Optimization Problems

Enhanced shifted Jacobi operational matrices of integrals: spectral algorithm for solving some types of ordinary and fractional differential equations

A Symmetric Multigrid-Preconditioned Krylov Subspace Solver for Stokes Equations

Learning to Relax: Setting Solver Parameters Across a Sequence of Linear System Instances