From: bangerth Date: Sat, 20 Sep 2008 03:24:05 +0000 (+0000) Subject: Move definition of the Rayleigh number into the EquationData namespace as well. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=18b028d981c86e30d7cb17ed8a24e63d1e17d4ea;p=dealii-svn.git Move definition of the Rayleigh number into the EquationData namespace as well. git-svn-id: https://svn.dealii.org/trunk@16872 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 e550f9b7c3..46821aaa83 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -107,6 +107,8 @@ namespace EquationData // define viscosity const double eta = 1; const double kappa = 1e-6; + const double Rayleigh_number = 10; + template class PressureBoundaryValues : public Function @@ -1277,10 +1279,6 @@ void BoussinesqFlowProblem::assemble_stokes_system () // Stokes system as a source // term in the momentum equation. // - // Then, we create a variable - // to hold the Rayleigh number, - // the measure of buoyancy. - // // The set of vectors we create // next hold the evaluations of // the basis functions that will @@ -1301,8 +1299,6 @@ void BoussinesqFlowProblem::assemble_stokes_system () std::vector old_temperature_values(n_q_points); - const double Rayleigh_number = 10; - std::vector > phi_u (dofs_per_cell); std::vector > grads_phi_u (dofs_per_cell); std::vector div_phi_u (dofs_per_cell); @@ -1392,7 +1388,7 @@ void BoussinesqFlowProblem::assemble_stokes_system () const Point gravity = ( (dim == 2) ? (Point (0,1)) : (Point (0,0,1)) ); for (unsigned int i=0; i class PressureBoundaryValues : public Function @@ -859,8 +861,6 @@ void BoussinesqFlowProblem::assemble_stokes_system () std::vector old_temperature_values(n_q_points); - const double Rayleigh_number = 10; - std::vector > phi_u (dofs_per_cell); std::vector > grads_phi_u (dofs_per_cell); std::vector div_phi_u (dofs_per_cell); @@ -915,7 +915,7 @@ void BoussinesqFlowProblem::assemble_stokes_system () const Point gravity = stokes_fe_values.quadrature_point(q) / stokes_fe_values.quadrature_point(q).norm(); for (unsigned int i=0; i