message: Adding support for mixed FiniteElement::INTEGRAL/VALUE in NodalFiniteElement::Project
diff: @@ -722,11 +722,11 @@ void NodalFiniteElement::Project(
  {
  if (fe.GetRangeType() == SCALAR)
  {
  - MFEM_ASSERT(map_type == fe.GetMapType(), "");
  -
  Vector shape(fe.GetDof());
  I.SetSize(dof, fe.GetDof());
  + if (map_type == fe.GetMapType())
  + {
  for (int k = 0; k < dof; k++)
  {
  fe.CalcShape(Nodes.IntPoint(k), shape);
  @@ -737,6 +737,23 @@ void NodalFiniteElement::Project(
  }
  }
  else
  + {
  + for (int k = 0; k < dof; k++)
  + {
  + Trans.SetIntPoint(&Nodes.IntPoint(k));
  + fe.CalcPhysShape(Trans, shape);
  + if (map_type == INTEGRAL)
  + {
  + shape *= Trans.Weight();
  + }
  + for (int j = 0; j < shape.Size(); j++)
  + {
  + I(k,j) = (fabs(shape(j)) < 1e-12) ? 0.0 : shape(j);
  + }
  + }
  + }
  + }
  + else
  {
  DenseMatrix vshape(fe.GetDof(), Trans.GetSpaceDim());
  