Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 26997)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 26998)
@@ -3197,5 +3197,5 @@
 	norm=sqrt(normal[0]*normal[0]+normal[1]*normal[1]+normal[2]*normal[2]);
 
-	for(int i=0;i<3;i++) normal[i]=normal[i]/norm;
+	for(int i=0;i<3;i++) normal[i]=normal[i]/(norm+1e-10);
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp	(revision 26997)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp	(revision 26998)
@@ -605,5 +605,5 @@
 	norm=sqrt(normal[0]*normal[0]+normal[1]*normal[1]+normal[2]*normal[2]);
 
-	for(int i=0;i<3;i++) normal[i]=normal[i]/norm;
+	for(int i=0;i<3;i++) normal[i]=normal[i]/(norm+1e-10);
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26997)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26998)
@@ -4679,6 +4679,6 @@
 	norm=sqrt(vector[0]*vector[0] + vector[1]*vector[1]);
 
-	normal[0]= + vector[1]/norm;
-	normal[1]= - vector[0]/norm;
+	normal[0]= + vector[1]/(norm+1e-10);
+	normal[1]= - vector[0]/(norm+1e-10);
 }
 /*}}}*/