message: Removing GetIntPoint/SetIntPoint calls
diff: @@ -193,12 +193,7 @@ void GradientGridFunctionCoefficient::SetGridFunction(GridFunction *gf)
  void GradientGridFunctionCoefficient::Eval(Vector &V, ElementTransformation &T,
  const IntegrationPoint &ip)
  {
  - // Perhaps we should avoid altering T so cache its IntegrationPoint
  - const IntegrationPoint & ipT = T.GetIntPoint();
  - T.SetIntPoint(&ip);
  GridFunc->GetGradient(T, V);
  - // Reset T's IntegrationPoint
  - T.SetIntPoint(&ipT);
  }
  void GradientGridFunctionCoefficient::Eval(
  @@ -224,12 +219,7 @@ void CurlGridFunctionCoefficient::SetGridFunction(GridFunction *gf)
  void CurlGridFunctionCoefficient::Eval(Vector &V, ElementTransformation &T,
  const IntegrationPoint &ip)
  {
  - // Perhaps we should avoid altering T so cache its IntegrationPoint
  - const IntegrationPoint & ipT = T.GetIntPoint();
  - T.SetIntPoint(&ip);
  GridFunc->GetCurl(T, V);
  - // Reset T's IntegrationPoint
  - T.SetIntPoint(&ipT);
  }
  DivergenceGridFunctionCoefficient::DivergenceGridFunctionCoefficient (
  @@ -241,13 +231,7 @@ DivergenceGridFunctionCoefficient::DivergenceGridFunctionCoefficient (
  double DivergenceGridFunctionCoefficient::Eval(ElementTransformation &T,
  const IntegrationPoint &ip)
  {
  - // Perhaps we should avoid altering T so cache its IntegrationPoint
  - const IntegrationPoint & ipT = T.GetIntPoint();
  - T.SetIntPoint(&ip);
  - double v = GridFunc->GetDivergence(T);
  - // Reset T's IntegrationPoint
  - T.SetIntPoint(&ipT);
  - return v;
  + return GridFunc->GetDivergence(T);
  }
  void VectorDeltaCoefficient::SetDirection(const Vector &_d)
  