Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
burgers_diffusion_operators.icc
Go to the documentation of this file.
1
25field lh (Float epsilon, Float t, const test& v) {
26#ifdef NEUMANN
27 return field (v.get_vf_space(), 0);
28#else // NEUMANN
29 size_t d = v.get_vf_space().get_geo().dimension();
30 size_t k = v.get_vf_space().degree();
31 Float beta = (k+1)*(k+d)/Float(d);
32 return epsilon*integrate ("boundary",
33 beta*penalty()*g(epsilon,t)*v
34 - g(epsilon,t)*dot(grad_h(v),normal()));
35#endif // NEUMANN
36}
37field gh (Float epsilon, Float t, const field& uh, const test& v) {
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",
42 compose (phi, normal(), uh, g(epsilon,t))*v);
43}
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
Definition cavity_dg.h:29
Definition cavity_dg.h:25
Definition phi.h:25
Float epsilon