High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code

DOI10.1016/J.CPC.2016.01.012zbMATH Open1348.35004arXiv1501.05508OpenAlexW1945958945MaRDI QIDQ339403FDOQ339403

L. Einkemmer

Publication date: 10 November 2016

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: The recently developed semi-Lagrangian discontinuous Galerkin approach is used to discretize hyperbolic partial differential equations (usually first order equations). Since these methods are conservative, local in space, and able to limit numerical diffusion, they are considered a promising alternative to more traditional semi-Lagrangian schemes (which are usually based on polynomial or spline interpolation). In this paper, we consider a parallel implementation of a semi-Lagrangian discontinuous Galerkin method for distributed memory systems (so-called clusters). Both strong and weak scaling studies are performed on the Vienna Scientific Cluster 2 (VSC-2). In the case of weak scaling, up to 8192 cores, we observe a parallel efficiency above 0.89 for both two and four dimensional problems. Strong scaling results show good scalability to at least 1024 cores (we consider problems that can be run on a single processor in reasonable time). In addition, we study the scaling of a two dimensional Vlasov--Poisson solver that is implemented using the framework provided. All of the simulation are conducted in the context of worst case communication overhead; i.e., in a setting where the CFL number increases linearly with the problem size. The framework introduced in this paper facilitates a dimension independent implementation (based on C++ templates) of scientific codes using both an MPI and a hybrid approach to parallelization. We describe the essential ingredients of our implementation.

Full work available at URL: https://arxiv.org/abs/1501.05508

zbMATH Keywords

parallelizationhigh-performance computingVlasov-Poisson equationssemi-Lagrangian discontinuous Galerkin methods

Mathematics Subject Classification ID

Parallel numerical computation (65Y05) Vlasov equations (35Q83) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Software, source code, etc. for problems pertaining to numerical analysis (65-04) Software, source code, etc. for problems pertaining to partial differential equations (35-04)

Cites Work

• Numerical study of blow up and stability of solutions of generalized Kadomtsev-Petviashvili equations
• Title not available (Why is that?)
• Title not available (Why is that?)
• Conservative semi-Lagrangian schemes for Vlasov equations
• Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: theoretical analysis and application to the Vlasov-Poisson system
• A parallel Vlasov solver based on local cubic spline interpolation on patches
• Numerical study of oscillatory regimes in the Kadomtsev-Petviashvili equation
• Discontinuous Galerkin semi-Lagrangian method for Vlasov-Poisson
• On the error propagation of semi-Lagrange and Fourier methods for advection problems
• Convergence Analysis of a Discontinuous Galerkin/Strang Splitting Approximation for the Vlasov--Poisson Equations
• A splitting approach for the Kadomtsev-Petviashvili equation

Cited In (13)

• The numerical flow iteration for the Vlasov-Poisson equation
• A Low-Rank Projector-Splitting Integrator for the Vlasov--Poisson Equation
• A hybrid semi-Lagrangian DG and ADER-DG solver on a moving mesh for Liouville's equation of geometrical optics
• A robust second-order low-rank BUG integrator based on the midpoint rule
• Dynamical Low-Rank Integrator for the Linear Boltzmann Equation: Error Analysis in the Diffusion Limit
• Efficient 6D Vlasov simulation using the dynamical low-rank framework \texttt{Ensign}
• Suppressing instability in a Vlasov-Poisson system by an external electric field through constrained optimization
• A Quasi-Conservative Dynamical Low-Rank Algorithm for the Vlasov Equation
• A low-rank projector-splitting integrator for the Vlasov-Maxwell equations with divergence correction
• A performance comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
• A split step Fourier/discontinuous Galerkin scheme for the Kadomtsev-Petviashvili equation
• Semi-Lagrangian Vlasov simulation on GPUs
• Exponential methods for solving hyperbolic problems with application to collisionless kinetic equations

Uses Software

• FFTW
• sldg
• GYSELA

This page was built for publication: High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q339403)