don't crash when trying to do HO meshes of componds - just do linear interp

CMakeLists.txt

@@ -94,7 +94,7 @@ opt(ZIPPER "Enable Zip file compression/decompression" OFF)
 
 set(GMSH_MAJOR_VERSION 2)
 set(GMSH_MINOR_VERSION 13)
-set(GMSH_PATCH_VERSION 0)
+set(GMSH_PATCH_VERSION 1)
 set(GMSH_EXTRA_VERSION "" CACHE STRING "Gmsh extra version string")
 
 set(GMSH_VERSION "${GMSH_MAJOR_VERSION}.${GMSH_MINOR_VERSION}")

Mesh/HighOrder.cpp

@@ -243,7 +243,8 @@ static void getEdgeVertices(GEdge *ge, MElement *ele, std::vector<MVertex*> &ve,
                             bool linear, int nPts = 1)
 {
   if(ge->geomType() == GEntity::DiscreteCurve ||
-     ge->geomType() == GEntity::BoundaryLayerCurve)
+     ge->geomType() == GEntity::BoundaryLayerCurve ||
+     ge->geomType() == GEntity::CompoundCurve)
     linear = true;
 
   std::vector<MVertex*> veOld;
@@ -272,7 +273,8 @@ static void getEdgeVertices(GFace *gf, MElement *ele, std::vector<MVertex*> &ve,
                             bool linear, int nPts = 1)
 {
   if(gf->geomType() == GEntity::DiscreteSurface ||
-     gf->geomType() == GEntity::BoundaryLayerSurface)
+     gf->geomType() == GEntity::BoundaryLayerSurface ||
+     gf->geomType() == GEntity::CompoundSurface)
     linear = true;
 
   for(int i = 0; i < ele->getNumEdges(); i++) {
@@ -598,7 +600,8 @@ static void getFaceVertices(GFace *gf, MElement *incomplete, MElement *ele,
                             faceContainer &faceVertices, bool linear, int nPts = 1)
 {
   if(gf->geomType() == GEntity::DiscreteSurface ||
-     gf->geomType() == GEntity::BoundaryLayerSurface)
+     gf->geomType() == GEntity::BoundaryLayerSurface ||
+     gf->geomType() == GEntity::CompoundSurface)
     linear = true;
 
   MFace face = ele->getFace(0);