ufl team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

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