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