Implement GVEC interface `template_gvec` and a new `kind_map` "cuboid slice".
Commits:
- Plot iota and q profiles as absolute values.
- Add
equilibrium_MHD.py
interface intemplate_gvec
. - Move diagnostics
GVEC_MHD_equilibrium.py
to teststest_GVEC_MHD_equilibrium.py
. - Modify
test_GVEC_MHD_equilibrium.py
into atest_
function. - Remove unused 1D profiles, and rename interface for toroidal and poloidal flux in
template_gvec
. - Implement new
kind_map
"cuboid slice". - Basic support of
source_domain
in GVEC'sequilibrium_MHD
interface. - Support
source_domain
transforms of geometric terms. - Rename all "profile" to "variable" in
gvec_to_python
package. - Provide derivatives of toroidal and poloidal fluxes in
template_gvec
. - Update
test_GVEC_MHD_equilibrium.py
. - Fix matrix multiplication in
template_gvec
'sequilibrium_MHD
. - Update GitLab CI's image from python:latest to python:3.9-bullseye.
- Create test for
template_gvec
. - Reference updated GVEC equilibrium data in code.
- Fix pyccel not recognizing tuple unpacking.
TODO:
- Fix blowing up MHD variables near
s=0
due to derivatives of splines (dphi/ds, dchi/ds) not clamped to zero, by (a) using analytic toroidal flux, or (b) read the derivatives directly as profiles, or (c) add clamp condition to solve for a new set of DoF. - Fix equilibrium current in
template_gvec
'sequilibrium_MHD
interface, after the above is fixed.