dolfin team mailing list archive

[Marie Rognes <660334@xxxxxxxxxxxxxxxxxx>]

On 14. okt. 2010 09:13, Gerd Wachsmuth wrote:
> Public bug reported:
>
> 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.
>
> ** Affects: dolfin
>      Importance: Undecided
>          Status: New
>
>   


Have you tried setting

    parameters["allow_extrapolation"] = True

?

--
Marie

-- 
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.