Efficient and accurate modeling of the wave equation is important for seismic exploration. We compare the performance of explicit time stepping with the discontinuous Galerkin method and with continuous augmented third and fourth-order mass-lumped elements for tetrahedra. There are two choices for the last, one with a more favorable CFL number. Numerical experiments illustrate the accuracy, usefulness, and versatility of these methods when solving 3D problems in inhomogeneous media.