message: Adding FiniteElement::INTEGRAL support to GridFunction::GetValues
diff: @@ -462,17 +462,28 @@ const
  fes->GetElementDofs(i, dofs);
  fes->DofsToVDofs(vdim-1, dofs);
  const FiniteElement *FElem = fes->GetFE(i);
  - MFEM_ASSERT(FElem->GetMapType() == FiniteElement::VALUE,
  - "invalid FE map type");
  int dof = FElem->GetDof();
  Vector DofVal(dof), loc_data(dof);
  GetSubVector(dofs, loc_data);
  + if (FElem->GetMapType() == FiniteElement::VALUE)
  + {
  for (int k = 0; k < n; k++)
  {
  FElem->CalcShape(ir.IntPoint(k), DofVal);
  vals(k) = DofVal * loc_data;
  }
  }
  + else
  + {
  + ElementTransformation *Tr = fes->GetElementTransformation(i);
  + for (int k = 0; k < n; k++)
  + {
  + Tr->SetIntPoint(&ir.IntPoint(k));
  + FElem->CalcPhysShape(*Tr, DofVal);
  + vals(k) = DofVal * loc_data;
  + }
  + }
  +}
  void GridFunction::GetValues(int i, const IntegrationRule &ir, Vector &vals,
  DenseMatrix &tr, int vdim)
  