Merge pull request #937 from nschloe/vtk-fix2

VTK fix

Author: nschloe

meshio/vtk/_vtk.py

@@ -470,16 +470,17 @@ def _read_scalar_field(f, num_data, split, is_ascii):
         raise ReadError()
 
     if is_ascii:
-        data = numpy.fromfile(f, count=num_data, sep=" ", dtype=dtype)
+        data = numpy.fromfile(f, count=num_data * num_comp, sep=" ", dtype=dtype)
     else:
         # Binary data is big endian, see
         # <https://www.vtk.org/Wiki/VTK/Writing_VTK_files_using_python#.22legacy.22>.
         dtype = dtype.newbyteorder(">")
-        data = numpy.fromfile(f, count=num_data, dtype=dtype)
+        data = numpy.fromfile(f, count=num_data * num_comp, dtype=dtype)
         line = f.readline().decode("utf-8")
         if line != "\n":
             raise ReadError()
 
+    data = data.reshape(-1, num_comp)
     return {data_name: data}
 
 
@@ -488,6 +489,7 @@ def _read_field(f, num_data, split, shape, is_ascii):
     data_type = split[2].lower()
 
     dtype = numpy.dtype(vtk_to_numpy_dtype_name[data_type])
+    # prod()
     # <https://stackoverflow.com/q/2104782/353337>
     k = reduce((lambda x, y: x * y), shape)
 

setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = meshio
-version = 4.2.1
+version = 4.2.2
 author = Nico Schlömer et al.
 author_email = [email protected]
 description = I/O for many mesh formats
@@ -24,6 +24,13 @@ classifiers =
     Programming Language :: Python :: 3.8
     Topic :: Scientific/Engineering
     Topic :: Utilities
+keywords =
+    mesh
+    file formats
+    scientific
+    engineering
+    fem
+    finite elements
 
 [options]
 packages = find:

test/meshes/vtk/gh-935.vtk

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:49d09f61cd621a26dcbbd715699462118b9b12cec4c73e64e50365ea0e351a07
+size 358

test/test_vtk.py

@@ -76,6 +76,7 @@ def test_reference_file(filename, ref_sum, ref_num_cells, binary):
         ("04_rectilinear.vtk", "quad", 27, 40),
         ("05_rectilinear.vtk", "quad", 27, 40),
         ("06_unstructured.vtk", "hexahedron", 12, 42),
+        ("gh-935.vtk", "triangle", 2, 6),
     ],
 )
 def test_structured(filename, ref_cells, ref_num_cells, ref_num_pnt):