This talk presents a strategy for computational wave propagation that consists in decomposing the solution wavefield onto a largely incomplete set of eigenfunctions of the weighted Laplacian, with eigenvalues chosen randomly. The recovery method is the ell-1 minimization of compressed sensing. For the mathematician, we establish three possibly new estimates for the wave equation that guarantee accuracy of the numerical method in one spatial dimension. For the engineer, the compressive strategy offers a unique combination of parallelism and memory savings that should be of particular relevance to applications in reflection seismology. Joint work with Gabriel Peyre.