As initial debug approach, I used the function local_assemble_system instead of assemble_cell_terms, which worked (left the other functions alone). Then I tried to assemble just the boundary terms, by adding local_assemble_boundaries, uncommenting assemble_Dirichlet_boundary_terms() and assemble_Neumann_boundary_terms() and adding MeshWorker::AssembleFlags::assemble_boundary_faces to the flags in the mesh_loop. Here it failed, but I was not able to find out why yet.
-- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.