Hi all,
1) Given a 3D mesh, I try to write a separate code to add boundary id,
using:
Mesh::element_iterator it_el =
mesh.active_local_elements_begin();
const Mesh::element_iterator it_last_el =
mesh.active_local_elements_end();
for ( ; it_el != it_last_el ; ++it_el) {
Elem* elem = *it_el;
// some code to detect null-neighbour
mesh.boundary_info->add_side(elem->id(), side, bc_code);
equation_systems.init();
2) then I do the check:
mesh.boundary_info->print_info();
mesh.boundary_info->print_summary();
mesh.write("withBC.xda");
3) The question is that:
For serial run, the output results are correct, but for parallel run
(np=2), the output results show fewer boundary ids, e.g. only part of
boundary info is added on to the mesh??
Is mesh.boundary_info->print_info() and mesh.write() only output data on
processor 0?
Why?
4) What is typical way to add boundary id on the mesh.
Thanks,
Walter
------------------------------------------------------------------------------
Learn Graph Databases - Download FREE O'Reilly Book
"Graph Databases" is the definitive new guide to graph databases and their
applications. Written by three acclaimed leaders in the field,
this first edition is now available. Download your free book today!
http://p.sf.net/sfu/NeoTech
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
