message: Rewrite ParNCL2FaceRestriction::ComputeGatherIndices.
diff: @@ -1099,6 +1099,7 @@ void ParNCL2FaceRestriction::ComputeGatherIndices(
  const FaceType type)
  {
  Mesh &mesh = *fes.GetMesh();
  +
  // Computation of gather_indices
  int f_ind = 0;
  for (int f = 0; f < mesh.GetNumFacesWithGhost(); ++f)
  @@ -1106,11 +1107,6 @@ void ParNCL2FaceRestriction::ComputeGatherIndices(
  Mesh::FaceInformation face = mesh.GetFaceInformation(f);
  if ( face.IsOfFaceType(type) )
  {
  - // I should probably shift by nfdofs on ghost faces
  - if (m==L2FaceValues::DoubleValued && face.IsGhost())
  - {
  - /* TODO: code */
  - }
  SetFaceDofsGatherIndices1(face,f_ind);
  if (m==L2FaceValues::DoubleValued &&
  type==FaceType::Interior &&
  @@ -1122,6 +1118,9 @@ void ParNCL2FaceRestriction::ComputeGatherIndices(
  }
  else
  {
  + // Contrary to the conforming case, there is no need to call
  + // PermuteFaceL2, the permutation is achieved by the
  + // interpolation operator for simplicity.
  SetFaceDofsGatherIndices2(face,f_ind);
  }
  }
  