an efficient C++ finite element environment
|
|
Go to the documentation of this file. 1 #include "neumann-laplace-assembly.h"
29 form m (Qh, Qh,
"mass");
33 form_diag dm (Qh,
"mass");
35 csr<Float>
c = neumann_laplace_assembly (
a.uu, mh.u);
38 vec<Float>
solve (
const vec<Float>& Mp)
const {
39 vec<Float> q1 =
fact_m.solve(Mp);
40 vec<Float> Mp_e (Mp.size()+1);
41 for (
size_t i = 0; i < Mp.size(); i++) Mp_e.at(i) = Mp.at(i);
42 Mp_e.at(Mp.size()) = 0;
43 vec<Float> q2_e =
fact_c.solve(Mp_e);
44 vec<Float> q2 (Mp.size());
45 for (
size_t i = 0; i < q2.size(); i++) q2.at(i) = q2_e.at(i);
46 vec<Float> q = q1 +
lambda*q2;
see the field 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
see the space 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 grad(const Expr &expr)
grad(uh): see the expression page for the full documentation
rheolef::details::is_vec dot
vec< Float > solve(const vec< Float > &Mp) const
see the test page for the full documentation
see the Float page for the full documentation
see the test page for the full documentation
cahouet_chabart(const space &Qh, Float lambda_1)