Dear Sir,
Thank you for your reply to my previous questions.
I have one more doubt regarding function space.
is it possible to produce the same answer if I use V=vectorfunctionspace(mesh, “CG”, 1) instead
Q=Finiteelement(“CG”, mesh, 1)
V=functionspace(mesh, mixedelement([Q, Q] )) …
advance thanks for answer