The Helmholtz problem on a surface – error analysis
using namespace std;
int main (
int argc,
char**argv) {
Float tol = (argc > 1) ? atof(argv[1]) : 1e+38;
const space& Wh = uh.get_space();
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;
}