From: bangerth Date: Wed, 13 Aug 2008 20:03:08 +0000 (+0000) Subject: Don't integrate by parts the terms on the temperature rhs. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=fbab6076e873f13428cb39c7ed98a9e7ebf598bf;p=dealii-svn.git Don't integrate by parts the terms on the temperature rhs. git-svn-id: https://svn.dealii.org/trunk@16536 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index cd1118a02b..2e7a91ea4f 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -1801,7 +1801,7 @@ void BoussinesqFlowProblem::assemble_rhs_T () n_face_q_points, Vector(dim+2)); - std::vector > > present_solution_grads( + std::vector > > old_solution_grads( n_q_points, std::vector >(dim+2)); @@ -1829,35 +1829,26 @@ void BoussinesqFlowProblem::assemble_rhs_T () fe_values.reinit (cell); fe_values.get_function_values (old_solution, old_solution_values); + fe_values.get_function_gradients (old_solution, old_solution_grads); fe_values.get_function_values (solution, present_solution_values); - fe_values.get_function_gradients (solution, present_solution_grads); for (unsigned int q=0; q old_grad_T = old_solution_grads[q][dim+1]; + Tensor<1,dim> present_u; for (unsigned int d=0; d grad_phi_i_T = fe_values[temperature].gradient (i, q); - + const double phi_i_T = fe_values[temperature].value (i, q); const Point p = fe_values.quadrature_point(q); - local_rhs(i) += (time_step * - old_T * - (present_u * - grad_phi_i_T - + - present_div_u * - phi_i_T) - + - old_T * phi_i_T + local_rhs(i) += (old_T * phi_i_T + - + time_step * + present_u * old_grad_T * phi_i_T + time_step * RightHandSide().value (p, dim+1)