For some strange reason the renumbering of node connectivity for VTK quads and hexahedrals is separated between VKT_interfaces.F90, vtkfortran.cpp and vtkmeshio.cpp. This leads to generating wrong vtus when creating a quad/hexahedral mesh from within Fluidity and then outputting it.
For example this results into the following bug if one decides to create a quad mesh and output it into a VTU
while the output mesh should be looking like this
I already have a fix for this and I will open a PR soon, but in general I think that libvtkfortran needs a revision and a refactoring to modernise the code (maybe a job for another time). I will make sure I add a unittest testing as many VTK elements I can.
For some strange reason the renumbering of node connectivity for VTK quads and hexahedrals is separated between VKT_interfaces.F90, vtkfortran.cpp and vtkmeshio.cpp. This leads to generating wrong vtus when creating a quad/hexahedral mesh from within Fluidity and then outputting it.
For example this results into the following bug if one decides to create a quad mesh and output it into a VTU
while the output mesh should be looking like this
I already have a fix for this and I will open a PR soon, but in general I think that libvtkfortran needs a revision and a refactoring to modernise the code (maybe a job for another time). I will make sure I add a unittest testing as many VTK elements I can.