message: Modify ex1 to work Element Assembly.
diff: @@ -65,6 +65,7 @@ int main(int argc, char *argv[])
  int order = 1;
  bool static_cond = false;
  bool pa = false;
  + bool ea = false;
  const char *device_config = "cpu";
  bool visualization = true;
  @@ -78,6 +79,8 @@ int main(int argc, char *argv[])
  "--no-static-condensation", "Enable static condensation.");
  args.AddOption(&pa, "-pa", "--partial-assembly", "-no-pa",
  "--no-partial-assembly", "Enable Partial Assembly.");
  + args.AddOption(&ea, "-ea", "--element-assembly", "-no-ea",
  + "--no-element-assembly", "Enable Element Assembly.");
  args.AddOption(&device_config, "-d", "--device",
  "Device configuration string, see Device::Configure().");
  args.AddOption(&visualization, "-vis", "--visualization", "-no-vis",
  @@ -167,6 +170,7 @@ int main(int argc, char *argv[])
  // domain integrator.
  BilinearForm *a = new BilinearForm(fespace);
  if (pa) { a->SetAssemblyLevel(AssemblyLevel::PARTIAL); }
  + else if (ea) { a->SetAssemblyLevel(AssemblyLevel::ELEMENT); }
  a->AddDomainIntegrator(new DiffusionIntegrator(one));
  // 10. Assemble the bilinear form and the corresponding linear system,
  @@ -183,7 +187,7 @@ int main(int argc, char *argv[])
  cout << "Size of linear system: " << A->Height() << endl;
  // 11. Solve the linear system A X = B.
  - if (!pa)
  + if (!pa && !ea)
  {
  #ifndef MFEM_USE_SUITESPARSE
  // Use a simple symmetric Gauss-Seidel preconditioner with PCG.
  