1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
#!/usr/bin/env python
import sys
for i in range(0, len(sys.argv)):
if sys.argv[i] == '-A' and i < len(sys.argv)-1:
sys.path = sys.path + [sys.argv[i+1]]
import vtk
from vtk.util.misc import vtkGetDataRoot
filename = vtkGetDataRoot() + '/Data/tetraMesh.vtk'
reader = vtk.vtkUnstructuredGridReader()
reader.SetFileName(filename)
a = vtk.vtkMeshQuality()
a.SetInputConnection(reader.GetOutputPort())
a.VolumeOn()
a.RatioOn()
a.Update()
mesh = a.GetOutput().GetFieldData().GetScalars()
for i in range(mesh.GetNumberOfTuples()):
print mesh.GetTuple2(i)
sys.exit(0)
|
