Maybe the definition of `hTinv` should be ``` hTinv = CellField(1 ./ collect(get_array(∫(1)dΩp) .^ (1/2)), Ωp) ``` i.e. with a power 1/2, since it is a 2D problem.
Maybe the definition of
hTinvshould bei.e. with a power 1/2, since it is a 2D problem.