From 0215d83181708c9f8197f236cbb5c37dff5bf6d1 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Fri, 9 Oct 2009 23:43:05 +0000 Subject: [PATCH] All header files should include base/config.h first. Here, we need a couple other includes files in addition as well. git-svn-id: https://svn.dealii.org/trunk@19793 0785d39b-7218-0410-832d-ea1e28bc413d --- .../include/numerics/mesh_worker_loop.h | 23 +++++++++++-------- 1 file changed, 14 insertions(+), 9 deletions(-) diff --git a/deal.II/deal.II/include/numerics/mesh_worker_loop.h b/deal.II/deal.II/include/numerics/mesh_worker_loop.h index c10f4990a4..1ab93d548f 100644 --- a/deal.II/deal.II/include/numerics/mesh_worker_loop.h +++ b/deal.II/deal.II/include/numerics/mesh_worker_loop.h @@ -13,6 +13,11 @@ #ifndef __deal2__mesh_worker_loop_h #define __deal2__mesh_worker_loop_h +#include +#include +#include + + DEAL_II_NAMESPACE_OPEN template class TriaActiveIterator; @@ -36,7 +41,7 @@ namespace MeshWorker return true; } } - + /** * The main work function of this namespace. Its action consists of two @@ -65,7 +70,7 @@ namespace MeshWorker { const_cast&>( begin->get_triangulation()).clear_user_flags(); - + // Loop over all cells for (ITERATOR cell = begin; cell != end; ++cell) { @@ -77,14 +82,14 @@ namespace MeshWorker cellinfo.reinit(cell); localworker.cell(cellinfo); } - + if (localworker.interior_fluxes || localworker.boundary_fluxes) for (unsigned int face_no=0; face_no < GeometryInfo::faces_per_cell; ++face_no) { typename ITERATOR::AccessorType::Container::face_iterator face = cell->face(face_no); // Treat every face only once if (face->user_flag_set ()) continue; - + if (cell->at_boundary(face_no)) { if (localworker.boundary_fluxes) @@ -100,7 +105,7 @@ namespace MeshWorker // Interior face typename ITERATOR::AccessorType::Container::cell_iterator neighbor = cell->neighbor(face_no); - + // Deal with // refinement edges // from the refined @@ -115,7 +120,7 @@ namespace MeshWorker { Assert(!cell->has_children(), ExcInternalError()); Assert(!neighbor->has_children(), ExcInternalError()); - + std::pair neighbor_face_no = cell->neighbor_of_coarser_neighbor(face_no); typename ITERATOR::AccessorType::Container::face_iterator nface @@ -142,18 +147,18 @@ namespace MeshWorker // internal face. if (internal::is_active_iterator(cell) && neighbor->has_children()) continue; - + unsigned int neighbor_face_no = cell->neighbor_of_neighbor(face_no); Assert (neighbor->face(neighbor_face_no) == face, ExcInternalError()); // Regular interior face faceinfo.reinit(cell, face, face_no); ngbrinfo.reinit(neighbor, neighbor->face(neighbor_face_no), neighbor_face_no); - + localworker.face(faceinfo, ngbrinfo); neighbor->face(neighbor_face_no)->set_user_flag (); } - + } } // faces // Execute this, if faces -- 2.39.5