]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Switch back to a more interesting case.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 19 Mar 2013 03:50:23 +0000 (03:50 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 19 Mar 2013 03:50:23 +0000 (03:50 +0000)
git-svn-id: https://svn.dealii.org/trunk@28940 0785d39b-7218-0410-832d-ea1e28bc413d

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

index 772635c629db504eba88ac631f03cdebb91a9ff6..91de3ff1a38d35b978e7ce933d699a1f46ce25d6 100644 (file)
@@ -98,9 +98,6 @@ namespace Step26
     Assert (component == 0, ExcInternalError());
     Assert (dim == 2, ExcNotImplemented());
 
-    return std::sin(p[0]*numbers::PI) * std::sin(p[1]*numbers::PI);
-
-
     const double time = this->get_time();
     const double point_within_period = (time/period - std::floor(time/period));
 
@@ -154,8 +151,7 @@ namespace Step26
   template<int dim>
   void HeatEquation<dim>::setup_system()
   {
-    //    GridGenerator::hyper_L (triangulation);
-    GridGenerator::hyper_cube (triangulation, 0, 1);
+    GridGenerator::hyper_L (triangulation);
     triangulation.refine_global (5);
 
     std::cout << "Number of active cells: " << triangulation.n_active_cells()
@@ -267,7 +263,7 @@ namespace Step26
                                                    boundary_values_u_function, boundary_values);
 
           matrix_u.copy_from(mass_matrix);
-          matrix_u.add(theta * time_step, laplace_matrix); // I omit here a time_step*theta
+          matrix_u.add(theta * time_step, laplace_matrix);
           MatrixTools::apply_boundary_values(boundary_values, matrix_u,
                                              solution_u, system_rhs);
         }

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.