From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Mon, 13 Jan 2020 23:08:23 +0000 (-0700)
Subject: Correctly name the variable that denotes gamma/h.
X-Git-Tag: v9.2.0-rc1~678^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9ca656ab3a4c27dc12bcce7d69eb409e0c9fee0c;p=dealii.git

Correctly name the variable that denotes gamma/h.
---

diff --git a/examples/step-71/step-71.cc b/examples/step-71/step-71.cc
index b8f2b44e15..36411e14bd 100644
--- a/examples/step-71/step-71.cc
+++ b/examples/step-71/step-71.cc
@@ -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)
                  ) *