message: Adding support for map type INTEGRAL in GridFunction::GetNodalValues
diff: @@ -416,9 +416,9 @@ void GridFunction::GetNodalValues(int i, Array<double> &nval, int vdim) const
  if (FElem->GetRangeType() == FiniteElement::SCALAR)
  {
  - MFEM_ASSERT(FElem->GetMapType() == FiniteElement::VALUE,
  - "invalid FE map type");
  Vector shape(dof);
  + if (FElem->GetMapType() == FiniteElement::VALUE)
  + {
  for (k = 0; k < n; k++)
  {
  FElem->CalcShape(ElemVert->IntPoint(k), shape);
  @@ -426,6 +426,17 @@ void GridFunction::GetNodalValues(int i, Array<double> &nval, int vdim) const
  }
  }
  else
  + {
  + ElementTransformation *Tr = fes->GetElementTransformation(i);
  + for (k = 0; k < n; k++)
  + {
  + Tr->SetIntPoint(&ElemVert->IntPoint(k));
  + FElem->CalcPhysShape(*Tr, shape);
  + nval[k] = shape * ((const double *)loc_data + dof * vdim);
  + }
  + }
  + }
  + else
  {
  ElementTransformation *Tr = fes->GetElementTransformation(i);
  DenseMatrix vshape(dof, FElem->GetDim());
  