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

Reply via email to