Hello everybody I'm trying to call the function DofMap::dof_indices on an element that is not in my processor and I get a exception thrown when it tries to call RemoteElem::n_nodes
My question is more about the details we need to take into account when trying to implement a function that can be run in distributed processors. For example, this function is going to check the neighbors for each element, and it will call the neighbors of the neighbors until a certain criteria is met. How would we implement this in parallel? Just calling the local iterators> MeshBase::const_element_iterator el = this >get_mesh().active_local_elements_begin(); const MeshBase::const_element_iterator end_el = this->get_mesh().active_local_elements_end(); Calling the neighbors of the neighbors could be a problem because they are too far from the processors border? Another question related with the parallelism in libmesh. If within a libmesh program, I call an external library, will I be calling it once per processor? For example, say that I have an optimization routine that I want to run it in only one processor, while I want the analysis in libmesh in parallel. How could I set this up? 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] ------------------------------------------------------------------------------ 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
