Our implementation is based on the constant-density acoustic wave equation

(1)

where is a pressure-field solution, are the model coordinates
( and for the two-dimensional case), and is time. Though the
implementation uses time-domain finite differences, it is more convenient to
express the equation in frequency and slowness :