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