def METHOD_NAME(p):
    node = np.array([
            (0.0, 0.0), (0.0, 1.0), (0.0, 2.0), (1.0, 0.0), (1.0, 1.0), (1.0, 2.0),
            (2.0, 0.0), (2.0, 1.0), (2.0, 2.0)], dtype=np.float64)
    cell = np.array([0, 3, 4, 4, 1, 0, 1, 4, 5, 2, 3, 6, 7, 4, 4, 7, 8, 5],
                dtype=np.int_)
    cellLocation = np.array([0, 3, 6, 10, 14, 18], dtype=np.int_)
    mesh = PolygonMesh(node, cell, cellLocation)
    space =  ConformingVectorVESpace2d(mesh, p=p)
    smldof = space.smspace.number_of_local_dofs()
    m = ScaledMonomialSpaceMassIntegrator2d()
    M = m.assembly_cell_matrix(space.smspace, p=p+1) 
    H1projector = ConformingVectorVEMH1Projector2d(M[:,:smldof,:smldof])
    PI1 = H1projector.assembly_cell_matrix(space)    
    L2projector = ConformingVectorVEML2Projector2d(M, PI1) 
    H = L2projector.assembly_cell_left_hand_side(space)