On 27 May 2013 22:06, Jan Blechta <[email protected]> wrote:
> How do I obtain FacetFunction
>   fd = mesh.domains().facet_domains()
> with redesigned MeshDomains?
>

  D = mesh.topology().dim()
  fd = mesh.domains().markers(D - 1)

'fd' is now not a MeshFunction (std::map from C++ and a dict on the
Python side). It was necessary to change the type to remove a circular
dependency that was making new developments exponentially more
complicated and to make memory management robust. A lot of boilerplate
code was removed.

If it turns out that a MeshFunction is required, we could add a
MeshFunction constructor to make it simple to make a MeshFunction.

Garth

> Jan
> _______________________________________________
> fenics mailing list
> [email protected]
> http://fenicsproject.org/mailman/listinfo/fenics
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics

Reply via email to