message: Fixed a memory leak.
  Found another bug in the DescentSolver and put a TODO there.
diff: @@ -326,7 +326,7 @@ double TMOPNewtonSolver::ComputeScalingFactor(const Vector &x,
  if (serial)
  {
  const SparseMatrix *cP = fes->GetConformingProlongation();
  - if (!cP) {x_out_loc.SetData(x_out.GetData());}
  + if (!cP) { x_out_loc = x_out; }
  else { cP->Mult(x_out, x_out_loc); }
  }
  #ifdef MFEM_USE_MPI
  @@ -507,6 +507,8 @@ double TMOPDescentNewtonSolver::ComputeScalingFactor(const Vector &x,
  for (int i = 0; i < NE; i++)
  {
  fes->GetElementVDofs(i, xdofs);
  + // TODO x_loc doesn't have valid values here!
  + MFEM_ABORT("This function has to be fixed!");
  x_loc.GetSubVector(xdofs, posV);
  for (int j = 0; j < nsp; j++)
  @@ -539,7 +541,7 @@ double TMOPDescentNewtonSolver::ComputeScalingFactor(const Vector &x,
  if (serial)
  {
  const SparseMatrix *cP = fes->GetConformingProlongation();
  - if (!cP) {x_loc.SetData(x_out.GetData());}
  + if (!cP) { x_loc = x_out; }
  else { cP->Mult(x_out,x_loc); }
  energy_out = nlf->GetGridFunctionEnergy(x_loc);
  }
  