Integrators - minimal working example for higher order integration methods

[dom-kufel]

Hi, I am interested in a minimal working example for using higher order integration methods. I cannot immediately see how does their implementation work. For instance, given a TDVP linear ODE in custom_simple_tdvp() from https://scipost.org/SciPostPhysCodeb.7/pdf how would you modify simple Euler rule to e.g., RK12 while utilizing methods from netket.experimental.dynamics?