From: bangerth Date: Wed, 17 Oct 2007 16:38:45 +0000 (+0000) Subject: Revert to using DG elements for the temperature field. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4b1f45f2235b83e28276459e1c6f394a3a127b3e;p=dealii-svn.git Revert to using DG elements for the temperature field. git-svn-id: https://svn.dealii.org/trunk@15339 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-22/step-22.cc b/deal.II/examples/step-22/step-22.cc index 13da6db65a..43690c31e7 100644 --- a/deal.II/examples/step-22/step-22.cc +++ b/deal.II/examples/step-22/step-22.cc @@ -38,6 +38,7 @@ #include #include +#include #include #include @@ -364,9 +365,9 @@ BoussinesqFlowProblem::BoussinesqFlowProblem (const unsigned int degree) degree (degree), fe (FE_Q(degree+1), dim, FE_Q(degree), 1, - FE_Q(degree), 1), + FE_DGQ(degree-1), 1), dof_handler (triangulation), - n_refinement_steps (6), + n_refinement_steps (5), time_step (0), rebuild_preconditioner (true) {} @@ -628,6 +629,7 @@ void BoussinesqFlowProblem::assemble_rhs_T () Vector local_rhs (dofs_per_cell); std::vector > old_solution_values(n_q_points, Vector(dim+2)); + std::vector > old_solution_values_face(n_face_q_points, Vector(dim+2)); std::vector > old_solution_values_face_neighbor(n_face_q_points, Vector(dim+2)); std::vector > present_solution_values(n_q_points, Vector(dim+2)); @@ -839,7 +841,7 @@ void BoussinesqFlowProblem::output_results () const DataOut::type_dof_data, data_component_interpretation); - data_out.build_patches (degree); + data_out.build_patches (degree+1); std::ostringstream filename; filename << "solution-" << Utilities::int_to_string(timestep_number, 4) << ".vtk";