From: Matthias Maier Date: Fri, 15 Jun 2018 23:55:17 +0000 (-0500) Subject: step-7: make bilinear forms more readable X-Git-Tag: v9.1.0-rc1~701^2~7 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=552d50402e149faecd1b486e20eded9e19cdd989;p=dealii.git step-7: make bilinear forms more readable --- diff --git a/examples/step-7/step-7.cc b/examples/step-7/step-7.cc index 49202e94a5..d6c30eab7b 100644 --- a/examples/step-7/step-7.cc +++ b/examples/step-7/step-7.cc @@ -655,14 +655,16 @@ namespace Step7 // The first thing that has changed is the bilinear form. It // now contains the additional term from the Helmholtz // equation: - cell_matrix(i, j) += ((fe_values.shape_grad(i, q_point) * - fe_values.shape_grad(j, q_point) + - fe_values.shape_value(i, q_point) * - fe_values.shape_value(j, q_point)) * + cell_matrix(i, j) += ((fe_values.shape_grad(i, q_point) * // + fe_values.shape_grad(j, q_point) // + + // + fe_values.shape_value(i, q_point) * // + fe_values.shape_value(j, q_point)) * // fe_values.JxW(q_point)); - cell_rhs(i) += (fe_values.shape_value(i, q_point) * - rhs_values[q_point] * fe_values.JxW(q_point)); + cell_rhs(i) += (fe_values.shape_value(i, q_point) * // + rhs_values[q_point] * // + fe_values.JxW(q_point)); } // Then there is that second term on the right hand side, the contour @@ -707,9 +709,9 @@ namespace Step7 fe_face_values.normal_vector(q_point)); for (unsigned int i = 0; i < dofs_per_cell; ++i) - cell_rhs(i) += - (neumann_value * fe_face_values.shape_value(i, q_point) * - fe_face_values.JxW(q_point)); + cell_rhs(i) += (neumann_value * // + fe_face_values.shape_value(i, q_point) * // + fe_face_values.JxW(q_point)); } }