From: Wolfgang Bangerth Date: Fri, 14 Aug 2009 17:57:52 +0000 (+0000) Subject: Call compress() on the rhs vector in project_initial_temperature. This X-Git-Tag: v8.0.0~7303 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=36aab63d26411e2e162bef0e6d4e2ebfbb7bcaed;p=dealii.git Call compress() on the rhs vector in project_initial_temperature. This didn't appear to matter so far because we always projected a zero function, but if the initial temperature is spatially variable we get very hard to track down errors without this call :-( git-svn-id: https://svn.dealii.org/trunk@19265 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index 0e499f3608..b284825b2a 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -1351,6 +1351,8 @@ void BoussinesqFlowProblem::project_temperature_field () rhs); } + rhs.compress (); + ReductionControl control(5*rhs.size(), 0., 1e-12, false, false); GrowingVectorMemory memory; SolverCG cg(control,memory);