Hi,
I need to solve the heterogeneous diffusion equation with the DG SWIP method. To this aim, I need to construct a new weighted mean operator on interior interfaces where the weighting coefficients depend on the (constant) value of the diffusivity defined on both sides of the interface. How to access the coefficient value of the neighboring element? Can I access it via the .Other() command like this:
k = CoefficientFunction(....)
w1=k/(k+k.Other())
w2=k.Other()/(k+k.Other())
wmean_gradu = w1*grad(u)+w2*grad(u.Other())
Any idea?
Regards
Vincent
