message: Take v-dobrev's remark into account in ex1.cpp
diff: @@ -121,19 +121,23 @@ int main(int argc, char *argv[])
  // 5. Define a finite element space on the mesh. Here we use continuous
  // Lagrange finite elements of the specified order. If order < 1, we
  // instead use an isoparametric/isogeometric space.
  - FiniteElementCollection *fec = nullptr;
  + FiniteElementCollection *fec;
  + bool delete_fec;
  if (order > 0)
  {
  fec = new H1_FECollection(order, dim);
  + delete_fec = true;
  }
  else if (mesh.GetNodes())
  {
  fec = mesh.GetNodes()->OwnFEC();
  + delete_fec = false;
  cout << "Using isoparametric FEs: " << fec->Name() << endl;
  }
  else
  {
  fec = new H1_FECollection(order = 1, dim);
  + delete_fec = true;
  }
  FiniteElementSpace fespace(&mesh, fec);
  cout << "Number of finite element unknowns: "
  @@ -183,7 +187,7 @@ int main(int argc, char *argv[])
  Vector B, X;
  a.FormLinearSystem(ess_tdof_list, x, b, A, X, B);
  - cout << "Size of linear system: " << a.Height() << endl;
  + cout << "Size of linear system: " << A->Height() << endl;
  // 11. Solve the linear system A X = B.
  if (!pa)
  @@ -236,6 +240,7 @@ int main(int argc, char *argv[])
  }
  // 15. Free the used memory.
  + if (delete_fec)
  delete fec;
  return 0;
  