From dee16a5d6a983bef0e9b282f8ef1535ab100ae63 Mon Sep 17 00:00:00 2001 From: bangerth Date: Wed, 10 Oct 2007 02:52:25 +0000 Subject: [PATCH] Remove pressure rhs, we won't need it here... git-svn-id: https://svn.dealii.org/trunk@15278 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-22/step-22.cc | 33 ++++------------------------- 1 file changed, 4 insertions(+), 29 deletions(-) diff --git a/deal.II/examples/step-22/step-22.cc b/deal.II/examples/step-22/step-22.cc index 21d24f6934..6edf2cab6e 100644 --- a/deal.II/examples/step-22/step-22.cc +++ b/deal.II/examples/step-22/step-22.cc @@ -91,28 +91,6 @@ class BoussinesqFlowProblem -template -class PressureRightHandSide : public Function -{ - public: - PressureRightHandSide () : Function(1) {} - - virtual double value (const Point &p, - const unsigned int component = 0) const; -}; - - - -template -double -PressureRightHandSide::value (const Point &/*p*/, - const unsigned int /*component*/) const -{ - return 0; -} - - - template class Buoyancy : public Function { @@ -130,8 +108,10 @@ void Buoyancy::vector_value (const Point &p, Vector &values) const { + Assert (values.size() == dim, ExcInternalError()); + values = 0; - values(dim-1) = p[dim-1]; + values(dim-1) = (p.distance (Point(.5,.5)) < .2 ? 1 : 0); } @@ -515,11 +495,9 @@ void BoussinesqFlowProblem::assemble_system () std::vector local_dof_indices (dofs_per_cell); - const PressureRightHandSide pressure_right_hand_side; const Buoyancy buoyancy; const PressureBoundaryValues pressure_boundary_values; - std::vector pressure_rhs_values (n_q_points); std::vector > buoyancy_values (n_q_points, Vector(dim)); std::vector boundary_values (n_face_q_points); @@ -539,8 +517,6 @@ void BoussinesqFlowProblem::assemble_system () fe_values.get_function_values (old_solution, old_solution_values); - pressure_right_hand_side.value_list (fe_values.get_quadrature_points(), - pressure_rhs_values); buoyancy.vector_value_list (fe_values.get_quadrature_points(), buoyancy_values); for (unsigned int q=0; q::assemble_system () } Assert (dim == 2, ExcInternalError()); - local_rhs(i) += (-phi_i_p * pressure_rhs_values[q] - +phi_i_u[0] * buoyancy_values[q](0) + local_rhs(i) += ( phi_i_u[0] * buoyancy_values[q](0) +phi_i_u[1] * buoyancy_values[q](1))* fe_values.JxW(q); } -- 2.39.5