I realized that my nodes id's are not continuous. Having 15 nodes in my
mesh, their id's are 0,1,2,3,4,5,6,8,10,12,14,16,18,20,22, hence why I get
that assertion error, which assumes that the nodes' id are within the
limits of the number of nodes I think. I'm going to try to get an example
that I can send and that reproduces the error.
Miguel
On Nov 18, 2015 12:15, "Roy Stogner" <royst...@ices.utexas.edu> wrote:
>
> 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
