← Back to team overview

dolfin team mailing list archive

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