MARS-M: When Variance Reduction Meets Matrices