def METHOD_NAME():
    A = Square('A',bounds1=(0, 0.5), bounds2=(0, 1))
    B = Square('B',bounds1=(0, 0.5), bounds2=(0, 1))
    M1 = AffineMapping('M1',2, c1=0.5, c2=0.,a11=1,  a22=1, a12=0, a21=0)
    M2 = AffineMapping('M2',2, c1=0.5, c2=0, a11=-1, a12=0, a21=0, a22=1)
    D1 = M1(A)
    D2 = M2(B)
    connectivity = [((0,0,-1),(1,0,-1))]
    patches = [D1,D2]
    domain = Domain.join(patches, connectivity, 'domain')
    x,y = domain.coordinates
    solution = sin(pi*x)*sin(pi*y)
    f        = 2*pi**2*solution
    l2_error, h1_error = run_poisson_2d(solution, f, domain, ncells=[2**2,2**2], degree=[2,2])
    expected_l2_error = 0.0020352296663934746
    expected_h1_error = 0.05679638799164659
    assert ( abs(l2_error - expected_l2_error) < 1e-7 )
    assert ( abs(h1_error - expected_h1_error) < 1e-7 )