Many thanks. I will look a those.
By the way, are there some need in libmesh to make some
partitioning with overloaping between the various subdomains .
I mean, if the idea is to make some flux volume method or any finite
element method
that has stencil which implies needs of information from the
neighbourhing tetrahedra
then overlapping is necessary.
SO, do the example that you are pointing out considering consider
those problems ?
Thanks again for your precious help
Tahar
--------------------------------------------
T. Amari
Centre de Physique Theorique
Ecole Polytechnique
91128 Palaiseau Cedex France
tel : 33 1 69 33 47 53
fax: 33 1 69 33 30 08
email: <mailto:[EMAIL PROTECTED]>
URL : http://www.cpht.polytechnique.fr/cpht/amari
Le 18 nov. 07 à 15:27, Benjamin Kirk a écrit :
Hello,
I am looking for the piece of MPI/ParMetis code that distribute the
Tetrahedral Mesh
over the processors and also the associated unknown to "update".
Please note that I am looking for the low level ones coded in Libmesh
and not
the high level ones.
Could you please let me know where I could find those ?
./contrib/parmetis/Lib is probably what you want.
There are number of partitioning schemes in there -- xyzpart.c uses a
space-filling curve key and sorts it in parallel to get a
partitioning of an
abitrary unstructured mesh.
Parmetis does not distribute the mesh for you, though. That is
onoing work
in libMesh, but MeshCommunication::delete_remote_elements() in
src/mesh/mesh_communication.C will give you an idea of how this
might work.
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
