The simulation of wavefield propagation in anisotropic media is essential to produce subsurface images.Low-rank approximation allows one to obtain a numerical approximation of wavefield evolution operator that honors the kinematics and approximates amplitudes in heterogeneous anisotropic media.Here we discuss our implementation of low-rank approximation for anisotropic wave propagation in 2D and 3D.The main feature of our investigation is the use of row action methods to considerably reduce storage when computing …