message: Same weird bug somewhere else.
diff: @@ -638,12 +638,13 @@ void FABilinearFormExtension::Assemble()
  const int elem_dofs = fes.GetFE(0)->GetDof();
  auto h_E = elem_begin.HostWrite();
  auto h_I = mat.HostWriteI();
  + auto h_nnz = elem_nnz.Read();
  size_t cpt = 0;
  //TODO take into account vd
  for (int e = 0; e < ne; e++)
  {
  h_E[e] = cpt + elem_dofs;
  - const int row_nnz = elem_nnz[e];
  + const int row_nnz = h_nnz[e];
  for (int dof = 0; dof < elem_dofs; dof++)
  {
  const int row = e*elem_dofs+dof;
  