I’m solving a 2D heat equation, and I am saving the solution both as vtk and as numpy array, by using the .vector().get_local() method. Is it possible that the ordering of the two arrays, the one contained in the vtk file and the one saved via numpy is not the same?
Thanks.
Saving as vtk and saving as an array does fundamentally different things. This has for instance been covered in: Mapping from mesh.coordinates() to u.vector() - #2 by dokken
Without a minimal reproducable example it is hard to give any more concrete advice.
1 Like