Guiding Center 4d Barely Passing with Radau IIA VPRK Methods

The time step and number of time steps are those of this orbit, set in src/guiding-center-4d.jl; the toroidal momentum is shown alongside the energy, as the second invariant of the guiding centre dynamics.

VPRKRadauIIA(2)

Plots

Plotting vprkradauIIA2toroidalmomentum_error failed: MethodError.

Plotting vprkradauIIA2toroidalmomentum_drift failed: MethodError.

VPRKRadauIIA(3)

Plots

Plotting vprkradauIIA3toroidalmomentum_error failed: MethodError.

Plotting vprkradauIIA3toroidalmomentum_drift failed: MethodError.