The Mossolov problem for a circular pipe – error analysis
int main(
int argc,
char**argv) {
Float tol_u = (argc > 1) ? atof(argv[1]) : 1e-15;
Float tol_s = (argc > 2) ? atof(argv[2]) : 1e-15;
const geo& omega = uh.get_geo();
const space& Xh = uh.get_space();
const space& Th = sigma_h.get_space();
space Xh1 (omega,
"P" +
itos(2*Xh.degree()));
space Th1 (omega,
"P" +
itos(2*Xh.degree()) +
"d",
"vector");
Float err_u_linf = euh.max_abs();
Float err_s_linf = esh.max_abs();
dout <<
"err_u_linf = " << err_u_linf << endl
<< "err_u_l2 = " << err_u_l2 << endl
<< "err_u_h1 = " << err_u_h1 << endl
<< "err_s_linf = " << err_s_linf << endl
<< "err_s_l2 = " << err_s_l2 << endl;
return (err_u_linf < tol_u) && (err_s_l2 < tol_s) ? 0 : 1;
}