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