Increasing the function-space order for Lagrange elements

[meteoguc]

Actually, I am using this mesh that you have coded for the most part:

2D structured mesh for hollow cylinder

Make sure you use the new_mesh, and not mesh the old mesh. You could start by creating a very coarse mesh and tabulating the dof coordinates, and visually inspect them.

But for FiniteElement(‘P’, triangle, p_h) for p_h > 3 it does not give me the results that I expect. I am assuming that my variational formulation is right. Even if it is wrong, that does not make sense to me that changing the order is not affecting the results at all. So I thought may be this is a common problem people are having with Lagrange elements when I found that topic.