  This paper presents HDGlab, an open source MATLAB implementation of the
hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a
detailed description of both the HDG method for elliptic problems and its
implementation available in HDGlab. Ultimately, this is expected to make this
relatively new advanced discretisation method more accessible to the
computational engineering community. HDGlab presents some features not
available in other implementations of the HDG method that can be found in the
free domain. First, it implements high-order polynomial shape functions up to
degree nine, with both equally-spaced and Fekete nodal distributions. Second,
it supports curved isoparametric simplicial elements in two and three
dimensions. Third, it supports non-uniform degree polynomial approximations and
it provides a flexible structure to devise degree adaptivity strategies.
Finally, an interface with the open-source high-order mesh generator Gmsh is
provided to facilitate its application to practical engineering problems.
