19 Jul
2016
19 Jul
'16
11:05 a.m.
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))