Extendible and Efficient Python Framework for Solving Evolution Equations with Stabilized Discontinuous Galerkin Methods

This paper discusses a Python interface for the recently published Dune-Fem-DG module which provides highly efficient implementations of the discontinuous Galerkin (DG) method for solving a wide range of nonlinear partial differential equations (PDEs). Although the C++ interfaces of Dune-Fem-DG are highly flexible and customizable, a solid knowledge of C++ is necessary to make use of this powerful tool. With this work, easier user interfaces based on Python and the unified form language are provided to open Dune-Fem-DG for a broader audience. The Python interfaces are demonstrated for both parabolic and first-order hyperbolic PDEs.