On Wed, 18 Nov 2015, Miguel Angel Salazar de Troya wrote:
> 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?
Not a clue. Can you boil this down into an example you can share?
> I tried to follow adjoint_ex3 and added
> mesh.renumber_nodes_and_elements () before calling the write, but it didn't
> make a difference.
I'm not sure why adjoints_ex3 is doing that, but it shouldn't be
necessary.
---
Roy
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
