pub fn vmec_fixed_boundary_solve(
boundary: &VmecBoundaryState,
config: &VmecSolverConfig,
pressure: &[f64],
iota: &[f64],
phi_edge: f64,
) -> FusionResult<VmecEquilibrium>Expand description
Solve VMEC fixed-boundary 3D equilibrium.
Given boundary shape (from VmecBoundaryState), pressure and rotational
transform profiles, finds force-balanced interior flux surface shapes.
The algorithm minimises the MHD energy functional W = ∫(B²/2μ₀ + p)dV via steepest descent on the interior Fourier coefficients R_mn^c(s), Z_mn^s(s) while holding the axis and boundary fixed.