message: Increase default quadrature in vector linear forms
  To ensure h^{p+1} convergence, the quadrature used for the
  right-hand side needs to be sufficiently accurate. This makes the
  vector versions of LinearFormIntegrator use the same degree of
  exactness as the scalar versions.
diff: @@ -181,7 +181,7 @@ void VectorDomainLFIntegrator::AssembleRHSElementVect(
  const IntegrationRule *ir = IntRule;
  if (ir == NULL)
  {
  - int intorder = el.GetOrder() + 1;
  + int intorder = 2*el.GetOrder();
  ir = &IntRules.Get(el.GetGeomType(), intorder);
  }
  @@ -240,7 +240,7 @@ void VectorBoundaryLFIntegrator::AssembleRHSElementVect(
  const IntegrationRule *ir = IntRule;
  if (ir == NULL)
  {
  - int intorder = el.GetOrder() + 1;
  + int intorder = 2*el.GetOrder();
  ir = &IntRules.Get(el.GetGeomType(), intorder);
  }
  @@ -275,7 +275,7 @@ void VectorBoundaryLFIntegrator::AssembleRHSElementVect(
  const IntegrationRule *ir = IntRule;
  if (ir == NULL)
  {
  - int intorder = el.GetOrder() + 1;
  + int intorder = 2*el.GetOrder();
  ir = &IntRules.Get(Tr.FaceGeom, intorder);
  }
  