Old LEoPart did not have callback for the velocity function on intermediate RK steps. Consider, for example, the implementation here wrapped for the python layer here.
This facilitates optimal rates of convergence of the RK method for time dependent velocity fields.
Old LEoPart did not have callback for the velocity function on intermediate RK steps. Consider, for example, the implementation here wrapped for the python layer here.
This facilitates optimal rates of convergence of the RK method for time dependent velocity fields.