From 533bb936cdb28ae5c75ae9a9c915e33cce562e31 Mon Sep 17 00:00:00 2001 From: bangerth Date: Wed, 21 May 2008 19:19:46 +0000 Subject: [PATCH] Remove adding hanging node constrained entries to the matrix -- it doesn't make a difference in accuracy or runtime. Also compute the predictor using a full time step instead of a half step. git-svn-id: https://svn.dealii.org/trunk@16165 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-33/step-33.cc | 9 +-------- 1 file changed, 1 insertion(+), 8 deletions(-) diff --git a/deal.II/examples/step-33/step-33.cc b/deal.II/examples/step-33/step-33.cc index 1dfc8440c1..1cfb7b7dd2 100644 --- a/deal.II/examples/step-33/step-33.cc +++ b/deal.II/examples/step-33/step-33.cc @@ -1760,13 +1760,6 @@ void ConservationLaw::setup_system () CompressedSparsityPattern sparsity_pattern (dof_handler.n_dofs(), dof_handler.n_dofs()); DoFTools::make_sparsity_pattern (dof_handler, sparsity_pattern); - - ConstraintMatrix hanging_node_constraints; - DoFTools::make_hanging_node_constraints (dof_handler, - hanging_node_constraints); - hanging_node_constraints.close (); - - hanging_node_constraints.condense (sparsity_pattern); sparsity_pattern.compress(); std::vector row_lengths (dof_handler.n_dofs()); @@ -3159,7 +3152,7 @@ void ConservationLaw::run () } predictor = current_solution; - predictor.sadd(3/2.0, -1/2.0, old_solution); + predictor.sadd (2.0, -1.0, old_solution); old_solution = current_solution; -- 2.39.5