@@ -244,19 +244,19 @@ def get_volume_subdivided(element_coordinates, center):
244244 #print(points)
245245
246246# debugging elementedgecouple based volume for maxtet
247- maxtet = 0
248- element_model (device = device , region = region , name = 'DerivedEdgeVolume' , equation = '1./6. * ElementEdgeCouple * EdgeLength' )
249- evol = get_element_model_values (device = device , region = region , name = 'DerivedEdgeVolume' )
250- print (evol [6 * maxtet :6 * maxtet + 6 ])
251- print (sum (evol [6 * maxtet :6 * maxtet + 6 ]))
252- print (element_node_volumes [6 * maxtet :6 * maxtet + 6 ])
253- print (sum (element_node_volumes [6 * maxtet :6 * maxtet + 6 ]))
254- print (actual_volumes [maxtet ])
255- print (tetrahedron_volumes [maxtet ])
256-
257- print (coordinate )
258- print (elements )
259- print ('start centers' )
260- for e in elements :
261- center , radius = get_center (get_element_coordinates (coordinate , e ))
262- print (center )
247+ # maxtet = 0
248+ # element_model(device=device, region=region, name='DerivedEdgeVolume', equation='1./6. * ElementEdgeCouple * EdgeLength')
249+ # evol = get_element_model_values(device=device, region=region, name='DerivedEdgeVolume')
250+ # print(evol[6*maxtet:6*maxtet+6])
251+ # print(sum(evol[6*maxtet:6*maxtet+6]))
252+ # print(element_node_volumes[6*maxtet:6*maxtet+6])
253+ # print(sum(element_node_volumes[6*maxtet:6*maxtet+6]))
254+ # print(actual_volumes[maxtet])
255+ # print(tetrahedron_volumes[maxtet])
256+ #
257+ # print(coordinate)
258+ # print(elements)
259+ # print('start centers')
260+ # for e in elements:
261+ # center, radius = get_center(get_element_coordinates(coordinate, e))
262+ # print(center)
0 commit comments