The cosinus product function – error analysis for the Poisson problem.
The cosinus product function – error analysis for the Poisson problem
int main(
int argc,
char**argv) {
Float err_u_linf_expected = (argc > 1) ? atof(argv[1]) : 1e+38;
geo omega = uh.get_geo();
space Xh = uh.get_space();
size_t k = Xh.degree();
size_t d = Xh.get_geo().dimension();
Float err_u_l2 = sqrt(integrate (omega, sqr(uh-
u_exact(
d)), iopt));
string opts = Xh.get_basis().option().stamp();
space Xh1 (omega,
"P"+to_string(k+1)+
"d"+opts);
Float err_u_linf = euh.max_abs();
Float err_u_h1 = sqrt(integrate (omega, norm2(grad_h(euh)), iopt));
dout << "err_u_l2 " << err_u_l2 << endl
<< "err_u_linf " << err_u_linf << endl
<< "err_u_h1 " << err_u_h1 << endl;
return (err_u_linf <= err_u_linf_expected) ? 0 : 1;
}
see the Float page for the full documentation
see the field page for the full documentation
see the geo page for the full documentation
see the environment page for the full documentation
see the integrate_option page for the full documentation
see the space page for the full documentation
This file is part of Rheolef.
rheolef - reference manual
The sinus product function.