Releases: hiddenSymmetries/simsopt
Graph Optimization framework for all of simsopt classes
This release extends the graph Optimizable framework to surface and MHD classes.
Changes Implemented:
- Enhancements to graph Optimizable class to keep order of the dependent functions
- ObjectiveFailure added to LeastSquaresProblem
- Tolerances can be specified as part of the solver
- Surface classes and the corresponding tests are separated into different modules
- Surface modules converted to graph framework
- MHD modules converted to graph framework
- New classes for finite differences based (both serial and MPI) Jacobian
- Weak references to enable garbage collection of out-of-scope optimizable objects
- SPEC CI issues fixed
New Coil-n-Currents, Derivatives, and stage II coil optimization
Courtesy @florianwechsung, this release introduces new Coil, Current, and Derivative classes and allows for stage-II optimization
-
A new
Coil
class that represents aCurve
and aCurrent
. A list ofCoil
objects is now passed toBiotSavart
(instead of separatecurves
andcurrents
in the past). These objects are built on top of the new graph framework -
A new
Derivative
object is introduced, that is essentially a dictionary containing the partial derivative of a scalar objective wrt to all the Optimizable objects (with >0 DOFs) that it depends on. -
Simsopt can now solve the Stage II coil optimisation problem.
Addition of adjoint methods, vmec residuals and improved plotting
This release adds adjoint methods (courtesy E J Paul), additional VMEC related residuals (courtesy M Landreman) and improved plotting options (courtesy C. Zhu).
v0.4.2 Released
This release introduces
-
The ability to compute fieldlines as well as particle trajectories using the guiding center approximation.
Credits to Florian Wechsung. -
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
v0.4.0 Released
This release adds conda package recipe for simsopt, toroidal surface using the parameterization in Henneberg, Helander, and Drevlak ( arXiv:2105.00768 (2021), and bug fixes.
v0.3.03 Released
This release fixes bugs and reorganizes examples.
- Earlier releases of simsopt required all optional packages to be installed. This release fixes issues in importing simsopt when optional packages are not installed.
- The examples are reorganized into simple, intermediate, advanced, and stellarator_benchmarks.
Docker and coverage features added to simsopt
This release eases the use of simsopt for end users. A docker container with simsopt and vmec preinstalled is available. It can be downloaded from DockerHub using the name hiddensymmetries/simsopt.
The coverage report is now avaialble and can be accessed from the codecov badge on the main README.me
v0.2.02
simsopt v0.2.0
This release introduces some additions and bug fixes for majority of modules in simsopt.
Changes introduced:
- New class for Reiman magnetic field (see section 5 of Reiman, Greenside (1986), "Calculation of three-dimensional MHD equilibria with islands and stochastic regions", Computer Physics Communications 43 (157-167)) is added by Rogerio. This field has the advantage of allowing an exact calculation of the island width which can be used for island width optimization (example: https://arxiv.org/abs/2102.04497).
- Faster Biot-Savart code by Florian
- A simplified logging interface is introduced by Bharat. Logging can be enabled by using the following two lines in your driver script:
from simsopt import initialize_logging
initialize_logging(filename="simsopt.log") # Filename option can be omitted in which case the default name is default.log
- Some modules are rearranged.
- Top-level imports for important classes are introduced.
- Automatic linting script implemented by Florian. Before you submit your additions/changes to simsopt as PR, please run the
run_autopep
script. Requires autopep8 and flake8. - A new exception class to handle failures in objective function (credits to Matt Landremann).
- Bug fixes.
- Improved documentation.
Boozer surfaces in simsopt
This release introduces boozer surface and computation of aspect ratio of surfaces and some major refactoring of surfaces-related code.
Major changes:
-
SurfaceXYZFourier surfaces are introduced, where (phi, theta) correspond to the Boozer angles on the surface, by solving a nonlinear least squares problem in boozersurface.py.
-
minimizeBoozerScalarizedLBFGS solves a scalarized constrained optimization problem using LBFGS, where
min 0.5* || f(x) ||^2_2 + 0.5 * constraint_weight * (label - labeltarget)^2
+0.5constraint_weight * (y(varphi=0,theta=0) -0)^2+0.5constraint_weight * (z(varphi=0,theta=0) -0)^2
where || f(x)||^2_2 is the sum of squares of the Boozer residual at
a set of quadrature points.
The final two terms in the objective function are boundary conditions to ensure that the parametrization of the surface is unique. -
minimizeBoozerScalarizedNewton solves the same scalarized constrained optimization problem as above, but instead of LBFGS, Newton's method is used.
-
minimizeBoozerConstrainedNewton solves the full constrained optimization problem using Newton's method.
min 0.5*|| f(x) ||^2_2
subject to
label - labeltarget = 0
y(varphi=0,theta=0) - 0 = 0
z(varphi=0,theta=0) - 0 = 0
The user is free to choose either a geometric constraint, e.g. area, or a toroidal flux constraint.
objectives.py is split into coilobjectives.py and surfaceobjectives.py -
Computation of aspect ratio of surfaces using the VMEC definition of aspect ratio is done in the Surface.aspect_ratio.
-
New capability of computing cross-sections of surfaces at a given cylindrical angle. This is done using bisection in the Surface.cross_section.
-
An implementation of SurfaceXYZFourier.to_RZFourier using Surface.cross_section. This is done by computing cross-sections of the SurfaceXYZFourier object and then fitting a SurfaceRZFourier to those cross-sections.