Merge pull request #111 from wenqing/fixing1D_DM

CFEMesh::GenerateHighOrderNodes(): Fixed a bug when 1D elements mixed with 3D elements
ufz · Sep 7, 2018 · e6bf66d · e6bf66d
2 parents feff81b + f91a0f8
commit e6bf66d
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/MSH/msh_mesh.cpp b/MSH/msh_mesh.cpp
@@ -979,6 +979,10 @@ void CFEMesh::GenerateHighOrderNodes()
 					if (ee == e)
 						continue;
 					thisElem = ele_vector[ee];
+
+					if (thisElem->GetElementType() == MshElemType::LINE)
+						continue;
+
 					nedges = thisElem->GetEdgesNumber();
 					// Edges of neighbors
 					for (ii = 0; ii < nedges; ii++)
@@ -1065,6 +1069,9 @@ void CFEMesh::GenerateHighOrderNodes()
 			for (size_t i = 0; i < thisElem0->GetFacesNumber(); i++)
 			{
 				thisElem = thisElem0->GetNeighbor(i);
+				if (!thisElem)
+					continue;
+
 				// look for adjacent solid elements
 				if (thisElem->GetElementType() == MshElemType::LINE)
 					continue;