Rheolef  7.2
an efficient C++ finite element environment
burgers_diffusion_operators.icc

The diffusive Burgers equation – operators.

The diffusive Burgers equation – operators

field lh (Float epsilon, Float t, const test& v) {
#ifdef NEUMANN
return field (v.get_vf_space(), 0);
#else // NEUMANN
size_t d = v.get_vf_space().get_geo().dimension();
size_t k = v.get_vf_space().degree();
Float beta = (k+1)*(k+d)/Float(d);
return epsilon*integrate ("boundary",
- g(epsilon,t)*dot(grad_h(v),normal()));
#endif // NEUMANN
}
field gh (Float epsilon, Float t, const field& uh, const test& v) {
return - integrate (dot(compose(f,uh),grad_h(v)))
+ integrate ("internal_sides",
compose (phi, normal(), inner(uh), outer(uh))*jump(v))
+ integrate ("boundary",
compose (phi, normal(), uh, g(epsilon,t))*v);
}
u_exact g
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
class rheolef::details::field_expr_v2_nonlinear_node_unary compose
rheolef::details::is_vec dot
std::enable_if< details::has_field_rdof_interface< Expr >::value, details::field_expr_v2_nonlinear_terminal_field< typenameExpr::scalar_type, typenameExpr::memory_type, details::differentiate_option::gradient > >::type grad_h(const Expr &expr)
grad_h(uh): see the expression page for the full documentation
details::field_expr_v2_nonlinear_terminal_function< details::normal_pseudo_function< Float > > normal()
normal: see the expression page for the full documentation
std::enable_if< details::is_field_expr_v2_nonlinear_arg< Expr >::value &&!is_undeterminated< Result >::value, Result >::type integrate(const geo_basic< T, M > &omega, const Expr &expr, const integrate_option &iopt, Result dummy=Result())
see the integrate page for the full documentation
Definition: integrate.h:211
details::field_expr_v2_nonlinear_terminal_function< details::penalty_pseudo_function< Float > > penalty()
penalty(): see the expression page for the full documentation
Float beta[][pmax+1]
Definition: cavity_dg.h:29
Definition: phi.h:25
Float epsilon