Only quadratic laplacian, but it works
The 3D profile is different with the two gauss points, it has to be equal