The Helmholtz problem on a surface – error analysis.
The Helmholtz problem on a surface – error analysis
int main (
int argc,
char**argv) {
Float tol = (argc > 1) ? atof(argv[1]) : 1e+38;
const space& Wh = uh.get_space();
form m = integrate (
u*v);
form a = integrate (dot(grad_s(
u),grad_s(v)));
size_t d = Wh.get_geo().dimension();
dout << "err_l2 " << sqrt(m(eh,eh)) << endl
<< "err_h1 " << sqrt(a(eh,eh)) << endl
<< "err_linf " << eh.max_abs() << endl;
return (eh.max_abs() < tol) ? 0 : 1;
}
see the Float page for the full documentation
see the field page for the full documentation
see the environment page for the full documentation
see the space page for the full documentation
see the test page for the full documentation
see the test page for the full documentation
This file is part of Rheolef.
rheolef - reference manual
The level set function for the sphere geometry.