message: Cleaning up ex1p.
diff: @@ -53,11 +53,8 @@ int main(int argc, char *argv[])
  // 2. Parse command-line options.
  const char *mesh_file = "../data/star.mesh";
  int order = 1;
  - int ref_levels = 0;
  bool static_cond = false;
  bool visualization = 1;
  - int seed = 0;
  - bool ptest = false, aniso = false;
  OptionsParser args(argc, argv);
  args.AddOption(&mesh_file, "-m", "--mesh",
  @@ -65,18 +62,11 @@ int main(int argc, char *argv[])
  args.AddOption(&order, "-o", "--order",
  "Finite element order (polynomial degree) or -1 for"
  " isoparametric space.");
  - args.AddOption(&ref_levels, "-r", "--ref-levels",
  - "Serial refinement levels.");
  - args.AddOption(&seed, "-s", "--seed", "Random seed.");
  args.AddOption(&static_cond, "-sc", "--static-condensation", "-no-sc",
  "--no-static-condensation", "Enable static condensation.");
  args.AddOption(&visualization, "-vis", "--visualization", "-no-vis",
  "--no-visualization",
  "Enable or disable GLVis visualization.");
  - args.AddOption(&ptest, "-pt", "--ptest", "-npt", "--no-ptest",
  - "Run P matrix test only.");
  - args.AddOption(&aniso, "-a", "--aniso", "-i", "--iso",
  - "Enable anisotropic refinement.");
  args.Parse();
  if (!args.Good())
  {
  @@ -98,18 +88,16 @@ int main(int argc, char *argv[])
  Mesh *mesh = new Mesh(mesh_file, 1, 1);
  int dim = mesh->Dimension();
  - mesh->EnsureNCMesh();
  -
  // 4. Refine the serial mesh on all processors to increase the resolution. In
  // this example we do 'ref_levels' of uniform refinement. We choose
  // 'ref_levels' to be the largest number that gives a final mesh with no
  // more than 10,000 elements.
  {
  - srand(seed);
  + int ref_levels =
  + (int)floor(log(10000./mesh->GetNE())/log(2.)/dim);
  for (int l = 0; l < ref_levels; l++)
  {
  - //mesh->UniformRefinement();
  - mesh->RandomRefinement(0.5, false);
  + mesh->UniformRefinement();
  }
  }
  @@ -119,7 +107,7 @@ int main(int argc, char *argv[])
  ParMesh *pmesh = new ParMesh(MPI_COMM_WORLD, *mesh);
  delete mesh;
  {
  - int par_ref_levels = 1;
  + int par_ref_levels = 2;
  for (int l = 0; l < par_ref_levels; l++)
  {
  pmesh->UniformRefinement();
  @@ -153,12 +141,6 @@ int main(int argc, char *argv[])
  cout << "Number of finite element unknowns: " << size << endl;
  }
  - if (ptest)
  - {
  - MPI_Finalize();
  - return 0;
  - }
  -
  // 7. Determine the list of true (i.e. parallel conforming) essential
  // boundary dofs. In this example, the boundary conditions are defined
  // by marking all the boundary attributes from the mesh as essential
  