Open navidcy opened 4 years ago
We can implement it with every stepper, though as you allude to changing the time step in ETDRK4
would require recomputing coefficients. I suppose what we need is to design a flexible abstraction that works for many types of problems.
Have you considered OrdinaryDiffEq.jl
for time-stepping? They have many solvers with adaptivity, and some with GPU-compatibility, though I haven't kept myself up-to-date recently.
It would be very nice if we could do that. How feasible is it? Or course it would, in practice, only work with the explicit one time-steppers (
ForwardEuler
,AB3
andRK4
)