From: bangerth Date: Thu, 14 Aug 2008 16:10:48 +0000 (+0000) Subject: Fix computation of lhs matrix: it needs to depend on the time step, but previously... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4479f8b1f881ed36055eaec33297ae376a54b867;p=dealii-svn.git Fix computation of lhs matrix: it needs to depend on the time step, but previously did not. git-svn-id: https://svn.dealii.org/trunk@16548 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 e1a0cfb6b7..a5fa1d70d1 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -1843,8 +1843,8 @@ void BoussinesqFlowProblem::assemble_rhs_T () // of diffusion (determined // impirically) to keep the // scheme stable - const double kappa = std::max (5e-4 * cell->diameter(), - 1e-6); + const double kappa = std::max (5e-3 * cell->diameter(), + 1e-5); for (unsigned int q=0; q::assemble_rhs_T () 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]; + Assert (false, ExcNotImplemented()); + } + else + { + 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 present_u; + for (unsigned int d=0; d