26 integrate_option iopt = integrate_option())
29 integrate (- dot(trans(grad_h(v))*w,
u) - 0.5*div_h(v)*dot(
u,w), iopt)
30 + integrate (
"internal_sides",
31 dot(average(
u),normal())*dot(jump(v),average(w))
32 + 0.5*dot(jump(v),normal())
33 *(dot(average(
u),average(w)) + 0.25*dot(jump(
u),jump(w))), iopt)
34 + integrate (
"boundary", 0.5*dot(v,normal())*dot(
u,w), iopt);