Defining weighted averages

[dokken]

If you consistently mark the domains \Omega_e with a cell function, the the restrictions s("+") and s("-") is consistently ordered, as shown in:

Integrating over an interior surface

Hi, consider using restriction using + and -. Here foo('+') restricts foo to the positive cell and FacetNormal('+') is the outer normal of a positive cell. Following the discussion here the +, - side can be changed by providing a form which has a dx measure term which has specified subdomain_data by some cell function g. Then a positive cell for a facet is the one for which value of g is higher. This allows for the following from dolfin import * from mshr import * omega1 = Rectangle(Point(0, 0…

and therefore you can define

def avg_gamma(u, gamma):
     return (0.5+gamma)*u("+") + (0.5-gamma)*u("-")