def METHOD_NAME():
    mesh = UnitSquareMesh(2, 1, quadrilateral=True)
    U = VectorFunctionSpace(mesh, 'DG', 1)
    v = TestFunction(U)
    n = FacetNormal(mesh)
    temp = assemble(jump(conj(v), n)*dS).dat.data
    assert not np.all(temp[:, 0] == 0.0)
    assert np.all(temp[:, 1] == 0.0)