Frequently updated Github repositories
Many of the codes below, and newer ones, can now be found in Github repositories that can be accessed via the following link:
https://github.com/acerfon
Exact axisymmetric plasma equilibria with Solov'ev profiles
The short MATLAB codes below compute exact toroidally axisymmetric equilibria with Solov'ev profiles as described in A.J. Cerfon and J.P. Freidberg,
"One size fits all" analytic solutions to the Grad-Shafranov equation,
Physics of Plasmas 17, 032502 (2010)
ECOM: Fast and high order accurate Grad-Shafranov solver
This is a FORTRAN code described in detail in the following articles:
A.Pataki, A.J. Cerfon, J.P. Freidberg, L.Greengard, M. O'Neil,
A fast, high-order solver for the Grad-Shafranov equation,
Journal of Computational Physics 243, 28 (2013)
J.P. Lee and A.J. Cerfon,
ECOM: a fast and accurate solver for toroidal axisymmetric MHD equilibria,
Computer Physics Communications 190, 72 (2015).
It is freely available in the following SVN repository:
https://sourceforge.net/projects/ecom-gssolver/
Fast free space Poisson solvers
The short MATLAB codes below compute Poisson's equation with free space boundary conditions with spectral accuracy, only using FFTs, as described in F. Vico, L. Greengard, and M. Ferrando,
Fast convolution with free-space Green's functions,
Journal of Computational Physics 323, 191 (2016)
Exact axisymmetric Taylor states for shaped plasmas
The short MATLAB codes below compute exact toroidally axisymmetric Taylor states profiles as described in A.J. Cerfon and M. O'Neil,
Exact axisymmetric Taylor states for shaped plasmas,
Physics of Plasmas 21, 064501 (2014)
Single particle stability and trajectories in a linear Paul trap