It's been a long time coming, but after the good work of Chris
Richardson on distributed mesh ghosting and a re-write of the dofmap
construction code, support for integration on interior facets (i.e. DG
methods) in parallel is now in the DOLFIN 'next' branch. All DG unit
tests and demos are now run in parallel as part of the test system.
To evaluate interior facet integrals in parallel, it's necessary to
construct a mesh with ghost cells. For now, this is somewhat clumsily
managed through the parameter system,
# Make mesh ghosted to permit evaluation of DG terms
# ("shared_facet", "shared_vertex" or "none")
parameters["ghost_mode"] = "shared_facet"
# Create mesh and define function space
mesh = UnitSquareMesh(32, 32)
Providing a better interface for mesh construction is the topic of an
ongoing discussion [1,2].
We'll move the new functionality into master soon, but it would useful
if it could get some testing first in next.
Garth
[1] https://bitbucket.org/fenics-project/dolfin/issue/67
[2]
https://bitbucket.org/fenics-project/dolfin/issue/67/move-building-of-foomesh-to-factory
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
