IncrementalMesh crash with Common generated surface.

I have two boxes touching at a part of the boundary. I can compute the intersection (BRepAlgoAPI_Common) but when I try to mesh the surface I get a crash. Perhaps some sees what I am doing wrong? Does the intersection perhaps return both intersecting faces from both boxed and I need to merge them somehow? Thanks.

Here is the code to reproduce the issue:

#include
#include
#include
#include

int main()
{
gp_Pnt gp1 = gp_Pnt(0,0,0);
gp_Pnt gp2 = gp_Pnt(1,1,1);
gp_Pnt gp3 = gp_Pnt(1,0.5,0.5);
gp_Pnt gp4 = gp_Pnt(2,1.5,1.5);
TopoDS_Shape s1 = BRepPrimAPI_MakeBox(gp1, gp2).Shape();
TopoDS_Shape s2 = BRepPrimAPI_MakeBox(gp3, gp4).Shape();
TopoDS_Shape shape = BRepAlgoAPI_Common(s1, s2).Shape();
BRepMesh_IncrementalMesh Mesh( shape, 0.01 );
Mesh.Perform();
return 0;
}