Create General RK Nystrom Stepper from Fine's 3(4) method

jonapost / field_propagation

Module for integrating a track's trajectory in a field, whether magnetic, electric, combined electromagnetic, or also including gravity or other forces.

6 stars 20 forks source link

Create General RK Nystrom Stepper from Fine's 3(4) method #1

Open jonapost opened 9 years ago

jonapost commented 9 years ago

Define an RKNG stepper that follows Fine's lower order method. ( Please add reference here. )

jsuagee commented 9 years ago

Wrote FineRKNG34 method with interpolation (From J.M. Fine's "Low Order Practical Runge-Kutta-Nystrom Methods," Computing 38, 1987). Interpolation is currently only for position, not velocity, but this is mostly because I haven't entered the coefficients yet for the first derivative interpolation. The most up to date version is in branch chawla-accepted-step-experimental.

jonapost commented 9 years ago

Good work!

Do you think that you can find the coefficients for the velocity in electronic form, so that you don't have to type them ?