dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #07830
Re: assembly on sub domains
On Wed, May 07, 2008 at 05:33:50PM +0100, Nuno David Lopes wrote:
> Ok i've got it, i think this shoud be the way:
>
> void dolfin::assemble(GenericTensor& A, Form& form, Mesh& mesh,
> const MeshFunction<dolfin::uint>& cell_domains,
> const MeshFunction<dolfin::uint>&exterior_facet_domains,
> const MeshFunction<dolfin::uint>& interior_facet_domains) ?
>
> But still, I have to pass 3 MeshFunction's ? even if we only need one,
> (but that isn't a problem).
Yes, that's unfortunate but I don't what the interface should look
for specifying just one of them.
If you use the Python interface, you can use keyword arguments to do
this (as I learnt myself from Martin earlier today...):
a = ...*dx
A = assemble(a, mesh, exterior_facet_domains=...)
--
Anders
Follow ups
References