stochastic_solvers calls transition_matrix without passing M_
transition_matrix optionally takes M_ as an additional argument. If it doesn't get the extra argument, then it uses the global M_. Given that stochastic_solvers takes M_ as an additional argument, it should also pass this to transition_matrix. Without this, user code that calls Dynare internals in parallel can be messed up.