message: Increasing integration order in ComputeLpError methods so they work with DG spaces
diff: @@ -1893,7 +1893,7 @@ double GridFunction::ComputeL2Error(
  fdof = fe->GetDof();
  transf = fes->GetElementTransformation(i);
  shape.SetSize(fdof);
  - intorder = 2*fe->GetOrder() + 1; // <----------
  + intorder = 2*fe->GetOrder() + 3; // <----------
  const IntegrationRule *ir;
  if (irs)
  {
  @@ -1948,7 +1948,7 @@ double GridFunction::ComputeL2Error(
  {
  if (elems != NULL && (*elems)[i] == 0) { continue; }
  fe = fes->GetFE(i);
  - int intorder = 2*fe->GetOrder() + 1; // <----------
  + int intorder = 2*fe->GetOrder() + 3; // <----------
  const IntegrationRule *ir;
  if (irs)
  {
  @@ -2143,7 +2143,7 @@ double GridFunction::ComputeMaxError(
  fdof = fe->GetDof();
  transf = fes->GetElementTransformation(i);
  shape.SetSize(fdof);
  - intorder = 2*fe->GetOrder() + 1; // <----------
  + intorder = 2*fe->GetOrder() + 3; // <----------
  const IntegrationRule *ir;
  if (irs)
  {
  @@ -2309,7 +2309,7 @@ double GridFunction::ComputeLpError(const double p, Coefficient &exsol,
  }
  else
  {
  - int intorder = 2*fe->GetOrder() + 1; // <----------
  + int intorder = 2*fe->GetOrder() + 3; // <----------
  ir = &(IntRules.Get(fe->GetGeomType(), intorder));
  }
  GetValues(i, *ir, vals);
  @@ -2378,7 +2378,7 @@ void GridFunction::ComputeElementLpErrors(const double p, Coefficient &exsol,
  }
  else
  {
  - int intorder = 2*fe->GetOrder() + 1; // <----------
  + int intorder = 2*fe->GetOrder() + 3; // <----------
  ir = &(IntRules.Get(fe->GetGeomType(), intorder));
  }
  GetValues(i, *ir, vals);
  