def METHOD_NAME():
    mesh = UnitSquareMesh(2, 1, quadrilateral=True)
    U = VectorFunctionSpace(mesh, 'DG', 0)
    V = FunctionSpace(mesh, 'RTCF', 1)
    W = U*V
    u1, u2 = TrialFunctions(W)
    v1, v2 = TestFunctions(W)
    pp = assemble(inner(u2('+'), v1('+'))*dS)
    pm = assemble(inner(u2('+'), v1('-'))*dS)
    mp = assemble(inner(u2('-'), v1('+'))*dS)
    mm = assemble(inner(u2('-'), v1('-'))*dS)
    assert not np.all(pp.M[0, 1].values == pm.M[0, 1].values)
    assert not np.all(pp.M[0, 1].values == mp.M[0, 1].values)
    assert not np.all(pp.M[0, 1].values == mm.M[0, 1].values)
    assert not np.all(pm.M[0, 1].values == mp.M[0, 1].values)
    assert not np.all(pm.M[0, 1].values == mm.M[0, 1].values)
    assert not np.all(mp.M[0, 1].values == mm.M[0, 1].values)