Before revision:  Moreover, our numerical experiments demonstrate that our new method is often several times faster than the standard optimization-based scheme.
After revision:  Moreover, a discrete entropy law can be enforced by modifying the time stepping scheme. Our numerical experiments demonstrate that our new method is often several times faster than the standard optimization-based scheme.