I guess the problem is that there is no rhs in E_du. To circumvent this, you could add:
inner(Constant(0,0), v)*dx
to E_du.
1 Like
I guess the problem is that there is no rhs in E_du. To circumvent this, you could add:
inner(Constant(0,0), v)*dx
to E_du.