The conventional finite-difference (FD) method often suffers from numerical dispersion when too few samples per wavelength are used or when models have large velocity contrast, or artefacts caused by source at grid points. ...