27 integrate_option iopt = integrate_option())
29 size_t k = Xh.degree();
30 size_t d = Xh.get_geo().dimension();
34 a = integrate (2*ddot(Dh(
u),Dh(v)), iopt)
35 + integrate (
"sides", beta*penalty()*dot(jump(
u),jump(v))
36 - 2*dot(jump(
u),average(Dh(v)*normal()))
37 - 2*dot(jump(v),average(Dh(
u)*normal())), iopt);
38 lh = integrate (dot(
f(),v), iopt)
39 + integrate (
"boundary", beta*penalty()*dot(
g(),v)
40 - 2*dot(
g(),Dh(v)*normal()), iopt);
41 b = integrate (dot(
u,grad_h(q)), iopt)
42 + integrate (
"internal_sides", - dot(average(
u),normal())*jump(q), iopt);
43 kh = integrate (
"boundary", dot(
g(),normal())*q, iopt);
44 c = integrate (
"internal_sides", h_local()*jump(
p)*jump(q), iopt);
45 mp = integrate (
p*q, iopt);
field lh(Float epsilon, Float t, const test &v)
see the Float page for the full documentation
see the field page for the full documentation
see the space page for the full documentation
see the test page for the full documentation
see the test page for the full documentation
void stokes_dirichlet_dg(const space &Xh, const space &Qh, form &a, form &b, form &c, form &mp, field &lh, field &kh, integrate_option iopt=integrate_option())