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] ------------------------------------------------------------------------------ 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
