ffc team mailing list archive

[Jake Ostien <tostien@xxxxxxxxx>]

Hi,

I'd like to write a form based on the lifting operator.  If the lifting 
operator (R) is defined as

   -\int_Gamma avg(A):jump(B) dS = \int_Omega A:R(B) dV

Is there any way I can define a form such as

   \int_Omega R(A):R(B) dV

Where e.g. A and B are BasisFunctions?

I can do this already for a known Function, say H, where I say

   V = TestFunction
   U = TrialFunction

   H = Function

   a = dot(V,U)*dx
   L = dot(avg(V),jump(H))*dS

Then U is the lifted H, projected onto the basis of V.

Should this same approach work for BasisFunctions?

Jake