The sinus product function – reconstruction for the hybrid high order method
int main(
int argc,
char**argv) {
Float err_us_s_h1_expected = (argc > 1) ? atof(argv[1]) : 1e+38;
geo omega = uhs.get_geo();
space Xhs = uhs.get_space();
size_t d = omega.dimension(), kp1 = Xhs.degree();
derr <<
"err_us_l2 = " << err_us_l2 << endl
<< "err_us_h1 = " << err_us_h1 << endl;
#ifdef TODO
derr <<
"err_us_s_l2 = " << err_us_s_l2 << endl
<< "err_us_s_h1 = " << err_us_s_h1 << endl;
#endif // TODO
return (err_us_h1 < err_us_s_h1_expected) ? 0 : 1;
}