Could someone point me in the direction of where to find info on how the time-stepping method for stiff SDE solvers is implemented? I am mainly interested in how the ImplicitEM solver does the adaptive stepsize control, which is the scheme I use. The documentation says that it “Uses a 1.0/1.5 heuristic for adaptive time stepping” for the order 0.5 and 1.0 schemes. I guess this is related to the paper: 10.3934/dcdsb.2017133 by Rackauckas et al., which as a physicist not specialized in SDE solvers, cannot see how it can be applied to the ImplicitEM scheme.
I also want to thank the SciML organization and community for all the amazing tools created!!