ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Jan 10, 2008 at 09:29:56AM -0500, Jake Ostien wrote:
> 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

I don't understand what you want to do. The only way I see to compute
a projection (or lifting) is to define a bilinear form and linear form
as you do and solve a linear system.

-- 
Anders