]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Remove adding hanging node constrained entries to the matrix -- it doesn't make a...
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 21 May 2008 19:19:46 +0000 (19:19 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 21 May 2008 19:19:46 +0000 (19:19 +0000)
git-svn-id: https://svn.dealii.org/trunk@16165 0785d39b-7218-0410-832d-ea1e28bc413d

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

index 1dfc8440c1d43e4da2bfe40ff4dcfd99efac79b2..1cfb7b7dd25752046d07bcd58fb58882de6162b8 100644 (file)
@@ -1760,13 +1760,6 @@ void ConservationLaw<dim>::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<int> row_lengths (dof_handler.n_dofs());
@@ -3159,7 +3152,7 @@ void ConservationLaw<dim>::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;
 

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.