message: Adding comments to skin_mesh
diff: @@ -120,6 +120,7 @@ Mesh *read_par_mesh(int np, const char *mesh_prefix)
  Mesh *skin_mesh(Mesh *mesh)
  {
  + // Determine mapping from vertex to boundary vertex
  Array<int> v2v(mesh->GetNV());
  v2v = -1;
  for (int i = 0; i < mesh->GetNBE(); i++)
  @@ -141,9 +142,11 @@ Mesh *skin_mesh(Mesh *mesh)
  }
  }
  + // Create a new mesh for the boundary
  Mesh * bmesh = new Mesh(mesh->Dimension() - 1, nbvt, mesh->GetNBE(),
  0, mesh->SpaceDimension());
  + // Copy vertices to the boundary mesh
  nbvt = 0;
  for (int i = 0; i < v2v.Size(); i++)
  {
  @@ -155,6 +158,7 @@ Mesh *skin_mesh(Mesh *mesh)
  }
  }
  + // Copy elements to the boundary mesh
  int bv[4];
  for (int i = 0; i < mesh->GetNBE(); i++)
  {
  @@ -185,10 +189,13 @@ Mesh *skin_mesh(Mesh *mesh)
  }
  bmesh->FinalizeTopology();
  + // Copy GridFunction describing nodes if present
  if (mesh->GetNodes())
  {
  FiniteElementSpace *fes = mesh->GetNodes()->FESpace();
  const FiniteElementCollection *fec = fes->FEColl();
  + if (dynamic_cast<const H1_FECollection*>(fec))
  + {
  FiniteElementCollection *fec_copy =
  FiniteElementCollection::New(fec->Name());
  FiniteElementSpace *fes_copy =
  @@ -209,7 +216,11 @@ Mesh *skin_mesh(Mesh *mesh)
  fes_copy->GetElementVDofs(i, bvdofs);
  bdr_nodes->SetSubVector(bvdofs, v);
  }
  -
  + }
  + else
  + {
  + cout << "\nDiscontinuous nodes not yet supported" << endl;
  + }
  }
  return bmesh;
  