]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Call compress() on the rhs vector in project_initial_temperature. This
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 14 Aug 2009 17:57:52 +0000 (17:57 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 14 Aug 2009 17:57:52 +0000 (17:57 +0000)
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

index 0e499f3608a3603146fe44d9bb69e0e2a4f2c545..b284825b2af8fe21c338468406c74c374ac8e419 100644 (file)
@@ -1351,6 +1351,8 @@ void BoussinesqFlowProblem<dim>::project_temperature_field ()
                                                            rhs);
       }
 
+  rhs.compress ();
+
   ReductionControl  control(5*rhs.size(), 0., 1e-12, false, false);
   GrowingVectorMemory<TrilinosWrappers::MPI::Vector> memory;
   SolverCG<TrilinosWrappers::MPI::Vector> cg(control,memory);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.