Whether the cumulative sum is stored separately inside the provided algorithm struct. In this case, a constructor is provided, which recieves the amount of propensietes during the simulation. An array of proper type is expanded to store the cumulative fold resuls, which can be search faster.
Whether the Sandmann enhancement is turned on
The algorithms should work with all kinds of floating numbers as well as with integers.
The convenience function to release a gillespie algorithm structure with some preset properties. In this case, the preset handles the memory allocation due using an internal buffer for storing passed propensities. Other properties do not affect memory allocations.
The floating type in use is defined either by the user on algorithm declaration or is chosen to be real, in case propensities come in form of integer values.
The module contains various physical time reaction propagation algorithms based on the Gillespie algorithm from 1977. In the current version this module contains - the default version of Gillespie - The modification introduced by Sandmann 2008, which makes possible faster sample reaction arrival times, in case all propensities of the reactions are known exactly. I. e. this modification can be applied, in case propensities were not estimated. - The modification introduced by Petzold et. al. 2006. At the expense of allocating an ancillary propensities array, where their cumulative sum is stored, the velocity of sampling the reaction index increases, dependent on the chosen search policy of assumedSorted ranges. - The module provides also a common testing routine, containing tests for all combinations of modifications.