ExodusReader only reads float data (not double)
Given a file with double
precision point coordinate data, e.g., rectangle_with_hole.e, the vtkExodusIIReader()
only reads the data in float
precision.
MWE:
import vtk
reader = vtk.vtkExodusIIReader()
reader.SetFileName('/tmp/rectangle_with_hole.e')
reader.Update()
out = reader.GetOutput()
mesh = out.GetBlock(0).GetBlock(0)
pt_data = mesh.GetPoints().GetData()
print(pt_data)
vtkFloatArray (0x1f29cf0)
[...]
Data type: float
[...]