]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Correctly name the variable that denotes gamma/h.
authorWolfgang Bangerth <bangerth@colostate.edu>
Mon, 13 Jan 2020 23:08:23 +0000 (16:08 -0700)
committerWolfgang Bangerth <bangerth@colostate.edu>
Tue, 14 Jan 2020 00:52:45 +0000 (17:52 -0700)
examples/step-71/step-71.cc

index b8f2b44e153e5a8aff973e6476f0183a7889d4fb..36411e14bd8c2ce007d05e0296929abca39e6696 100644 (file)
@@ -504,7 +504,7 @@ namespace Step71
       // further adjusted if one were to use hanging nodes resulting from
       // adaptive mesh refinement.
       const unsigned int p = fe.degree;
-      const double       gamma =
+      const double       gamma_over_h =
         std::max((1.0 * p * (p + 1) /
                   cell->extent_in_direction(
                     GeometryInfo<dim>::unit_normal_direction[f])),
@@ -561,7 +561,7 @@ namespace Step71
                      - av_hessian_j_dot_n_dot_n      // - {grad^2 u n n }
                          * jump_grad_i_dot_n         // [grad v n]
                      +                               // +
-                     gamma *                         // gamma
+                     gamma_over_h *                  // gamma/h
                        jump_grad_i_dot_n *           // [grad v n]
                        jump_grad_j_dot_n) *          // [grad u n]
                     fe_interface_values.JxW(qpoint); // dx
@@ -611,7 +611,7 @@ namespace Step71
       // face (as we are on the boundary), the computation of the penalty
       // factor $\gamma$ is substantially simpler:
       const unsigned int p = fe.degree;
-      const double       gamma =
+      const double       gamma_over_h =
         (1.0 * p * (p + 1) /
          cell->extent_in_direction(
            GeometryInfo<dim>::unit_normal_direction[face_no]));
@@ -654,7 +654,7 @@ namespace Step71
                      - av_hessian_j_dot_n_dot_n // - {grad^2 u n n}
                          * jump_grad_i_dot_n    //   [grad v n]
                                                 //
-                     + gamma                    //  gamma
+                     + gamma_over_h             //  gamma/h
                          * jump_grad_i_dot_n    // [grad v n]
                          * jump_grad_j_dot_n    // [grad u n]
                      ) *
@@ -665,7 +665,7 @@ namespace Step71
                 (-av_hessian_i_dot_n_dot_n *       // - {grad^2 v n n }
                    (exact_gradients[qpoint] * n)   //   (grad u_exact . n)
                  +                                 // +
-                 gamma                             //  gamma
+                 gamma_over_h                      //  gamma/h
                    * jump_grad_i_dot_n             // [grad v n]
                    * (exact_gradients[qpoint] * n) // (grad u_exact . n)
                  ) *

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.