def METHOD_NAME():
    mesh = UnitSquareMesh(1, 2, quadrilateral=True)
    U = VectorFunctionSpace(mesh, 'DG', 0)
    u = TrialFunction(U)
    v = TestFunction(U)
    n = FacetNormal(mesh)
    temp = assemble(avg(inner(dot(u, n), dot(v, n)))*dS)
    assert temp.M.values[0, 0] == 0.0
    assert temp.M.values[1, 1] != 0.0
    assert temp.M.values[2, 2] == 0.0
    assert temp.M.values[3, 3] != 0.0