dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, Apr 08, 2011 at 11:34:27AM +0200, Neilen Marais wrote:
> Hi,
>
> While experimenting with dolfin in python for my own education, I
> instrumented a SubDomain() subclass to count the number of calls to
> inside(). What I found is that for a tet mesh with 29 cells, the
> inside() function is called 44 times when doing:
>
> subdomain.mark(mesh_function, value)
>
> Just out of interest, what are the extra calls to inside() for?

All vertices + the midpoint is checked. That doesn't seem to add up
with your numbers, but take a look in

  dolfin/mesh/SubDomain.cpp

and you should be able to figure out how it works.

--
Anders