def METHOD_NAME():
    A = Square('A',bounds1=(0, 0.5), bounds2=(0, 1))
    B = Square('B',bounds1=(0.5, 1.), bounds2=(0, 1))
    connectivity = [((0,0,1),(1,0,-1))]
    patches = [A,B]
    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.002035229666394183
    expected_h1_error = 0.056796387991647795
    assert ( abs(l2_error - expected_l2_error) < 1e-7 )
    assert ( abs(h1_error - expected_h1_error) < 1e-7 )