Dual-pairing summation by parts finite difference methods for large scale elastic wave simulations in 3D complex geometries.