Commit 02081973 authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

fix fluid vtk export for paraview 5.8

parent cd7d96e4
Pipeline #7054 passed with stage
in 2 minutes and 19 seconds
......@@ -267,7 +267,7 @@ class FluidProblem :
field_data = [(b"Boundary %s"%(name), nodes) for name,nodes in self._mesh_boundaries().items()]
connectivity = self.elements()
types = np.repeat([5 if self._dim == 2 else 10],connectivity.shape[0])
offsets = np.cumsum(np.repeat([self._dim+1],connectivity.shape[0]))
offsets = np.cumsum(np.repeat([self._dim+1],connectivity.shape[0])).astype(np.int32)
VTK.write(output_dir+"/fluid",it,t,(types,connectivity,offsets),self.coordinates(),data,field_data,cell_data)
def import_vtk(self, filename) :
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment