ufl team mailing list archive

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

Hello,

Should the following form be legal in UFL?

element = FiniteElement("Lagrange", triangle, 1)
element1 = VectorElement("Lagrange", triangle, 1)

v = TestFunction(element)
u = TrialFunction(element)

V = TestFunction(element1)
U = TrialFunction(element1)

f = Function(element)

a = inner(V, U)*dx + f*v*u*dx

It passes the validate_form() in checks.py

However, extracting basisfunctions yield:
[BasisFunction(VectorElement('Lagrange', Cell('triangle', 1), 1, 2), -2),
BasisFunction(VectorElement('Lagrange', Cell('triangle', 1), 1, 2), -1)]

redefining the form:

a = f*inner(V, U)*dx + v*u*dx

The extract basis yield:

[BasisFunction(FiniteElement('Lagrange', Cell('triangle', 1), 1), -2),
BasisFunction(FiniteElement('Lagrange', Cell('triangle', 1), 1), -1)]


Kristian