On Mon, 27 May 2013 23:10:47 +0100
"Garth N. Wells" <[email protected]> wrote:
> 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.

Indeed, MeshFunction is needed to feed DirichletBC or Measure.

Jan

> 
> Garth
> 
> > Jan
> > _______________________________________________
> > fenics mailing list
> > [email protected]
> > http://fenicsproject.org/mailman/listinfo/fenics
> _______________________________________________
> 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