Julia package FractionalTimeDG.jl: Generate coefficient arrays needed for discontinuous Galerkin time-stepping of fractional diffusion problems. Paper: Implementation of high-order, discontinuous Galerkin time stepping for fractional diffusion problems. The discontinuous Galerkin dG method provides a robust and flexible technique for the time integration of fractional diffusion problems. However, a practical implementation uses coefficients defined by integrals that are not easily evaluated. We describe specialised quadrature techniques that efficiently maintain the overall accuracy of the dG method. In addition, we observe in numerical experiments that known superconvergence properties of dG time stepping for classical diffusion problems carry over in a modified form to the fractional-order setting.

Keywords for this software

Anything in here will be replaced on browsers that support the canvas element