Couple more questions, this time regarding weakly imposed BCs
In the undocumented FEniCS
dg-poisson example, they have this weak formulation:
a = dot(grad(v), grad(u))*dx \
- dot(avg(grad(v)), jump(u, n))*dS \
- dot(jump(v, n), avg(grad(u)))*dS \
+ alpha/h_avg*dot(jump(v, n), jump(u, n))*dS \
- dot(grad(v), u*n)*ds(1) \
- dot(v*n, grad(u))*ds(1) \
+ (gamma/h)*v*u*ds(1)
L = v*f*dx - u0*dot(grad(v), n)*ds(1) + (gamma/h)*u0*v*ds(1) + g*v*ds(2)
I am attempting to emulate this in Firedrake but with a diffusion tensor (let's say D). From what I understand in the literature, I know I have to add D to the first three terms in the bilinear form.