Motivation / Current Behaviour
As of now SDE models have to implement checks for whether flows/derivatives result in overflows in the compartments (negative compartments or compartments greater than population size). The current implementation can also result in undefined behavior by floating point operations (see Issue #1009)
Enhancement description
Adapting the EulerCore (potentially writing a EulerMaruyamaCore) one can use vectorized functions to add those safeguards and other modifications into the core. Potential modifications are:
- Clipping compartments to the intervall [0, population_size]
- Setting compartments below some threshhold to zero
- Rescaling compartments in case clipping changed total population size
Additional context
No response
Checklist
Motivation / Current Behaviour
As of now SDE models have to implement checks for whether flows/derivatives result in overflows in the compartments (negative compartments or compartments greater than population size). The current implementation can also result in undefined behavior by floating point operations (see Issue #1009)
Enhancement description
Adapting the EulerCore (potentially writing a EulerMaruyamaCore) one can use vectorized functions to add those safeguards and other modifications into the core. Potential modifications are:
Additional context
No response
Checklist