Constraint degeneracy

[Alan-Goodman]

When I have this in my optimisation:

R_modes = np.vstack(
        (
            [0, 0, 0],
            eq.surface.R_basis.modes[
                np.max(np.abs(eq.surface.R_basis.modes), 1) > k, :
            ],
        )
    )
    Z_modes = eq.surface.Z_basis.modes[
        np.max(np.abs(eq.surface.Z_basis.modes), 1) > k, :
    ]
    # next we create the constraints, using the mode number arrays just created
    # if we didn't pass those in, it would fix all the modes (like for the profiles)
    constraints = (
        ForceBalance(eq=eq),
        # try to keep the aspect ratio between 7 and 9
        AspectRatio(eq=eq, bounds=(0, 10)),
        # similarly, try to keep it from getting too elongated
        Elongation(eq=eq, bounds=(0, 6)),
        # Keep volume the same as the initial volume
        #Volume(eq=eq0, target=eq0.compute("V")["V"]),
        # target for average iota
        RotationalTransform(eq=eq, target=1.1, loss_function="mean"),
        # fix major radius
        FixBoundaryR(eq=eq, modes=[0, 0, 0]),
        FixBoundaryR(eq=eq, modes=R_modes),
        # fix vacuum profiles
        FixPressure(eq=eq),
        FixCurrent(eq=eq),
        FixPsi(eq=eq),
    )
    optimizer = Optimizer("lsq-auglag")

...there is a bug, because of FixBoundaryR(0,0,0) and the FixBoundaryR(modes), which both fix the (000) mode. This should presumably not bug out, I assume?

[dpanici]

Correct, thanks for catching this Alan!