]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove the project_back_temperature function as we don't need it here (unlike for...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 12 Oct 2007 18:02:17 +0000 (18:02 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 12 Oct 2007 18:02:17 +0000 (18:02 +0000)
git-svn-id: https://svn.dealii.org/trunk@15307 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-22/step-22.cc

index c22d20a1def5b9d8b30f6e79d67cb2a9d9164729..a02291d0f258753526c00f3fa5744291bdeb2559 100644 (file)
@@ -65,7 +65,6 @@ class BoussinesqFlowProblem
     void assemble_rhs_T ();
     double get_maximal_velocity () const;
     void solve ();
-    void project_back_temperature ();
     void output_results () const;
     
     const unsigned int   degree;
@@ -357,7 +356,7 @@ BoussinesqFlowProblem<dim>::BoussinesqFlowProblem (const unsigned int degree)
                     FE_Q<dim>(degree), 1,
                     FE_Q<dim>(degree), 1),
                 dof_handler (triangulation),
-                n_refinement_steps (5),
+                n_refinement_steps (6),
                 time_step (0)
 {}
 
@@ -781,8 +780,6 @@ void BoussinesqFlowProblem<dim>::solve ()
       }
        
                 
-    project_back_temperature ();
-        
     std::cout << "   "
               << solver_control.last_step()
               << " CG iterations for temperature."
@@ -798,7 +795,7 @@ void BoussinesqFlowProblem<dim>::solve ()
 template <int dim>
 void BoussinesqFlowProblem<dim>::output_results ()  const
 {
-  if (timestep_number % 10 != 0)
+  if (timestep_number % 25 != 0)
     return;
   
   std::vector<std::string> solution_names;
@@ -839,21 +836,6 @@ void BoussinesqFlowProblem<dim>::output_results ()  const
 
 
 
-
-template <int dim>
-void
-BoussinesqFlowProblem<dim>::project_back_temperature ()
-{
-  for (unsigned int i=0; i<solution.block(2).size(); ++i)
-    if (solution.block(2)(i) < 0)
-      solution.block(2)(i) = 0;
-    else
-      if (solution.block(2)(i) > 1)
-        solution.block(2)(i) = 1;
-}
-
-
-
 template <int dim>
 double
 BoussinesqFlowProblem<dim>::get_maximal_velocity () const

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.