Discontinuous Raviart-Thomas L2 projection – error analysis
int main(
int argc,
char**argv) {
Float err_p_l2_valid = (argc > 1) ? atof(argv[1]) : 1;
Float err_p_div_l2_valid = (argc > 2) ? atof(argv[2]) : 1;
const geo& omega = p_Vh_u.get_geo();
const space& Vh = p_Vh_u.get_space();
size_t d = omega.dimension();
size_t k = Vh.get_basis().family_index();
dout <<
"err_p_l2 = " << err_p_l2 << endl
<< "err_i_l2 = " << err_i_l2 << endl
<< "err_p_div_l2 = " << err_p_div_l2 << endl
<< "err_i_div_l2 = " << err_i_div_l2 << endl;
return (err_p_l2 <= err_p_l2_valid && err_p_div_l2 <= err_p_div_l2_valid) ? 0: 1;
}