def METHOD_NAME():
    m = UnitIntervalMesh(1)
    mesh = ExtrudedMesh(m, layers=2)
    U = VectorFunctionSpace(mesh, 'DG', 1)
    v = TestFunction(U)
    n = FacetNormal(mesh)
    temp = assemble(jump(conj(v), n)*dS_h).dat.data
    assert np.all(temp[:, 0] == 0.0)
    assert not np.all(temp[:, 1] == 0.0)