Re: [Libmesh-users] mesh.write()

2015-11-18 Thread Miguel Angel Salazar de Troya
That was the problem. I changed my mesh file cantilever.xda to follow the syntax in H_channel_quads.xda in adjoint_ex3 and that made it work. I remember reading a pdf document about xda files in libmesh, can't find it now. Could anybody explain me why my original mesh file wasn't working? Thanks M

Re: [Libmesh-users] mesh.write()

2015-11-18 Thread Miguel Angel Salazar de Troya
I think my problem must be in my original mesh format. This simple code has the same problem: int main (int argc, char** argv) { // Initialize libraries. LibMeshInit init (argc, argv); // Skip adaptive examples on a non-adaptive libMesh build #ifndef LIBMESH_ENABLE_AMR

Re: [Libmesh-users] mesh.write()

2015-11-18 Thread Miguel Angel Salazar de Troya
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 c

Re: [Libmesh-users] mesh.write()

2015-11-18 Thread Roy Stogner
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 ref

[Libmesh-users] mesh.write()

2015-11-18 Thread Miguel Angel Salazar de Troya
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