def METHOD_NAME(self):
    """Verify the available definitions."""
    surf = SurfaceRZFourier.from_vmec_input(filename)
    ntheta = len(surf.quadpoints_theta)
    nphi = len(surf.quadpoints_phi)
    ncoils = 3
    base_curves = create_equally_spaced_curves(
        ncoils, surf.nfp, stellsym=surf.stellsym, R0=1.0, R1=0.5, order=6
    )
    base_currents = [Current(1e5) for i in range(ncoils)]
    coils = coils_via_symmetries(base_curves, base_currents, surf.nfp, surf.stellsym)
    bs = BiotSavart(coils)
    # Test definition = "quadratic flux":
    target = np.ones(surf.gamma().shape[0:2])
    J = SquaredFlux(surf, bs, target, definition="quadratic flux").J()
    bs.set_points(surf.gamma().reshape((-1, 3)))
    B = bs.B()
    normal = surf.normal().reshape((-1, 3))
    norm_normal = np.sqrt(normal[:, 0]**2 + normal[:, 1]**2 + normal[:, 2]**2)
    B_dot_n = np.sum(B * surf.unitnormal().reshape((-1, 3)), axis=1)
    should_be = 0.5 * sum((B_dot_n - target.reshape((-1,)))**2 * norm_normal) / (ntheta * nphi)
    np.testing.assert_allclose(J, should_be)
    # Test definition = "normalized":
    J2 = SquaredFlux(surf, bs, target, definition="normalized").J()
    mod_B_squared = np.sum(B * B, axis=1)
    numerator = 0.5 * sum(
        (B_dot_n - target.reshape((-1,)))**2 * norm_normal
    ) / (ntheta * nphi)
    denominator = sum(mod_B_squared * norm_normal) / (ntheta * nphi)
    np.testing.assert_allclose(J2, numerator / denominator)
    # Test definition = "local":
    J3 = SquaredFlux(surf, bs, target, definition="local").J()
    should_be3 = 0.5 * sum(
        (B_dot_n - target.reshape((-1,)))**2 / mod_B_squared * norm_normal
    ) / (ntheta * nphi)
    np.testing.assert_allclose(J3, should_be3)
    with self.assertRaises(ValueError):
        J4 = SquaredFlux(surf, bs, target, definition="foobar")