more debug

tcaduser · tcaduser · commit 8b5fc1d708f3 · 2023-11-17T17:17:56.000-06:00
diff --git a/tetgenmesh/check_tetrahedra.py b/tetgenmesh/check_tetrahedra.py
@@ -253,3 +253,10 @@ def get_volume_subdivided(element_coordinates, center):
     print(sum(element_node_volumes[6*maxtet:6*maxtet+6]))
     print(actual_volumes[maxtet])
     print(tetrahedron_volumes[maxtet])
+
+    print(coordinate)
+    print(elements)
+    print('start centers')
+    for e in elements:
+        center, radius = get_center(get_element_coordinates(coordinate, e))
+        print(center)
