message: PA mass diagonal works for multiple elements, does not work in parallel.
diff: @@ -19,10 +19,13 @@ namespace assemblediagonalpa
  TEST_CASE("massdiag")
  {
  - std::cout << "Testing partial assembly mass diagonal: 1 element" << std::endl;
  + for (int ne = 1; ne < 3; ++ne)
  + {
  + std::cout << "Testing partial assembly mass diagonal: "
  + << ne*ne << " elements." << std::endl;
  for (int order = 1; order < 5; ++order)
  {
  - Mesh mesh(1, 1, Element::QUADRILATERAL, 1, 1.0, 1.0);
  + Mesh mesh(ne, ne, Element::QUADRILATERAL, 1, 1.0, 1.0);
  FiniteElementCollection *h1_fec = new H1_FECollection(order, 2);
  FiniteElementSpace h1_fespace(&mesh, h1_fec);
  BilinearForm paform(&h1_fespace);
  @@ -48,5 +51,6 @@ TEST_CASE("massdiag")
  delete h1_fec;
  }
  }
  +}
  } // namespace assemblediagonalpa
  