The Ehrenfest+R method 10 was originally developed to simulate the electronic quantum subsystem coupled to the classical electromagnetic field in order to accurately describes spontaneous emission processes. It effectively captures Lindblad dynamics with a deterministic change of the magnitude of the quantum coefficients, and stochastic changes of the phases. While the L-MFE method is similar to (and largely inspired by) the Ehren-fest+R method, there is a key difference between Ehren-fest+R and L-MFE in the off-diagonal reduced density matrix element decay procedure. To facilitate the theoretical comparison, we briefly summarize the Ehren-fest+R method 10 in Appendix A.