As with the STELLOPT and ROSE codes, the ALPOpthttps://github.com/ejpaul/ALPOpt tool optimizes the boundary of VMEC [1]} equilibria. The VMEC code obtains solutions of () under the assumption of nested toroidal magnetic surfaces. The plasma boundary is described by a set of Fourier coefficients of the cylindrical coordinates, \(\lbrace R_{m,n}^c,Z_{m,n}^s\rbrace \) ,
\(R(\theta ,\phi ) &= \sum _{m,n} R_{m,n}^c \cos (m \theta - n N_P \phi ) \\Z(\theta ,\phi ) &= \sum _{m,n} Z_{m,n}^s \sin (m \theta - nN_P \phi ),\) 
