The discontinuous Galerkin method in the Nmesh program

We present an introduction to the new Nmesh code, which uses a discontinuous Galerkin (DG) method that is intended to parallelize well. We explain our implementation of the DG method. In particular, we discuss a simplification in the normalization of the domain boundary normals. We are able to normalize them with a flat metric, even when the physical metric is far from flat. We also present recent results from single neutron star evolutions with and without a perturbation. Since the star matter is not smooth across the neutron star surface, this is challenging for a DG method. To address this problem we have developed a positivity limiter that together with exponential filters is enough to stabilize the evolution.