Changeset 26998

Timestamp:

05/08/22 04:50:49 (3 years ago)

Author:

vverjans

Message:

CHG adding +1e-10 to norm in denominator to avoid norm==0 in case of mixed NaN and fixed levelset spc

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Penta.cpp (modified) (1 diff)
• Tetra.cpp (modified) (1 diff)
• Tria.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -3197 +3197 @@
         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);
 }
 /*}}}*/

issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp

@@ -605 +605 @@
         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);
 }
 /*}}}*/

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -4679 +4679 @@
         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);
 }
 /*}}}*/