"""Fixed-boundary Grad-Shafranov on MAST: FEM Q2 and PINN, same model, compared.

The equilibrium of the EQDSK file, recomputed from its profiles FF' and P'
with the wall flux as Dirichlet datum -- the PINN example
``pinns/stationary_pdes/elliptic_pdes/grad_shafranov_with_xpoints.py mast``
(same physics, boundary weight and network retuned, see ``_fem_vs_pinn.py``), next to the same model solved by finite elements
(:class:`~scimba_jax.physical_models.elliptic_pde.grad_shafranov_fem.
GradShafranovFullPhysicsWeakModel`). Both recompute the magnetic axis and the
X-points of the current iterate before every step
(``pre_computation_without_diff``); the FEM finds them cell by cell and
iterates Picard to ``||F|| < 1e-10``.

Prints the axis, X-points and boundary flux of both, the gaps FEM-PINN,
FEM-EQDSK and PINN-EQDSK at the FEM nodes, and plots the three fluxes and the
three gaps. The setup is in ``_fem_vs_pinn.py``. Needs the optional ``mesh``
extra (gmsh); the mesh is cached next to this file.
"""

from _fem_vs_pinn import run

run("mast")