From 3f1efc80dd6d583e4359c161c7980a364999a776 Mon Sep 17 00:00:00 2001
From: Oliver Gloth <ogloth@engits.com>
Date: Tue, 12 Jul 2011 14:57:59 +0200
Subject: [PATCH] added error message if surface is not water-tight

---
 src/libengrid/insertpoints.cpp | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/libengrid/insertpoints.cpp b/src/libengrid/insertpoints.cpp
index 11d1104..55c43f2 100644
--- a/src/libengrid/insertpoints.cpp
+++ b/src/libengrid/insertpoints.cpp
@@ -216,6 +216,10 @@ char InsertPoints::getNewNodeType(stencil_t S)
   
   vtkIdType id_node1 = S.p1;
   vtkIdType id_node2 = S.p2;
+
+  if (S.id_cell.size() != 2) {
+    EG_ERR_RETURN("The surface mesh is not water-tight");
+  }
   
   EG_VTKDCN(vtkCharArray, node_type, m_Grid, "node_type");
   if( node_type->GetValue(id_node1)==VTK_SIMPLE_VERTEX || node_type->GetValue(id_node2)==VTK_SIMPLE_VERTEX ) {
-- 
2.11.4.GIT