A modified next reaction method for simulating chemical systems with time dependent propensities and delays
Click to view
The above plot compares the speeds of the rejection method, direct method, and next reaction method for systems with delays as the percentage of timesteps that are rejected in the rejection method, as parametrized by , increases. For different values of , each method was used to simulate the system (18) times. The plot above gives the CPU time needed for each method as a function of . We see that the rejection and direct methods are nearly equivalent while the next reaction method for systems with delays is significantly more efficient than both for all .
Article metrics loading...
Full text loading...