Abstract We develop a new finite-difference scheme that reduces the number of grid points per wavelength required in frequency-domain elastic modeling. Our approach computes weighted averages of the spatial second-order derivative and the mass acceleration terms using a 25-point computational stenci...