ffc team mailing list archive

[Jehanzeb Hameed <jhameed2@xxxxxxxxxxxx>]

> Use one single mesh and define the integral over the interface using
> the standard dS and jump(), but when you assemble make sure to use the
> argument interior_facet_domains. You need to create a
> MeshFunction<uint> over interior facets to specify which facets are
> part of the interface.
>

Can you please tell me what the individual values for MeshFunction
should be like? E.g., say I can mark facets corresponding to
interface, will these be 1 in the mesh_function, and other facets set
to 0 ? I tried looking at the source code for assembler, which had
this line corresponding to interior_facet_domains, but I couldnt
understand what exactly it is doing:

if (domain < ufc.form.num_interior_facet_integrals())
       integral = ufc.interior_facet_integrals[domain];

Also, by standard dS, I take it you mean no integrals over subdomains
in the definition of the form (i.e. no dS1 etc).

I appreciate your patient answering of my questions,

Thanks,
-Jehanzeb