Can you put the actual simulation into a `if __name__ == '__main__'`  so we can import the simulation functions?
I had made this change locally so I can just import the functions from this module and simulate a specific example without running the entire simulations.