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