diff --git a/Geo/GModelIO_NEU.cpp b/Geo/GModelIO_NEU.cpp
index f937ae4340c2967b758fd9baeb593c6ce5b200dc..022b5db4e662224b17bbb6dfa25436aa0774e21b 100644
--- a/Geo/GModelIO_NEU.cpp
+++ b/Geo/GModelIO_NEU.cpp
@@ -117,6 +117,11 @@ namespace {
       }
     }
 
+    // because we use a set_intersection later on, the vectors of vertmap should be sorted
+    for (auto &it: vertmap) {
+      std::sort(it.second.begin(), it.second.end());
+    }
+
     // determine which faces belong to which tetrahedra by comparing vertices
     IDTetFaceMap tetfacemap;
     for(GModel::riter it = gm->firstRegion(); it != gm->lastRegion(); ++it) {