dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

I'm interesting in accessing the vector of a Function after mesh-refinement.
At the moment, I'm getting:  

   RuntimeError: *** Error: You are attempting to access a non-const
   vector from a sub-Function.

Is this inevitable, a bug or a trigger-happy error message?

   from dolfin import *
   mesh = UnitSquare(2,2)
   V = FunctionSpace(mesh, "CG", 1)
   u = Function(V)
   mesh.refine()
   u.vector()


--
Marie