The tensorial transport benchmark – error computation
int main(
int argc,
char**argv) {
Float tol = (argc > 1) ? atof(argv[1]) : 1e-10;
size_t k = sigma_h.get_space().degree();
space Th1 (sigma_h.get_geo(),
"P"+
itos(3*(k+1))+
"d",
"tensor");
dout <<
"err_l2 = " << err_l2 << endl
<< "err_linf = " << err_linf << endl;
return (err_l2 < tol) ? 0 : 1;
}