Found the bug and fixed it. First reaction method works for non-exponential distributions now because enabling times are clearer.