nijawa

I use std::min and std::max instead of std::clamp to ensure that flows are in the correct bonds. Usage of std::clamp gives a "invalid bounds arguments passed to std::clamp" error which...

@mknaranja At `t = 21.5` we have the compartment InfectedV1V2 at `0.014175747933863577` with an ingoing flow of `0` and an outgoing flow of `0.14175747933863578`. With `step size = 0.1` that...

Clamping/maximizing with `compartment / step_size - eps` might work

my seir2v model is not yet in main yet so if I branch from main that model will be missing I think? There also still is a bit I have...

Clamping with a tolerance also is problematic if compartments are 0. To fix this I could either use a multiplicative tolerance (something like 1-1e-10) or use std::max(std::min to enforce non-negative...

> I suggest clamping to zero exactly, and only use tolerances on the upper bound. Then you don't need a relative tolerance, at least I would've just subtracted `tol`. You...

@xsaschako @mknaranja Attached you can find the error message when I try to install with Python 3.10.11 and VS [fehler_pycode.txt](https://github.com/SciCompMod/memilio/files/13907026/fehler_pycode.txt)

I have proposed a possible "new" EulerCore in the linked branch, before opening a pull request I would like to discuss whether this should be moved into a new core

Hi René thank you for your insight. I will implement a model with 3 variants and 4 different immunity states as an example.

I have implemented two models sde_ui3 and sde_ui10, reflecting a SIS model with virus mutation (the first with 3 variants and 4 immunity states, the second with 10 variants and...