dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

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.

Ok, thanks!

However,

(a) Why is f not a discrete function? (It is defined on a finite element 
space?)
(b) I would very much like this not to give an error, but rather work ;) 
I would use something like this all the time when checking convergence.

--
Marie


> 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