From: Wolfgang Bangerth Date: Mon, 18 May 2020 23:54:37 +0000 (-0600) Subject: Rename the penalty parameter. X-Git-Tag: v9.2.0-rc3~1^2~4 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b3f3022344850ea9c78fe4d6e236baa1a70c2375;p=dealii.git Rename the penalty parameter. --- diff --git a/examples/step-70/step-70.cc b/examples/step-70/step-70.cc index 0c524d61d7..0f19d3f9f2 100644 --- a/examples/step-70/step-70.cc +++ b/examples/step-70/step-70.cc @@ -1325,7 +1325,8 @@ namespace Step70 fluid_fe->dofs_per_cell); dealii::Vector local_rhs(fluid_fe->dofs_per_cell); - const auto k = 1.0 / GridTools::minimal_cell_diameter(fluid_tria); + const auto penalty_parameter = + 1.0 / GridTools::minimal_cell_diameter(fluid_tria); // We loop over all the local particles. Although this could be achieved // directly by looping over all the cells, this would force us @@ -1380,11 +1381,11 @@ namespace Step70 fluid_fe->system_to_component_index(j).first; if (comp_i == comp_j) local_matrix(i, j) += - k * par.penalty_term * + penalty_parameter * par.penalty_term * fluid_fe->shape_value(i, ref_q) * fluid_fe->shape_value(j, ref_q) * JxW; } - local_rhs(i) += k * par.penalty_term * + local_rhs(i) += penalty_parameter * par.penalty_term * solid_velocity.value(real_q, comp_i) * fluid_fe->shape_value(i, ref_q) * JxW; }