Chris has been working on adding a layer of ghost cells/mesh entities to a
DOLFIN Mesh when running in parallel. The branch is at:
https://bitbucket.org/fenics-project/dolfin/branch/chris/ghost-mesh
The motivations for this are:
1. Support DG methods in parallel.
2. Communication-less assembly in parallel with MPI (GenericTensor::apply()
will involve no communication).
3. Simpler construction of dofmaps (less MPI communication).
4. Simple detection of entities on a domain boundary.
5. Eliminates the need for a boundary mesh to build a distributed mesh.
In the first instance, the plan is that the mesh iterators will iterate over
the same entities they do now, ignoring ghost entities. Existing parallel code
should just run as it does now. Longer term, it would be nice to support
iterators over different mesh views. There are quite a few problems for which
iteration over mesh views would be handy (in serial and in parallel).
Please send a message if you have any comments or suggestions on the above. I
anticpate the code being ready to merge not long after the 1.3 release of
DOLFIN.
Garth
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
