Antoine Cerfon

Dear all, I do not think symplectic integration would lead to faster field line tracing. For a given accuracy, a non-symplectic ODE solver will be faster than a symplectic one...

Hi Chris, I fully agree with your single turn/interpolation idea. I also agree that RK4/5 is pretty bad. But then again, do not discard Spectral Deferred Correction. The best one...

For what it's worth, we have been very happy with the Flatiron's NUFFT, which we have used with much success for a new Particle in Fourier scheme: https://www.sciencedirect.com/science/article/pii/S0021999124006387#:~:text=Journal%20of%20Computational%20Physics.%20Volume I highly...

As another 2-cents, for the singular integrals you are talking about, I very highly recommend the Kapur-Rokhlin quadrature schemes, which we have also found easy to use, and super accurate,...