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 _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
