From d2f1dbbd2b7a7313b91d6c67e748bc3a864e24cb Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Wed, 28 Mar 2018 21:40:49 -0600
Subject: [PATCH] Adapt step-12 some further.

In #5988, I added a function beta(...) that computes the advection direction
in step-12, but I forgot to update one other place where we computed this
direction as well. Use the new function in there as well.
---
 examples/step-12/step-12.cc | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/examples/step-12/step-12.cc b/examples/step-12/step-12.cc
index e82e3f9beb..41cd7213d4 100644
--- a/examples/step-12/step-12.cc
+++ b/examples/step-12/step-12.cc
@@ -351,16 +351,13 @@ namespace Step12
     // the current point.
     for (unsigned int point=0; point<fe_values.n_quadrature_points; ++point)
       {
-        Point<dim> beta;
-        beta(0) = -fe_values.quadrature_point(point)(1);
-        beta(1) = fe_values.quadrature_point(point)(0);
-        beta /= beta.norm();
+        const Tensor<1,dim> beta_at_q_point = beta (fe_values.quadrature_point(point));
 
         // We solve a homogeneous equation, thus no right hand side shows up
         // in the cell term.  What's left is integrating the matrix entries.
         for (unsigned int i=0; i<fe_values.dofs_per_cell; ++i)
           for (unsigned int j=0; j<fe_values.dofs_per_cell; ++j)
-            local_matrix(i,j) += -beta *
+            local_matrix(i,j) += -beta_at_q_point *
                                  fe_values.shape_grad(i,point) *
                                  fe_values.shape_value(j,point) *
                                  JxW[point];
-- 
2.39.5