27 return field (v.get_vf_space(), 0);
29 size_t d = v.get_vf_space().get_geo().dimension();
30 size_t k = v.get_vf_space().degree();
32 return epsilon*integrate (
"boundary",
34 -
g(
epsilon,t)*dot(grad_h(v),normal()));
38 return - integrate (dot(compose(
f,uh),grad_h(v)))
39 + integrate (
"internal_sides",
40 compose (
phi, normal(), inner(uh), outer(uh))*jump(v))
41 + integrate (
"boundary",
field lh(Float epsilon, Float t, const test &v)
field gh(Float epsilon, Float t, const field &uh, const test &v)
see the Float page for the full documentation
see the field page for the full documentation
see the test page for the full documentation