Problematic results for a simple reversible autocatalytic reaction

[PENG-Zhen]

Consider a simple reversible autocatalytic reaction: A+B<=>2A, forward_stochastic_reaction_constant = reverse_stochastic_reaction_constant = 0.01. Let this reaction run in a closed reactor, starting with 999 B's and 1 A, and the expected result should be that both A and B counts oscillate around 500. However, the actual result is that A count oscillates around 670 while B count oscillates around 330. It is likely that there is something wrong with the module that deals with stoichiometry.
I also suspect that reaction systems with similar reactions, such as the Brusselator (e.g., 2X+Y=>3X) and the Oregonator (e.g., A+X=>2X+2Z), will face similar issues.
I am gonna use my own version of GA script for now, which seems not as fast as yours but the results look correct. Hopefully the issue will be fixed in your next update.