an efficient C++ finite element environment
|
|
Go to the documentation of this file.
27 : We(0),
alpha(8./9),
a(1), Re(1), delta_t(0.025), tol(1e-6), max_iter(500),
28 Th(), Xh(), Qh(),
b(),
c(),
d(), mt(), inv_mt(),
mu(), mp(), th(), thb(),
29 theta(),
lambda(),
eta(),
nu(), c1(),
c2(), c3(), c4(), c5(), stokes() {}
32 Th =
space (omega,
"P1d",
"tensor");
33 Xh = P::velocity_space (omega,
"P2");
34 Qh =
space (omega,
"P1d");
36 lambda = Re/(theta*delta_t);
37 eta = ((1 -
alpha)*We + theta*delta_t)/(We + theta*delta_t);
38 nu = 1/((1-2*theta)*delta_t);
39 c1 = We/(We + theta*delta_t);
40 c2 = - We*theta*delta_t/(We + theta*delta_t);
41 c3 =
alpha*theta*delta_t/(We + theta*delta_t);
42 c4 = 1/((1-2*theta)*delta_t) - 1/We;
45 test v (Xh), xi (Th), q (Qh);
49 integrate_option iopt;
56 stokes.set_metric (mp);
std::enable_if< details::is_field_convertible< Expr >::value,details::field_expr_v2_nonlinear_terminal_field< typename Expr::scalar_type,typename Expr::memory_type,details::differentiate_option::divergence >>::type div(const Expr &expr)
div(uh): see the expression page for the full documentation
Float alpha[pmax+1][pmax+1]
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
see the problem_mixed page for the full documentation
see the space page for the full documentation
void reset(const geo &omega)
T ddot(const tensor_basic< T > &a, const tensor_basic< T > &b)
ddot(x,y): see the expression page for the full documentation
rheolef::details::is_vec dot
see the test page for the full documentation
std::enable_if< details::is_field_convertible< Expr >::value,details::field_expr_v2_nonlinear_terminal_field< typename Expr::scalar_type,typename Expr::memory_type,details::differentiate_option::gradient >>::type D(const Expr &expr)
D(uh): see the expression page for the full documentation.
see the test page for the full documentation
see the geo page for the full documentation