Before revision:  Several numerical experiments based on typical benchmark problems illustrate the effectiveness of the new matrix-oriented setting, also for coupled systems of nonlinear matrix differential equations .
After revision:  By maintaining a matrix-oriented reduction, we are able to employ first order exponential integrators at negligible costs. Numerical experiments on benchmark problems illustrate the effectiveness of the new matrix-oriented setting, also for coupled systems of nonlinear matrix differential equations .