Sam Isaacson

Algorithm 3.4.1 gives pseudo-code.

@Vilin97 will merge when I get tests passing. @TorkelE expanded this a bit to allow not needing to pass `SSAStepper` too, i.e. with this PR this should now work ```julia...

@TorkelE actually MTK and Catalyst need appropriate dispatches for JumpProblem to make this work, so not quite there yet for being usable in Catalyst.

@ChrisRackauckas if you could give a quick look to ensure I haven't added a dispatch better handled elsewhere (or shadowing one elsewhere) that would be helpful.

@Vilin97 thanks for your work on this and sorry for the long long delay!

Yes, the docs need updates in places now. Using the defaults in the most introductory places, but explaining how to select solvers and such as they currently do in more...

We should probably document this more clearly though (i.e. we tell people in the FAQs to use `remake` but don't explain not to directly mutate).

Well, see the Slack you were just posting on.

@Vilin97 has an open PR that I haven’t had a chance to review yet which addresses default aggregators.

One related thing that needs to be added at the MTK level is setting up a simpler workflow when there are VariableRateJumps.