A hybrid numerical scheme that simultaneously retrieves single-crystal elastic constants, C-ij, and wave-normal directions is applied to determine the elasticity of natural samples of both magnesite and dolomite, as measured by Brillouin spectroscopy at a