This work pertains to numerical aspects of a finite element method based discontinuous functions. Our study focuses on the Interior Penalty Discontinuous Galerkin method (IPDGM) because of its high-level of flexibility for solving the full wave equation in heterogeneous media. We assess the performance of IPDGM through a comparison study with a spectral element method (SEM). We show that IPDGM is as accurate as SEM. In addition, we illustrate the efficiency of IPDGM when employed in a seismic imaging process by considering two-dimensional problems involving the Reverse Time Migration.
