dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Anders Logg wrote:
> On Wed, Oct 07, 2009 at 09:08:53AM +0100, Garth N. Wells wrote:
> > We seem to have a problem with Python and interpolating functions. When
> > I try
> >
> >    P1 = FunctionSpace(mesh_post, "CG", 1)
> >    u_post = interpolate(u, P1);
> >
> > I get the below error.
> >
> > Garth
> >
> >
> > Traceback (most recent call last):
> >    File "demo.py", line 128, in <module>
> >      u_post = interpolate(u, P1);
> >    File
> > "/home/garth/code/fenics/dolfin/dolfin-all/local/lib/python2.6/site-packages/dolfin/interpolate.py",
> > line 28, in interpolate
> >      Pv.interpolate(v)
> > RuntimeError: *** Error: Treating Function as UFC function. Is that correct?
>
> Yes, it is correct.

OK, so it's a DOLFIN error.

> I assume you are interpolating u from another space, possibly DG?

Yes, it's another space. I want to interpolate a high-order Lagrange 
function onto a finer mesh with P1 elements.

> Then we should get to the else{} clause in the restrict() which
> triggers the error. I added the error there temporarily so that we can
> check that we only get to else{} when we should. Before the hash fix
> yesterday, we got there all the time which works but is very
> inefficient.
>
> I have not had time yet to run more than a couple of demos. Feel free
> to remove the error if you have checked that all demos run except the
> few (doing interpolation) where the error should be triggered.

I'll take a look if I have time.

Garth


> --
> Anders
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev