message: Slight update to vtk output of Point meshes.
  Now supports Point meshes with nodal grid functions.
diff: @@ -7269,7 +7269,12 @@ void Mesh::PrintVTK(std::ostream &out)
  }
  out << "CELLS " << NumOfElements << ' ' << size << '\n';
  const char *fec_name = Nodes->FESpace()->FEColl()->Name();
  - if (!strcmp(fec_name, "Linear") ||
  +
  + if (!strcmp(fec_name, "H1_0D_P1"))
  + {
  + order = 0;
  + }
  + else if (!strcmp(fec_name, "Linear") ||
  !strcmp(fec_name, "H1_2D_P1") ||
  !strcmp(fec_name, "H1_3D_P1"))
  {
  @@ -7291,6 +7296,12 @@ void Mesh::PrintVTK(std::ostream &out)
  {
  Nodes->FESpace()->GetElementDofs(i, dofs);
  out << dofs.Size();
  + if (order == 0)
  + {
  + MFEM_ASSERT(dofs.Size() == 1,
  + "Point meshes should have a single dof per element");
  + out << " " << dofs[0];
  + }
  if (order == 1)
  {
  for (int j = 0; j < dofs.Size(); j++)
  @@ -7325,6 +7336,13 @@ void Mesh::PrintVTK(std::ostream &out)
  for (int i = 0; i < NumOfElements; i++)
  {
  int vtk_cell_type = 5;
  + if (order == 0)
  + {
  + switch (elements[i]->GetGeometryType())
  + {
  + case Geometry::POINT: vtk_cell_type = 1; break;
  + }
  + }
  if (order == 1)
  {
  switch (elements[i]->GetGeometryType())
  