This release introduces
1) The ability to compute fieldlines as well as particle trajectories using the guiding center approximation.
**Credits to Florian Wechsung.**
2) A quadratic-flux minimizing surface solver, defined as the minimizer of the objective function
f = int d^2 x (B \cdot n)^2/int d^2 x B^2
at fixed flux surface label (e.g., volume, area, or toroidal flux). This optimization problem can be solved using a
penalty formulation with LBFGS or with the SLSQP constrained optimization algorithm. Much of the code was
modeled after the Boozer surfaces.
**Credits to Elizabeth J Paul**