There is a problem with buildmesh when doing domain decomposition: by default Gfem swap the diagonals at the corners of the domain if the triangle has two boundary edges. This will lead to bad domain decomposition at the sub-domain interfaces.
To solve this, there is a new flag for buildmesh which is optional:
buildmesh(<max_number_of_vertices>, <flag>)
where <flag> =