Throw negative Jacobian error

mrp089 · Mar 3, 2024 · e05e3f9 · e05e3f9
1 parent 2ec754d
commit e05e3f9
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/Code/Source/svFSI/gr_equilibrated.cpp b/Code/Source/svFSI/gr_equilibrated.cpp
@@ -59,6 +59,8 @@ void stress_tangent_(const grModelType &grM, const double Fe[3][3],
 
   // determinant of the deformation gradient
   const double J = F.det();
+  if (J < 0.0)
+    throw std::runtime_error("[gr_equilibrated] Negative Jacobian");
 
   // set example
   enum Example { none, aneurysm, tortuosity, stenosis };