]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Use of new n_children function.
authorRalf Hartmann <Ralf.Hartmann@dlr.de>
Mon, 6 Mar 2006 09:46:21 +0000 (09:46 +0000)
committerRalf Hartmann <Ralf.Hartmann@dlr.de>
Mon, 6 Mar 2006 09:46:21 +0000 (09:46 +0000)
git-svn-id: https://svn.dealii.org/trunk@12538 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-12/step-12.cc
deal.II/examples/step-14/step-14.cc
deal.II/examples/step-21/step-21.cc

index 148d9701209d3b2ad148b1d378737049ddd23b70..2a84391d2020badad03f082ccdb6a7535785d2ab 100644 (file)
@@ -952,8 +952,7 @@ void DGMethod<dim>::assemble_system1 ()
                                                   // We loop over
                                                   // subfaces
                  for (unsigned int subface_no=0;
-                      subface_no<GeometryInfo<dim>::subfaces_per_face;
-                      ++subface_no)
+                      subface_no<face->n_children(); ++subface_no)
                    {
                                                       // and set the
                                                       // cell
@@ -1314,8 +1313,7 @@ void DGMethod<dim>::assemble_system2 ()
                    cell->neighbor_of_neighbor(face_no);
                  
                  for (unsigned int subface_no=0;
-                      subface_no<GeometryInfo<dim>::subfaces_per_face;
-                      ++subface_no)
+                      subface_no<face->n_children(); ++subface_no)
                    {
                      typename DoFHandler<dim>::cell_iterator neighbor_child
                         = cell->neighbor_child_on_subface (face_no, subface_no);
index ca5e7877d0a80bf96e1fd26a83ede52f64857813..ec822783c78165305420d1006912db33b7147a60 100644 (file)
@@ -3471,6 +3471,8 @@ namespace LaplaceSolver
     const unsigned int
       n_q_points = face_data.fe_face_values_cell.n_quadrature_points;
 
+    const typename DoFHandler<dim>::face_iterator
+      face = cell->face(face_no);    
     const typename DoFHandler<dim>::cell_iterator
       neighbor = cell->neighbor(face_no);    
     Assert (neighbor.state() == IteratorState::valid,
@@ -3495,8 +3497,7 @@ namespace LaplaceSolver
                                     // for one face for all the
                                     // sub-faces now:
     for (unsigned int subface_no=0;
-        subface_no<GeometryInfo<dim>::subfaces_per_face;
-        ++subface_no)
+        subface_no<face->n_children(); ++subface_no)
       {
                                         // Start with some checks
                                         // again: get an iterator
@@ -3584,10 +3585,8 @@ namespace LaplaceSolver
                                     // computed and do not carry an
                                     // invalid value.
     double sum = 0;
-    typename DoFHandler<dim>::face_iterator face = cell->face(face_no);
     for (unsigned int subface_no=0;
-        subface_no<GeometryInfo<dim>::subfaces_per_face;
-        ++subface_no) 
+        subface_no<face->n_children(); ++subface_no) 
       {
        Assert (face_integrals.find(face->child(subface_no)) !=
                face_integrals.end(),
index 9fabd1168c139ef495a8cfd187fbd44073a23dee..3700e7f1804fb71c9ea9b95f4e2d4716f23b53c1 100644 (file)
@@ -980,8 +980,7 @@ void DGMethod<dim>::assemble_system1 ()
                                                   // We loop over
                                                   // subfaces
                  for (unsigned int subface_no=0;
-                      subface_no<GeometryInfo<dim>::subfaces_per_face;
-                      ++subface_no)
+                      subface_no<face->n_children(); ++subface_no)
                    {
                                                       // and set the
                                                       // cell
@@ -1378,8 +1377,7 @@ void DGMethod<dim>::assemble_system2 ()
                    cell->neighbor_of_neighbor(face_no);
                  
                  for (unsigned int subface_no=0;
-                      subface_no<GeometryInfo<dim>::subfaces_per_face;
-                      ++subface_no)
+                      subface_no<face->n_children(); ++subface_no)
                    {
                      typename hp::DoFHandler<dim>::cell_iterator neighbor_child=
                        neighbor->child(GeometryInfo<dim>::child_cell_on_face(

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.