Christopher Rackauckas
Christopher Rackauckas
I think that could make sense.
Yup this is just a missing part of the interface.
Related https://github.com/SciML/NonlinearSolve.jl/pull/206 @frankschae
A good case to test this on would be the double pendulum https://sciml.github.io/ModelingToolkitCourse/dev/lectures/lecture8/#Bonus-Demo
> The policy is to never mutate these directly though, and always provide functions for this. Not that there is anything that can be done to stop people, but in...
``` help?> ODEProblem search: ODEProblem RODEProblem SplitODEProblem DynamicalODEProblem ODEAdjointProblem IncrementingODEProblem RODEAdjointProblem Defines an ordinary differential equation (ODE) problem. Documentation Page: https://docs.sciml.ai/DiffEqDocs/stable/types/ode_types/ Mathematical Specification of an ODE Problem ============================================ To define...
> I think the reason was that (long time ago) the idea was to neatly combine Gillespie simulations with ODEs? I cannot speak for everyone, but my impression is that...
> Finally, I have never used anything but SSAStepper() in `solve. Is there a reason that this is not a default argument, given how common it is? That should be...
> (This is something that comes up whenever I show people how to do systems bio in Julia. "This is how to do ODE", "This is how to do SDE",...
> So, I agree that having hybrid simulations being the default should be the norm and possible (and have a good interface). I still think for tracking people about CRNs...