“…Seismic wavefield modeling is a key component of seismic imaging, e.g., reverse-time migration [3,27,33,49,50,52,53,60], and inversion, e.g., full-waveform inversion [36,37,43,48]. Due to the heterogeneity of the Earth, the simulation of seismic wavefields is achieved most commonly using numerical methods, such as the finite-element method [8,30,35], the finite-difference method [1,25,44,45], spectral-element method [23,24], or discontinuous Galerkin [6,7,11,15,21].…”