The yield slip problem on a circle – error computation
int main(
int argc,
char**argv) {
Float tol = (argc > 1) ? atof(argv[1]) : 1e-15;
const geo& omega = uh.get_geo();
const space& Xh = uh.get_space();
dout <<
"err_linf = " << err_linf << endl
<< "err_l2 = " << err_l2 << endl
<< "err_h1 = " << err_h1 << endl
<< "err_b = " << err_b << endl
<< "err_lb = " << err_lb << endl;
return (err_linf < tol) ? 0 : 1;
}