2 Aug
2018
2 Aug
'18
8:19 a.m.
On 1 Aug 2018, at 18:41, Amireh Mousavi <amireh.mousavi@math.iut.ac.ir> wrote:
sigma, u = w.split() Re = assemble((div(sigma)+f)*(div(sigma)+f)*dx).array()
Here you assemble a functional. So it is just a number. If you want the cellwise value of this expression write: P0 = FunctionSpace(mesh, "DG", 0) Re = project(expr, P0) Where expr is the integrals of your expression above. Thanks, Lawrence