Hi,

I am trying to implement an HDG+ method. I am having trouble with a term involving the L^2projection of an unknown.

For example, for an element K [tex]u_h \in \mathcal P_{p+1}(K),\quad\mu\in \Pi_{F\in \partial K}\mathcal P_p(F)[/tex]. I would like to compute

[tex]\int_{\partial K} P_{M} uh \mu ds[/tex]

where P_M is the L^2 projection onto the space [tex]\mathcal P_{p+1}(K)[/tex].

Is possible to implement L^2 projection as operators?

I would appreciate any help.

Thanks.