pyvista 0.45 temporary fix

With Pyvista 0.45

pv::from_meshio

had added

# Convert cell sets to cell data

if mesh.cell_sets:

mesh.cell_sets_to_data()

But apparently this function is only working, when the cell_set (e.g. "bottom", "top", "rect") are disjunct data sets. This does not hold for data we prepare for OGS

The problem is here:

arr = np.full(len(self.cells[k]), default_value, dtype=int)

for i, cc in enumerate(c):

if cc is None:

continue

arr[cc] = i

intfun.append(arr)

where the index cc, is too large for the size of arr.

self.cells[k] is working on the local/subset/set while .cell_sets.values() is associated to the global/all.

• Clarify if report to pyvista OR if this is really something very specific to OGS
• meanwhile workaround with this MR