Please see the Jupyter notebooks for the main use of the code to generate the figures found in the paper.
The main algorithms are implemented in simulate.py and inverse_moments.py.
For further questions, feel free to email me (Pat Laub), I'll be happy to discuss anything related to this codebase.
Reference:
Young Lee, Patrick J. Laub, Thomas Taimre, Hongbiao Zhao, Jiancang Zhuang (2021), Exact simulation of extrinsic stress-release processes, Journal of Applied Probability, to appear in 59(1) in 2022.