]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Be uniform in the order of things we start the assembly loop with. 16955/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Thu, 2 May 2024 18:09:49 +0000 (23:39 +0530)
committerWolfgang Bangerth <bangerth@colostate.edu>
Thu, 2 May 2024 18:09:49 +0000 (23:39 +0530)
examples/step-20/step-20.cc
examples/step-22/step-22.cc
examples/step-4/step-4.cc
examples/step-40/step-40.cc
examples/step-5/step-5.cc
examples/step-6/step-6.cc
examples/step-7/step-7.cc
examples/step-8/step-8.cc

index d90e8d493dd001dc738dc41ecc441d2bd3399470..63285a2705c0c2d822212f90d4b8e5eb8c6ac890 100644 (file)
@@ -488,6 +488,7 @@ namespace Step20
     for (const auto &cell : dof_handler.active_cell_iterators())
       {
         fe_values.reinit(cell);
+
         local_matrix = 0;
         local_rhs    = 0;
 
index 453f322966d544bfb7292837460309c06aa2fec9..bb3b712056e6eeed6854a3a7eb84371e89f8871c 100644 (file)
@@ -659,6 +659,7 @@ namespace Step22
     for (const auto &cell : dof_handler.active_cell_iterators())
       {
         fe_values.reinit(cell);
+
         local_matrix                = 0;
         local_preconditioner_matrix = 0;
         local_rhs                   = 0;
index d0e622d7c6d86a0fa24dd48c7620b0b3c1c196c7..4e3f6c9a24d898427919bb755ca698d6f62ddc37 100644 (file)
@@ -354,6 +354,7 @@ void Step4<dim>::assemble_system()
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
       fe_values.reinit(cell);
+
       cell_matrix = 0;
       cell_rhs    = 0;
 
index 622700234bb94b3d3e8600a3e9f1f117ed3e314f..a0d3e51b6758aee836997d2fd91825af8c6f7aad 100644 (file)
@@ -399,11 +399,11 @@ namespace Step40
     for (const auto &cell : dof_handler.active_cell_iterators())
       if (cell->is_locally_owned())
         {
+          fe_values.reinit(cell);
+
           cell_matrix = 0.;
           cell_rhs    = 0.;
 
-          fe_values.reinit(cell);
-
           for (unsigned int q_point = 0; q_point < n_q_points; ++q_point)
             {
               const double rhs_value =
index bba78459f6ab6bec91f2e998a14587e9edbfe797..394e83629dda765aa39840e5816d06a14555e8a4 100644 (file)
@@ -178,11 +178,11 @@ void Step5<dim>::assemble_system()
   // coefficient value at each quadrature point.
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
+      fe_values.reinit(cell);
+
       cell_matrix = 0.;
       cell_rhs    = 0.;
 
-      fe_values.reinit(cell);
-
       for (const unsigned int q_index : fe_values.quadrature_point_indices())
         {
           const double current_coefficient =
index 3dd3c7153c91b90c70fb68a3bb4f741842337350..9522232af0299f23c642531bdfbb5b6f7f1963e2 100644 (file)
@@ -271,11 +271,11 @@ void Step6<dim>::assemble_system()
 
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
+      fe_values.reinit(cell);
+
       cell_matrix = 0;
       cell_rhs    = 0;
 
-      fe_values.reinit(cell);
-
       for (const unsigned int q_index : fe_values.quadrature_point_indices())
         {
           const double current_coefficient =
index b13c306c5cf4fca40c5d486c5d76d5c107207a1f..5be9442e223a5352273347054bc5ce5f5a3ae58c 100644 (file)
@@ -583,11 +583,11 @@ namespace Step7
     // previous examples, so we only comment on the things that have changed.
     for (const auto &cell : dof_handler.active_cell_iterators())
       {
+        fe_values.reinit(cell);
+
         cell_matrix = 0.;
         cell_rhs    = 0.;
 
-        fe_values.reinit(cell);
-
         right_hand_side.value_list(fe_values.get_quadrature_points(),
                                    rhs_values);
 
index 517f37bda5a4d34220622c2fd7e4c844461219ff..f9c29a94806b720ce247407cca16733577d69634 100644 (file)
@@ -318,11 +318,11 @@ namespace Step8
     // Now we can begin with the loop over all cells:
     for (const auto &cell : dof_handler.active_cell_iterators())
       {
+        fe_values.reinit(cell);
+
         cell_matrix = 0;
         cell_rhs    = 0;
 
-        fe_values.reinit(cell);
-
         // Next we get the values of the coefficients at the quadrature
         // points. Likewise for the right hand side:
         lambda.value_list(fe_values.get_quadrature_points(), lambda_values);

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.