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 list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake