Hello,
I am trying to get to run my libmesh code in parallel. Thanks to the
help of you guys most of the code already works. However, I still have
some small problems where I still don't have enough understanding of
libmesh to solve them alone.
One of them is to do mesh refinement in parallel:
Mesh new_mesh;
TetGenIO TETGEN(new_mesh);
{ // read TETGEN file
OStringStream inmesh;
inmesh<<"geometry_wings/export_pov/out/"<<ReadMesh[actual_pulse]<<"_box.1.ele";
TETGEN.read(inmesh.str().c_str());
}
set_subdomains(&TETGEN, new_mesh); // works
set_boundary_ids(new_mesh); // works
MeshRefinement mesh_refinement (new_mesh); // see note below!!!
mesh_refinement.uniformly_refine(no_unirefine); // see note below!!!
new_mesh.prepare_for_use();
EquationSystems new_equation_systems (new_mesh);
The mesh.uniformly_refine(1) takes some minutes when I do it in serial
on my PC. However, on the BlueGene/P it doesn't come to an end (I
aborted after 15 min). Do I have to prepare my mesh in a special way for
use in parallel? One thing I have tried was to call the
prepare_for_use() before the refinement. I thought this makes sense
because the mesh is partitioned before. But, the problem remains.
Mesh Information:
Mesh Information:
mesh_dimension()=3
spatial_dimension()=3
n_nodes()=447273
n_local_nodes()=447273
n_elem()=2994336
n_local_elem()=2994336
n_active_elem()=2661632
n_subdomains()=2
n_processors()=1
processor_id()=0
Robert
------------------------------------------------------------------------------
Special Offer -- Download ArcSight Logger for FREE!
Finally, a world-class log management solution at an even better
price-free! And you'll get a free "Love Thy Logs" t-shirt when you
download Logger. Secure your free ArcSight Logger TODAY!
http://p.sf.net/sfu/arcsisghtdev2dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
