From dc5e1cde8f05a539869f9fda882a038d17c79732 Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 14 Aug 2009 17:57:52 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-32/step-32.cc | 2 ++ 1 file changed, 2 insertions(+) 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); -- 2.39.5