I present the first known numerical solver for fixed high-frequency wave propagation, which is at the same time fast (linear runtime) and accurate (second-order pollution-free), when the wave source and velocity background both have geophysically relevant singularities.