]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Default to cube geometry in 3d. Use proper direction of gravity.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 14 Aug 2008 18:59:51 +0000 (18:59 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 14 Aug 2008 18:59:51 +0000 (18:59 +0000)
git-svn-id: https://svn.dealii.org/trunk@16554 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-31/step-31.cc

index ff1431057bf9b11ca16f579b621e15ea1807ae1b..9e7a0e84eecf832b7e518d44ba2c7e62fe49b7f3 100644 (file)
@@ -1542,7 +1542,7 @@ void BoussinesqFlowProblem<dim>::assemble_system ()
                                       * fe_values.JxW(q);
 
              const Point<dim> gravity = ( (dim == 2) ? (Point<dim> (0,1)) : 
-                                                       (Point<dim> (0,1,0)) );
+                                                       (Point<dim> (0,0,1)) );
 
              local_rhs(i) += (Rayleigh_number *
                               gravity * phi_u[i] * old_temperature)*
@@ -2334,13 +2334,15 @@ void BoussinesqFlowProblem<dim>::run ()
 
       case 3:
       {
-       GridGenerator::hyper_shell (triangulation,
-                                   Point<dim>(), 0.5, 1.0);
+//     GridGenerator::hyper_shell (triangulation,
+//                                 Point<dim>(), 0.5, 1.0);
 
-       static HyperShellBoundary<dim> boundary;
-       triangulation.set_boundary (0, boundary);
+//     static HyperShellBoundary<dim> boundary;
+//     triangulation.set_boundary (0, boundary);
+
+       GridGenerator::hyper_cube (triangulation);
 
-       triangulation.refine_global (2);
+       triangulation.refine_global (3);
 
        break;
       }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.