def METHOD_NAME(space, mesh):
    """
    Verify that the 1-form Hodge Laplacian with strong Dirichlet
    boundary conditions has kernel of dimension equal to the 1st Betti
    number of the annulus mesh, i.e. 1.
    """
    V0tag, V1tag, V2tag = space
    if len(V0tag) == 2:
        V0 = FunctionSpace(mesh, V0tag[0], V0tag[1])
    else:
        V0a = FiniteElement(V0tag[0], "triangle", V0tag[1])
        V0b = FiniteElement(V0tag[2], "triangle", V0tag[3])
        V0 = FunctionSpace(mesh, V0a + V0b)
    V1 = FunctionSpace(mesh, V1tag[0], V1tag[1])
    W = V0*V1
    sigma, u = TrialFunctions(W)
    tau, v = TestFunctions(W)
    L = assemble((inner(sigma, tau) - inner(u, rot(tau)) + inner(rot(sigma), v)
                  + inner(div(u), div(v))) * dx)
    bc0 = DirichletBC(W.sub(0), Constant(0.0), 9)
    bc1 = DirichletBC(W.sub(1), Constant((0.0, 0.0)), 9)
    L0 = assemble((inner(sigma, tau) - inner(u, rot(tau)) + inner(rot(sigma), v)
                   + inner(div(u), div(v))) * dx, bcs=[bc0, bc1])
    dV0 = V0.dof_count
    dV1 = V1.dof_count
    A = numpy.zeros((dV0+dV1, dV0+dV1), dtype=utils.ScalarType)
    A[:dV0, :dV0] = L.M[0, 0].values
    A[:dV0, dV0:dV0+dV1] = L.M[0, 1].values
    A[dV0:dV0+dV1, :dV0] = L.M[1, 0].values
    A[dV0:dV0+dV1, dV0:dV0+dV1] = L.M[1, 1].values
    u, s, v = linalg.svd(A)
    nharmonic = sum(s < 1.0e-5)
    assert nharmonic == 1
    dV0 = V0.dof_count
    dV1 = V1.dof_count
    A0 = numpy.zeros((dV0+dV1, dV0+dV1), dtype=utils.ScalarType)
    A0[:dV0, :dV0] = L0.M[0, 0].values
    A0[:dV0, dV0:dV0+dV1] = L0.M[0, 1].values
    A0[dV0:dV0+dV1, :dV0] = L0.M[1, 0].values
    A0[dV0:dV0+dV1, dV0:dV0+dV1] = L0.M[1, 1].values
    u, s, v = linalg.svd(A0)
    nharmonic = sum(s < 1.0e-5)
    assert nharmonic == 1