From b041faeb187de8cb11e6416269600095620de986 Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 14 Aug 2008 23:04:21 +0000 Subject: [PATCH] Remove the function that outputs the artificial viscosity parameter. git-svn-id: https://svn.dealii.org/trunk@16558 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-31/step-31.cc | 84 ----------------------------- 1 file changed, 84 deletions(-) diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index 25c9c23216..5ed0f8a9b1 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -325,7 +325,6 @@ class BoussinesqFlowProblem double get_maximal_velocity () const; double get_maximal_temperature () const; void solve (); - void test (Vector &) const; void output_results () const; void refine_mesh (); @@ -2004,83 +2003,6 @@ void BoussinesqFlowProblem::assemble_rhs_T () -template -void -BoussinesqFlowProblem::test (Vector &viscosity) const -{ - QGauss quadrature_formula(degree+2); - FEValues fe_values (fe, quadrature_formula, - update_values | update_gradients | - update_hessians | - update_quadrature_points | update_JxW_values); - - const unsigned int n_q_points = quadrature_formula.size(); - - std::vector > present_solution_values (n_q_points, - Vector(dim+2)); - std::vector > old_solution_values(n_q_points, - Vector(dim+2)); - std::vector > old_old_solution_values(n_q_points, - Vector(dim+2)); - std::vector > > - old_solution_grads(n_q_points, - std::vector >(dim+2)); - std::vector > > - old_old_solution_grads(n_q_points, - std::vector >(dim+2)); - std::vector > > old_solution_hessians( - n_q_points, - std::vector >(dim+2)); - std::vector > > old_old_solution_hessians( - n_q_points, - std::vector >(dim+2)); - RightHandSide right_hand_side; - std::vector gamma_values (n_q_points); - - const FEValuesExtractors::Scalar temperature (dim+1); - - const double global_u_infty = get_maximal_velocity(); - const double global_T_infty = get_maximal_temperature(); - const double global_Omega_diameter = GridTools::diameter (triangulation); - - typename DoFHandler::active_cell_iterator - cell = dof_handler.begin_active(), - endc = dof_handler.end(); - - for (unsigned int index=0; cell!=endc; ++cell, ++index) - { - fe_values.reinit (cell); - - fe_values.get_function_values (solution, present_solution_values); - fe_values.get_function_values (old_solution, old_solution_values); - fe_values.get_function_values (old_old_solution, old_old_solution_values); - - fe_values.get_function_gradients (old_solution, old_solution_grads); - fe_values.get_function_gradients (old_old_solution, old_old_solution_grads); - - fe_values.get_function_hessians (old_solution, old_solution_hessians); - fe_values.get_function_hessians (old_old_solution, old_old_solution_hessians); - - const double kappa = 1e-6; - right_hand_side.value_list (fe_values.get_quadrature_points(), - gamma_values, dim+1); - - viscosity(index) = - (timestep_number == 0 - ? - 5e-3 * cell->diameter() - : - compute_viscosity (present_solution_values, old_solution_values, - old_old_solution_values, old_solution_grads, old_old_solution_grads, - old_solution_hessians, old_old_solution_hessians, gamma_values, - kappa, global_u_infty, global_T_infty, - global_Omega_diameter, cell->diameter(), - old_time_step)); - } -} - - - // @sect4{BoussinesqFlowProblem::solve} template @@ -2199,9 +2121,6 @@ void BoussinesqFlowProblem::output_results () const solution_names.push_back ("p"); solution_names.push_back ("T"); - Vector viscosity (triangulation.n_active_cells()); - test (viscosity); - DataOut data_out; data_out.attach_dof_handler (dof_handler); @@ -2216,9 +2135,6 @@ void BoussinesqFlowProblem::output_results () const data_out.add_data_vector (solution, solution_names, DataOut::type_dof_data, data_component_interpretation); - - data_out.add_data_vector (viscosity, "viscosity"); - data_out.build_patches (); std::ostringstream filename; -- 2.39.5