From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Fri, 3 Apr 2015 15:13:50 +0000 (-0500)
Subject: Minor edits to a bilinear form.
X-Git-Tag: v8.3.0-rc1~322^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4ec6ed4c6f63b26e324ab218098e0dfbc9c38e35;p=dealii.git

Minor edits to a bilinear form.

The existing code is a bit hard to read because it groups terms in ways that
do not match the way it is shown in the documentation. Change the grouping so
that documentation and code match 1:1.
---

diff --git a/examples/step-15/step-15.cc b/examples/step-15/step-15.cc
index ba2c2e9b0f..d76f131925 100644
--- a/examples/step-15/step-15.cc
+++ b/examples/step-15/step-15.cc
@@ -293,17 +293,17 @@ namespace Step15
               {
                 for (unsigned int j=0; j<dofs_per_cell; ++j)
                   {
-                    cell_matrix(i, j) += (fe_values.shape_grad(i, q_point)
-                                          * coeff
-                                          * (fe_values.shape_grad(j, q_point)
-                                             -
-                                             coeff * coeff
+                    cell_matrix(i, j) += ( ((fe_values.shape_grad(i, q_point)
+                                             * coeff
+                                             * fe_values.shape_grad(j, q_point))
+                                            -
+                                            (fe_values.shape_grad(i, q_point)
+                                             * coeff * coeff * coeff
                                              * (fe_values.shape_grad(j, q_point)
                                                 *
                                                 old_solution_gradients[q_point])
-                                             * old_solution_gradients[q_point]
-                                            )
-                                          * fe_values.JxW(q_point));
+                                             * old_solution_gradients[q_point]))
+                                           * fe_values.JxW(q_point));
                   }
 
                 cell_rhs(i) -= (fe_values.shape_grad(i, q_point)