Strang Splitting with SplitODEProblem

Hi all,

I am interested in implementing a split ODE problem,
I am solving the wave equation with a Finite Volume scheme such as MUSCL with a discontinuous source term.

I want to implement a fractional step method and SciML’s SplitODEProblem seems to be suitable. But as I see most methods are only first order in the split. The second order split would be Strang’s split, which is symmetric.

How could I implement that or something higher order?

I don’t think we have that currently implemented and IIRC we have it open as a GSoC project :sweat_smile:

ah ok! thank you. I am working on something towards this, if it works I will do a pr so that you can check it out.