Numerical Issues in the Implementation of High Order Polynomial Multi-Domain Penalty Spectral Galerkin Methods for Hyperbolic Conservation Laws