Based on a landmark paper by Pao and Gajewski, some novel developments of the method of generalized ray integrals are discussed. The expansion of the dynamic Green's function of the infinite space into plane waves allows benchmark 3-D solutions in the layered half-space and even enters the background formulation of elastic-viscoplastic wave propagation. New developments of software of combined symbolic-numerical manipulation and parallel computing make the method a competitive solution technique.