Hello I'm trying to output my mesh into a file with Mesh::write ("mesh.xda"). Because my mesh is parallel I call Mesh::allgather () first. Then I call Mesh::write("mesh.xda") which doesn't give me problems if I don't make any uniform refinement, however if I do, I obtain this assertion
Assertion `(3*local_idx+2) < coords.size()' failed. (3*local_idx+2) = 38 coords.size() = 33 At the moment of the assertion, my mesh consists of eight elements in a 4x2 grid, coming from a 2x1 mesh before one single uniform refinement. Any idea? I tried to follow adjoint_ex3 and added mesh.renumber_nodes_and_elements () before calling the write, but it didn't make a difference. Thanks Miguel ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users