Make the FE code work for the H.pseudo.gaussian