dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

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?

--
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