The old implementation used to result in stochastic order error in
debugmode. After many attempts to solve it, I decided that it would be
better and faster just to rewrite it. This new implementation does not
suffer from any bug (i.e. all tests pass in debug mode).