- bugfix

Author: janbender

Utils/IndexedFaceMesh.h

@@ -187,6 +187,9 @@ namespace Utilities
 
 				m_normals[i] = v1.cross(v2);
 				m_normals[i].normalize();
+				// fix normals of degenerate triangles that can become zero vectors
+				if (m_normals[i].squaredNorm() < 1e-6f)
+					m_normals[i] = Vector3r::UnitX();
 			}
 		}
 	}

Utils/VolumeIntegration.cpp

@@ -140,7 +140,8 @@ void VolumeIntegration::face_integrals(unsigned int f)
 
   w = m_weights[f];
   n = m_face_normals[f];
-  k1 = 1 / n[C]; k2 = k1 * k1; k3 = k2 * k1; k4 = k3 * k1;
+  k1 = (n[C] == 0) ? 0 : 1 / n[C];
+  k2 = k1 * k1; k3 = k2 * k1; k4 = k3 * k1;
 
   Fa = k1 * Pa;
   Fb = k1 * Pb;