An MSSS-preconditioned matrix equation approach for the time-harmonic elastic wave equation at multiple frequencies