I ended up adding the boundary conditions in the file instead of inside the program. More handy I guess.
On Fri, Aug 29, 2014 at 9:19 PM, Miguel Angel Salazar de Troya < [email protected]> wrote: > I've realized that the error comes from here > > dof_id_type last_elem = mesh.n_active_elem()-1; > mesh.boundary_info->add_side(last_elem,2,0); > > Basically I want to add a boundary condition to the top of the last > element. Can I only do it calling it from its processor? > > > On Fri, Aug 29, 2014 at 8:45 PM, Miguel Angel Salazar de Troya < > [email protected]> wrote: > >> I have this simple mesh in .xda format: >> >> DEAL 003:003 >> 4 # Num. Elements >> 10 # Num. Nodes >> 16 # Length of connectivity vector (Num elements * Num Nodes) >> 0 # Num. Boundary Conds. >> 65536 # String Size (ignore) >> 1 # Num. Element Types. We just have quads, no triangles >> 5 # Element types in each block. >> 4 # Num. of elements in each block. >> Id String >> Title String >> 0 1 4 5 >> 1 2 3 4 >> 5 4 6 7 >> 7 6 8 9 >> 0. 0. 0. >> 1. 0. 0. >> 2. 0. 0. >> 2. 1. 0. >> 1. 1. 0. >> 0. 1. 0. >> 1. 2. 0. >> 0. 2. 0. >> 1. 3. 0. >> 0. 3. 0. >> >> I can read it in serial, but in parallel (tried from 2 to 4) I get this >> error >> >> Assertion `_elements[i]' failed >> >> I think it must be something simple and wrong with my .xda file, but I >> have no idea. Any clue? >> >> Thanks >> Miguel >> >> -- >> *Miguel Angel Salazar de Troya* >> Graduate Research Assistant >> Department of Mechanical Science and Engineering >> University of Illinois at Urbana-Champaign >> (217) 550-2360 >> [email protected] >> >> > > > -- > *Miguel Angel Salazar de Troya* > Graduate Research Assistant > Department of Mechanical Science and Engineering > University of Illinois at Urbana-Champaign > (217) 550-2360 > [email protected] > > -- *Miguel Angel Salazar de Troya* Graduate Research Assistant Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign (217) 550-2360 [email protected] ------------------------------------------------------------------------------ Slashdot TV. Video for Nerds. Stuff that matters. http://tv.slashdot.org/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
