From 52c13be7c6958e9eedcc3e72578123669cc72795 Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 14 Aug 2008 16:03:33 +0000 Subject: [PATCH] Fuse two loops. git-svn-id: https://svn.dealii.org/trunk@16547 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-31/step-31.cc | 53 ++++++++++++++--------------- 1 file changed, 26 insertions(+), 27 deletions(-) diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index ec62724e9c..e1a0cfb6b7 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -1765,7 +1765,7 @@ void BoussinesqFlowProblem::assemble_rhs_T () { // TODO: right now, always do explicit // Euler - const bool is_first_timestep = (timestep_number == 0 ? true : true); + const bool use_bdf2_scheme = (timestep_number == 0 ? false : false); system_matrix.block(2,2) = 0; @@ -1853,41 +1853,40 @@ void BoussinesqFlowProblem::assemble_rhs_T () grad_phi_T[k] = fe_values[temperature].gradient(k,q); phi_T[k] = fe_values[temperature].value (k, q); } + + const Point p = fe_values.quadrature_point(q); + const double gamma = RightHandSide().value (p, dim+1); + for (unsigned int i=0; i old_grad_T = old_solution_grads[q][dim+1]; + for (unsigned int i=0; i old_grad_T = old_solution_grads[q][dim+1]; - Tensor<1,dim> present_u; - for (unsigned int d=0; d p = fe_values.quadrature_point(q); - - 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) - * phi_i_T) - * - fe_values.JxW(q); - } - + Tensor<1,dim> present_u; + for (unsigned int d=0; dget_dof_indices (local_dof_indices); for (unsigned int i=0; i