Abstract:This paper presents the parallelization of two widely used implicit numerical solvers for the solution of partial differential equations on structured meshes, namely, the ADI (Alternating-Direction Implicit) solver for tridiagonal linear systems and the SIP (Strongly Implicit Procedure) solver for the penta-diagonal systems. Both solvers were parallelized using CUDA (Computer Unified Device Architecture) Fortran on GPGPUs (General-Purpose Graphics Processing Units). The parallel ADI solver (P-ADI) is based on the Parallel Cyclic Reduction (PCR) algorithm, while the parallel SIP solver (P-SIP) uses the wave front method (WF) following a diagonal line calculation strategy. To map the solution schemes onto the hierarchical block-threads framework of the CUDA on the GPU, the P-ADI solver adopted two mapping methods, one block thread with iterations (OBM-it) and multi-block threads (MBMs), while the P-SIP solver also used two mappings, one conventional mapping using effective WF lines (WF-e) with matrix coefficients and solution variables defined on original computational mesh, and a newly proposed mapping using all WF mesh (WF-all), on which matrix coefficients and solution variables are defined. Both the P-ADI and the P-SIP have been integrated into a two-dimensional (2D) hydrodynamic model, the CCHE2D (Center of Computational Hydroscience and Engineering) model, developed by the National Center for Computational Hydroscience and Engineering at the University of Mississippi. This study for the first time compared these two parallel solvers and their efficiency using examples and applications in complex geometries, which can provide valuable guidance for future uses of these two parallel implicit solvers in computational fluids dynamics (CFD). Both parallel solvers demonstrated higher efficiency than their serial counterparts on the CPU (Central Processing Unit): 3.73~4.98 speedup ratio for flow simulations, and 2.166~3.648 speedup ratio for sediment transport simulations. In general, the P-ADI solver is faster than but not as stable as the P-SIP solver; and for the P-SIP solver, the newly developed mapping method WF-all significantly improved the conventional mapping method WF-e.

The effective use of BLAS interface for implementation of finite-element ADER-DG and finite-volume ADER-WENO methods

Efficient implementation of ADER discontinuous Galerkin schemes for a scalable hyperbolic PDE engine

High order ADER-DG method with local DG predictor for solutions of differential-algebraic systems of equations

Arbitrary high order ADER-DG method with local DG predictor for solutions of initial value problems for systems of first-order ordinary differential equations

Developing a BLAS library for the AMD AI Engine

On improving the efficiency of ADER methods

Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell WENO finite-volume limiting for simulation of non-stationary compressible multicomponent reactive flows

Matrix-free approaches for GPU acceleration of a high-order finite element hydrodynamics application using MFEM, Umpire, and RAJA

An Efficient Implementation of Automatic Differentiation in Interior Point Optimal Power Flow

FT-BLAS: A Fault Tolerant High Performance BLAS Implementation on x86 CPUs

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Interface for Sparse Linear Algebra Operations

A Semi-Lagrangian Adaptive-Rank (SLAR) Method for Linear Advection and Nonlinear Vlasov-Poisson System

A compact simple HWENO scheme with ADER time discretization for hyperbolic conservation laws I: Structured meshes

Implementation of strong implicit procedure for the energy equations in subchannel code athas

Toward a Simple and Accurate Lagrangian-based Error Estimator for the BDF Algorithms and Adaptive Time-Stepping

Improved ADI Iterative Algorithm With Gauss-Seidel Ideology for Efficient WLP-FDTD Method in 3-D Cylindrical Coordinate System

Using AI libraries for Incompressible Computational Fluid Dynamics