From: Wolfgang Bangerth Date: Thu, 29 Mar 2018 03:40:49 +0000 (-0600) Subject: Adapt step-12 some further. X-Git-Tag: v9.0.0-rc1~269^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F6117%2Fhead;p=dealii.git 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. --- 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