message: Make memory allocation precise to avoid futures issues.
diff: @@ -141,18 +141,21 @@ int main(int argc, char *argv[])
  // linear forms fform and gform for the right hand side. The data
  // allocated by x and rhs are passed as a reference to the grid functions
  // (u,p) and the linear forms (fform, gform).
  - BlockVector x(block_offsets), rhs(block_offsets);
  + MemoryType mt = device.GetMemoryType();
  + BlockVector x(block_offsets, mt), rhs(block_offsets, mt);
  LinearForm *fform(new LinearForm);
  fform->Update(R_space, rhs.GetBlock(0), 0);
  fform->AddDomainIntegrator(new VectorFEDomainLFIntegrator(fcoeff));
  fform->AddBoundaryIntegrator(new VectorFEBoundaryFluxLFIntegrator(fnatcoeff));
  fform->Assemble();
  + fform->SyncAliasMemory(rhs);
  LinearForm *gform(new LinearForm);
  gform->Update(W_space, rhs.GetBlock(1), 0);
  gform->AddDomainIntegrator(new DomainLFIntegrator(gcoeff));
  gform->Assemble();
  + gform->SyncAliasMemory(rhs);
  // 9. Assemble the finite element matrices for the Darcy operator
  //
  