Skip to content

Commit

Permalink
remove some debugging stuff
Browse files Browse the repository at this point in the history
  • Loading branch information
tcaduser committed Nov 18, 2023
1 parent 8b5fc1d commit 24f1280
Showing 1 changed file with 16 additions and 16 deletions.
32 changes: 16 additions & 16 deletions tetgenmesh/check_tetrahedra.py
Original file line number Diff line number Diff line change
Expand Up @@ -244,19 +244,19 @@ def get_volume_subdivided(element_coordinates, center):
#print(points)

# debugging elementedgecouple based volume for maxtet
maxtet = 0
element_model(device=device, region=region, name='DerivedEdgeVolume', equation='1./6. * ElementEdgeCouple * EdgeLength')
evol = get_element_model_values(device=device, region=region, name='DerivedEdgeVolume')
print(evol[6*maxtet:6*maxtet+6])
print(sum(evol[6*maxtet:6*maxtet+6]))
print(element_node_volumes[6*maxtet:6*maxtet+6])
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)
# maxtet = 0
# element_model(device=device, region=region, name='DerivedEdgeVolume', equation='1./6. * ElementEdgeCouple * EdgeLength')
# evol = get_element_model_values(device=device, region=region, name='DerivedEdgeVolume')
# print(evol[6*maxtet:6*maxtet+6])
# print(sum(evol[6*maxtet:6*maxtet+6]))
# print(element_node_volumes[6*maxtet:6*maxtet+6])
# 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)

0 comments on commit 24f1280

Please sign in to comment.