message: Tweak the implementation of Mesh::PrintVTK.
diff: @@ -7265,22 +7265,23 @@ void Mesh::PrintVTK(std::ostream &out)
  for (int i = 0; i < NumOfElements; i++)
  {
  Nodes->FESpace()->GetElementDofs(i, dofs);
  + MFEM_ASSERT(Dim != 0 || dofs.Size() == 1,
  + "Point meshes should have a single dof per element");
  size += dofs.Size() + 1;
  }
  out << "CELLS " << NumOfElements << ' ' << size << '\n';
  const char *fec_name = Nodes->FESpace()->FEColl()->Name();
  - if (!strcmp(fec_name, "H1_0D_P1"))
  - {
  - order = 0;
  - }
  - else if (!strcmp(fec_name, "Linear") ||
  + if (!strcmp(fec_name, "Linear") ||
  + !strcmp(fec_name, "H1_0D_P1") ||
  + !strcmp(fec_name, "H1_1D_P1") ||
  !strcmp(fec_name, "H1_2D_P1") ||
  !strcmp(fec_name, "H1_3D_P1"))
  {
  order = 1;
  }
  else if (!strcmp(fec_name, "Quadratic") ||
  + !strcmp(fec_name, "H1_1D_P2") ||
  !strcmp(fec_name, "H1_2D_P2") ||
  !strcmp(fec_name, "H1_3D_P2"))
  {
  @@ -7296,12 +7297,6 @@ 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++)
  @@ -7314,6 +7309,7 @@ void Mesh::PrintVTK(std::ostream &out)
  const int *vtk_mfem;
  switch (elements[i]->GetGeometryType())
  {
  + case Geometry::SEGMENT:
  case Geometry::TRIANGLE:
  case Geometry::SQUARE:
  vtk_mfem = vtk_quadratic_hex; break; // identity map
  @@ -7336,13 +7332,6 @@ 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())
  