dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

Garth N. Wells wrote:
> Marie Rognes wrote:
> > I have a slightly interesting issue, possibly with interpolate:
> >
> > Say I have a RT_0 function u. I want to interpolate this function 
> > onto a DG_1 function Pi_u.
> > Since RT_0 \subset DG_1, I would expect to get the same function. 
> > This does not seem
> > to be the case at the moment.
>
> Need to be careful interpolating in discontinuous spaces because you 
> will pick up the 'last' evaluated term, i.e. you can't be sure which 
> on side the function will be evaluated on. This could be the issue. 
> Have a look inside the interpolate functions in Function.cpp.


Ah. This makes the use of interpolate in errornorm pretty unpredictable ...


--
Marie