> On 5 Jul 2016, at 09:41, Colin Cotter <colin.cotter@imperial.ac.uk> wrote:
>
> Actually, I'm now less sure where the problem is.
>
> If I do this, in the state object
> ...
> self.V_elt[2] = V2t_elt
> ...
>
> and then later,
> Vtdg_elt = BrokenElement(state.V_elt[2])
> Vtdg = FunctionSpace(mesh, Vtdg_elt)
>
> then I get the same error.
Did you try my suggestion?
Lawrence