Yes, it then just hit another use of ref_el lower down. I can fix this of course, but I wonder if I am doing something weird since this is just the usual usage of BrokenElement and hasn't been an issue in other examples. --cjc On 5 July 2016 at 10:00, Mitchell, Lawrence < lawrence.mitchell@imperial.ac.uk> wrote:

...

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

-- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916