An Efficient Implementation of the Divergence Free Constraint in a Discontinuous Galerkin Method for Magnetohydrodynamics on Unstructured Meshes