← Back to team overview

ufl team mailing list archive

Restrictions at facets

 

FFC used to throw an error if functions were integrated on interior facets but not restricted to the '+' or '-' side, but UFL doesn't. If I do

    element = FiniteElement("Discontinuous Lagrange", triangle, 1)
    v = TestFunction(element)
    u = TrialFunction(element)

    # Integrate on interior facet
    a = v*u*dS

what result should I expect? I had anticipated an error.

Garth


Follow ups