From d2f1dbbd2b7a7313b91d6c67e748bc3a864e24cb Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth 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 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