Selectively under-integrate coordinates: still oscillations in J

mrp089 · Jan 31, 2024 · e5ed2de · e5ed2de
1 parent 2070435
commit e5ed2de
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/Code/Source/svFSI/nn_elem_gip.h b/Code/Source/svFSI/nn_elem_gip.h
@@ -73,6 +73,11 @@ GetElementGausIntMapType get_element_gauss_int_data = {
     xi(0,7) = t; 
     xi(1,7) = s; 
     xi(2,7) = s;
+    for (int i = 0; i < 8; i++) {
+      // xi(0,i) = 0.0;
+      xi(1,i) = 0.0;
+      xi(2,i) = 0.0;
+    }
     }
   },
 
@@ -292,6 +297,11 @@ SetElementGausIntMapType set_element_gauss_int_data = {
       mesh.xi(0,5) = s; mesh.xi(1,5) = t; mesh.xi(2,5) = s;
       mesh.xi(0,6) = s; mesh.xi(1,6) = s; mesh.xi(2,6) = s;
       mesh.xi(0,7) = t; mesh.xi(1,7) = s; mesh.xi(2,7) = s;
+      for (int i = 0; i < 8; i++) {
+        // mesh.xi(0,i) = 0.0;
+        mesh.xi(1,i) = 0.0;
+        mesh.xi(2,i) = 0.0;
+      }
     }
   },