dolfin team mailing list archive

[Gerd Wachsmuth <gerd.wachsmuth@xxxxxxxxxxxxxxxxxxxxxxxxx>]

** Attachment added: "python script with the bug"
   https://bugs.launchpad.net/bugs/660334/+attachment/1692078/+files/demo.py

-- 
Assembling over subdomain fails
https://bugs.launchpad.net/bugs/660334
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.

Status in DOLFIN: New

Bug description:
Following setup:

I have a mesh (with FunctionSpace V) and a submesh (with FunctionSpace W). Assembling a linear form on V over this subdomain involving a coefficient from W will fail.

I think the problem is that Dolfin would like to evaluate the function from W at the boundary of the subdomain. Maybe due to rounding errors this point lies outside the subdomain and the evaluation fails.

It does not happen, if W is a CG or DG space.