How can I load a polygonal 2D mesh (VTK format) into JuAFEM?

Hi, I want to solve a reaction-diffusion equation on a 2D polygonal mesh, where each cell has a different parameter value for the reaction-diffusion equation. I have the mesh in the VTK format. Now I wonder whether I should simply write a script to convert the VTK data into a JuAFEM grid or whether there exists some IO functionality for that in the Julia package universe. Thanks for any hints.

One option is to make use of meshio if you don’t want to implement it yourself.
It can be easily called from julia via PyCall.

2 Likes