From 918b40de88b4536bd51efd8f3e685b8a5348d68b Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 14 Aug 2008 18:59:51 +0000 Subject: [PATCH] Default to cube geometry in 3d. Use proper direction of gravity. git-svn-id: https://svn.dealii.org/trunk@16554 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-31/step-31.cc | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index ff1431057b..9e7a0e84ee 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -1542,7 +1542,7 @@ void BoussinesqFlowProblem::assemble_system () * fe_values.JxW(q); const Point gravity = ( (dim == 2) ? (Point (0,1)) : - (Point (0,1,0)) ); + (Point (0,0,1)) ); local_rhs(i) += (Rayleigh_number * gravity * phi_u[i] * old_temperature)* @@ -2334,13 +2334,15 @@ void BoussinesqFlowProblem::run () case 3: { - GridGenerator::hyper_shell (triangulation, - Point(), 0.5, 1.0); +// GridGenerator::hyper_shell (triangulation, +// Point(), 0.5, 1.0); - static HyperShellBoundary boundary; - triangulation.set_boundary (0, boundary); +// static HyperShellBoundary boundary; +// triangulation.set_boundary (0, boundary); + + GridGenerator::hyper_cube (triangulation); - triangulation.refine_global (2); + triangulation.refine_global (3); break; } -- 2.39.5