From: frohne Date: Tue, 30 Oct 2012 16:10:07 +0000 (+0000) Subject: set all contraints to zero for solve X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=37edb30ed14a431ceecdcead93715cc960117df8;p=dealii-svn.git set all contraints to zero for solve git-svn-id: https://svn.dealii.org/trunk@27246 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc index 1370581865..919838848b 100644 --- a/deal.II/examples/step-42/step-42.cc +++ b/deal.II/examples/step-42/step-42.cc @@ -13,7 +13,7 @@ // @sect3{Include files} // We are using the the same - // include files as in step-42: + // include files as in step-41: #include #include @@ -752,7 +752,7 @@ namespace Step42 cell->get_dof_indices (local_dof_indices); constraints.distribute_local_to_global (cell_matrix, cell_rhs, local_dof_indices, - system_matrix_newton, system_rhs_newton, true); + system_matrix_newton, system_rhs_newton, false); }; system_matrix_newton.compress (); @@ -1063,7 +1063,8 @@ namespace Step42 TrilinosWrappers::MPI::Vector distributed_solution (system_rhs_newton); distributed_solution = solution; - constraints_hanging_nodes.set_zero (distributed_solution); + // constraints_hanging_nodes.set_zero (distributed_solution); + constraints.set_zero (distributed_solution); // Solving iterative @@ -1080,10 +1081,10 @@ namespace Step42 MPI_Barrier (mpi_communicator); t.restart(); - // ReductionControl reduction_control (10000, 1e-15, 1e-4); - // SolverCG - // solver (reduction_control, mpi_communicator); - // solver.solve (system_matrix_newton, distributed_solution, system_rhs_newton, preconditioner_u); + // ReductionControl solver_control (10000, 1e-15, 1e-4); + // SolverCG + // solver (solver_control, mpi_communicator); + // solver.solve (system_matrix_newton, distributed_solution, system_rhs_newton, preconditioner_u); PrimitiveVectorMemory mem; TrilinosWrappers::MPI::Vector tmp (system_rhs_newton);