dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Mon, Nov 22, 2010 at 05:45:07PM +0100, Andre Massing wrote:
> Den 22.11.10 17.40, skrev Johan Hake:
> >On Monday November 22 2010 07:36:25 Andre Massing wrote:
> >>Hi!
> >>
> >>I am not entirely sure, whether there was a similar discussion about that,
> >>   but I just recognized that using a ufl file like
> >>
> >>element = FiniteElement("Lagrange","tetrahedron", 1)
> >>w = Coefficient(element)
> >>M = w*w*dx(0)
> >>
> >>and a meshfunction with 0, 1, 2 values
> >>results into the assembly of a function on *entire* domain and not only
> >>cells that a marked with 0.
> >>Can anybody confirm this behaviour?

Take a look in Assembler.cpp in the function assemble_cells. It should
be quite easy to find what goes wrong with some suitable print
statements. Look for the if-statements starting with

  if (domain && domains->size() > 0)

--
Anders