]> https://gitweb.dealii.org/ - dealii.git/commitdiff
meshworker: only assemble boundary terms for own cells
authorTimo Heister <timo.heister@gmail.com>
Fri, 3 Jan 2014 13:47:57 +0000 (13:47 +0000)
committerTimo Heister <timo.heister@gmail.com>
Fri, 3 Jan 2014 13:47:57 +0000 (13:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@32152 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/meshworker/loop.h

index 8b362bdf7f206a4dc8c319964313a09bc5412c80..22110233a2c0d1c9891742a42f1d26ca7c31bc8c 100644 (file)
@@ -220,7 +220,8 @@ namespace MeshWorker
           typename ITERATOR::AccessorType::Container::face_iterator face = cell->face(face_no);
           if (cell->at_boundary(face_no))
             {
-              if (integrate_boundary)
+              // only integrate boundary faces of own cells
+              if (integrate_boundary && own_cell)
                 {
                   dof_info.interior_face_available[face_no] = true;
                   dof_info.interior[face_no].reinit(cell, face, face_no);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.