From 9ca656ab3a4c27dc12bcce7d69eb409e0c9fee0c Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 13 Jan 2020 16:08:23 -0700 Subject: [PATCH] Correctly name the variable that denotes gamma/h. --- examples/step-71/step-71.cc | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) 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::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::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) ) * -- 2.39.5