As discussed with @MichelJuillard, we should try to get rid of globals in - `stoch_simul` - `simult_` as most of their callers and subsequently called functions already use local instances of the variables