Guiding Center 4d Barely Passing with Symmetric SRK3 VPRK Method and Projections
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.
SymplecticSRK3
Simulation crashed after 4304 of 12500 time steps: solver error – non-finite direction vector.
Plotting vprksrk3pnonetoroidalmomentum_error failed: MethodError.
Plotting vprksrk3pnonetoroidalmomentum_drift failed: MethodError.
SymplecticSRK3 with standard projection
Plotting vprksrk3pstandardtoroidalmomentum_error failed: MethodError.
Plotting vprksrk3pstandardtoroidalmomentum_drift failed: MethodError.
SymplecticSRK3 with symmetric projection
Plotting vprksrk3psymmetrictoroidalmomentum_error failed: MethodError.
Plotting vprksrk3psymmetrictoroidalmomentum_drift failed: MethodError.
SymplecticSRK3 with symplectic projection
Plotting vprksrk3psymplectictoroidalmomentum_error failed: MethodError.
Plotting vprksrk3psymplectictoroidalmomentum_drift failed: MethodError.
SymplecticSRK3 with midpoint projection
Plotting vprksrk3pmidpointtoroidalmomentum_error failed: MethodError.
Plotting vprksrk3pmidpointtoroidalmomentum_drift failed: MethodError.
SymplecticSRK3 with internal stage projection
Simulation crashed after 0 of 12500 time steps: MethodError.