Multigrid wave ray algorithm, proposed by Brandt and Livshits in 1997, exhibits fast convergence and good scalability for the Helmholtz equations with large constant wave numbers, and until lately was deemed not extendable to problems with variable, in particular discontinuous, coefficients. In this talk we present first results of such extension, both in algebraic and geometric frameworks.