Ah, I should have spotted that! thanks Lawrence --cjc On 19 July 2016 at 12:21, Bendall, Thomas <t.bendall15@imperial.ac.uk> wrote:

Fantastic, that seems to have worked! Thank you so much.

On 19/07/16 12:05, Lawrence Mitchell wrote:

On 19/07/16 11:47, T. M. Bendall wrote:

Qmean.assign((m * Qmean) + Q0) Qmean.assign(Qmean / (m + 1))

This line (and the psimean guy below it) is causing problems (at least). m is just an integer and this assign statement cause a new kernel to be generated with the value of the integer pasted in. Use m = Constant(m) m.assign(...)

There may be other issues.

Lawrence

Q1.assign(Q0) psi_solver.solve()

psimean.assign((m * psimean) + psi0) psimean.assign(psimean / (m + 1))

_______________________________________________ firedrake mailing listfiredrake@imperial.ac.ukhttps://mailman.ic.ac.uk/mailman/listinfo/firedrake

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