thanks for the suggestion Dokken, I will do that. To compute the integral analytically, I would need to know the form of a basis function. I asked a question a long time ago trying to go down this vein : Documentation on definition of basis functions over elements for CG and DG for arbitary degree of polynomial? . Do the links you listed have the functional form of an n-th degree polynomial over an element, or is there FENICS documentation?