Hello. Need help for the details of Brettel's simulation algorithm.

[wilixx]

I implement the algorithm using python as described in your readme file. But the results seems to be different. Could you give me more details about Brettel's algorithm? Actually I have got several algorithms claimed to be Brettel's. Many of them are proved to be wrong, hopefully you are right. Thank you~