Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aezo_v1_0.tar.gz(5591 Kbytes)
Manuscript Title: High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
Authors: Lukas Einkemmer
Program title: sldg
Catalogue identifier: AEZO_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 202(2016)326
Programming language: C++03.
Computer: PC, HPC systems.
Operating system: POSIX compatible (extensively tested on various Linux systems). In fact only the timing class requires POSIX routines; all other parts of the program can be run on any system where a C++ compiler, Boost, and FFTW is available.
RAM: Megabytes to Terabytes (depending on the problem size)
Keywords: Advection equations, High dimensional problems, C++, MPI, OpenMP, Vlasov solver.
PACS: 02.60.Lj, 02.60.-x, 52.65.-y.
Classification: 6.5, 4.3, 19.8.

External routines: Boost, FFTW

Nature of problem:
An approximate solution of the advection equation is computed using the semi-Lagrangian discontinuous Galerkin method. This procedure is used as the fundamental building block in a splitting based Vlasov solver.

Solution method:
The semi-Lagrangian discontinuous Galerkin approach is employed. For the Vlasov-Poisson problem we, in addition, use a splitting approach in time and the FFT to solve Poisson's equation.

We assume that the hyperbolic part of the problem under consideration can be (for example, using a splitting approach) decomposed into a sequence of one-dimensional advections.

Unusual features:
The framework is independent of the dimension of the problem under consideration. This is accomplished using template techniques.

Running time:
From minutes to weeks (depending on the problem size).