Subject: [UFL-dev] form integration

From: Kristian Oelgaard <[email-masked]>
Date: Fri, 05 Dec 2008 17:39:23 +0100

Hi,

I stumbled over the following which seems peculiar to me:

P2 = VectorElement("Lagrange", "triangle", 2)

DG = FiniteElement("Discontinuous Lagrange", "triangle", 0)

ME = MixedElement([P2, DG])

(v, q) = TestFunctions(ME)

(u, p) = TrialFunctions(ME)

a = q*p*dx

Now, FFC insists that the degree of the above form is 4 and it will therefore
use 3x3 integration points. For tensor representation this will lead to an
increase in FFC compile time, for quadrature an increase in runtime.

Is this a bug? Or something that happens for a good reason and which is hard to
do anything about? What will UFL do?

Kristian