dolfin team mailing list archive

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

Johan Hake wrote:
> On Thursday 02 July 2009 12:58:25 Garth N. Wells wrote:
> > Marie Rognes wrote:
> > > The following code gives r = 0.0. It is not supposed to be.
> > >
> > > The problem seems to be that f's vector is still all zeros at the call
> > > to interpolate. Could this be easily fixed?
> > This example should have led to an error message since f is not a
> > discrete function. I'll take a look.
>
> Why should f be discrete? 

It is "discrete" in that it has a finite element type, but it is not 
"discrete" in the the sense that it has a vector of dofs.

Garth

f is interpolated into Q_h. Pi_f is the one that has
> to be discrete. Isn't f's eval used in the interpolation when it is not 
> discrete?
>
> ... or have I got lost in the interpolation jungle?
>
> Johan
>
> > Garth
> >
> > > --
> > > from dolfin import *
> > >
> > > mesh = UnitSquare(1,1)
> > > V_h = FunctionSpace(mesh, "CG", 1)
> > > f = Function(V_h, "1.0")
> > >
> > > Q_h = FunctionSpace(mesh, "DG", 1)
> > > Pi_f = interpolate(f, Q_h)
> > >
> > > r = norm(Pi_f)
> > > print "r = ", r
> > > ---
> > >
> > > --
> > > Marie
> > >
> > >
> > > _______________________________________________
> > > DOLFIN-dev mailing list
> > > DOLFIN-dev@xxxxxxxxxx
> > > http://www.fenics.org/mailman/listinfo/dolfin-dev
> > _______________________________________________
> > DOLFIN-dev mailing list
> > DOLFIN-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/dolfin-dev