message: Cleaning up comments and style
diff: // Extruder Miniapp: Extrude a low dimensional mesh into a higher dimension
  // ------------------------------------------------------------------------
  //
  -// This miniapp performs multiple levels of adaptive mesh refinement to resolve
  -// the interfaces between different "materials" in the mesh, as specified by the
  -// given material() function. It can be used as a simple initial mesh generator,
  -// for example in the case when the interface is too complex to describe without
  -// local refinement. Both conforming and non-conforming refinements are supported.
  +// This miniapp creates higher dimensional meshes from lower dimensional meshes
  +// by extrusion. Simple coordinate transformations can also be applied if
  +// desired. The initial mesh can be 1D or 2D. 1D meshes can be extruded in the
  +// y-direction first and then in the z-direction. 2D meshes can be triangular,
  +// quadrilateral, or contain both element types. The initial mesh can also be
  +// curved although NURBS meshes are not supported.
  +//
  +// The resulting mesh will be displayed with GLVis (unless disabled) and also
  +// written to the file "extruder.mesh".
  //
  // Compile with: make extruder
  //
  // extruder -m ../../data/star.mesh -nz 3
  // extruder -m ../../data/star-mixed.mesh -nz 3
  // extruder -m ../../data/square-disc.mesh -nz 3
  -// extruder -m ../../data/square-disc.mesh -nz 3
  // extruder -m ../../data/inline-segment.mesh -ny 8 -wy 2 -trans
  // extruder -m ../../data/inline-segment.mesh -ny 8 -wy 2 -nz 12 -hz 3 -trans
  -// extruder -m ../../data/star.mesh -nz 3 -trans
  +// extruder -m ../../data/square-disc-p2.mesh -nz 16 -hz 2 -trans
  #include "mfem.hpp"
  #include <fstream>
  @@ -81,13 +84,13 @@ int main(int argc, char *argv[])
  int dim = mesh->Dimension();
  // Determine the order to use for a transformed mesh
  - if ( order < 0 && trans )
  - {
  int meshOrder = 1;
  if (mesh->GetNodalFESpace() != NULL)
  {
  - meshOrder = mesh->GetNodalFESpace()->GetORder(0);
  + meshOrder = mesh->GetNodalFESpace()->GetOrder(0);
  }
  + if (order < 0 && trans)
  + {
  order = meshOrder;
  }
  @@ -109,8 +112,11 @@ int main(int argc, char *argv[])
  mesh = mesh2d;
  dim = 2;
  if (trans)
  + {
  + if (order != meshOrder)
  {
  mesh->SetCurvature(order, false, 2, Ordering::byVDIM);
  + }
  mesh->Transform(trans2D);
  }
  newMesh = true;
  @@ -125,8 +131,11 @@ int main(int argc, char *argv[])
  mesh = mesh3d;
  dim = 3;
  if (trans)
  + {
  + if (order != meshOrder)
  {
  mesh->SetCurvature(order, false, 3, Ordering::byVDIM);
  + }
  mesh->Transform(trans3D);
  }
  newMesh = true;
  