a b s t r a c tThis is the first in a series of papers on implementing a discontinuous Galerkin (DG) method as an open source MATLAB/GNU Octave toolbox. The intention of this ongoing project is to provide a rapid prototyping package for application development using DG methods. The implementation relies on fully vectorized matrix/vector operations and is carefully documented; in addition, a direct mapping between discretization terms and code routines is maintained throughout. The present work focuses on a two-dimensional time-dependent diffusion equation with space/time-varying coefficients. The spatial discretization is based on the local discontinuous Galerkin formulation. Approximations of orders zero through four based on orthogonal polynomials have been implemented; more spaces of arbitrary type and order can be easily accommodated by the code structure.
This is the second in a series of papers on implementing a discontinuous Galerkin (DG) method as an open source MATLAB / GNU Octave toolbox. The intention of this ongoing project is to offer a rapid prototyping package for application development using DG methods. The implementation relies on fully vectorized matrix / vector operations and is comprehensively documented. Particular attention was paid to maintaining a direct mapping between discretization terms and code routines as well as to supporting the full code functionality in GNU Octave. The present work focuses on a two-dimensional time-dependent linear advection equation with space / time-varying coefficients, and provides a general order implementation of several slope limiting schemes for the DG method.
The third paper in our series on open source MATLAB / GNU Octave implementation of the discontinuous Galerkin (DG) method(s) focuses on a hybridized formulation. The main aim of this ongoing work is to develop rapid prototyping techniques covering a range of standard DG methodologies and suitable for small to medium sized applications. Our FESTUNG package relies on fully vectorized matrix / vector operations throughout, and all details of the implementation are fully documented. Once again, great care is taken to maintain a direct mapping between discretization terms and code routines as well as to ensure full compatibility to GNU Octave. The current work formulates a hybridized DG scheme for a linear advection problem, describes hybrid approximation spaces on the mesh skeleton, and compares the performance of this discretization to the standard (element-based) DG method for different polynomial orders.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.