dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, Feb 12, 2010 at 04:23:51PM +0100, Marie Rognes wrote:
>
> 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()

Looks like a bug. The function is not being refined but the mesh, the
dofmap and the functionspace is refined. Might be related to the
recent work on refinement.

--
Anders

Attachment: signature.asc
Description: Digital signature