From: bangerth Date: Tue, 19 Mar 2013 03:50:23 +0000 (+0000) Subject: Switch back to a more interesting case. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=91f4f3ef1bfc8c00c84ecc669d9886e8c6726647;p=dealii-svn.git Switch back to a more interesting case. git-svn-id: https://svn.dealii.org/trunk@28940 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-26/step-26.cc b/deal.II/examples/step-26/step-26.cc index 772635c629..91de3ff1a3 100644 --- a/deal.II/examples/step-26/step-26.cc +++ b/deal.II/examples/step-26/step-26.cc @@ -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 void HeatEquation::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); }