From: Wolfgang Bangerth Date: Fri, 18 Feb 2011 17:02:50 +0000 (+0000) Subject: Need to also distribute system_rhs. X-Git-Tag: v8.0.0~4308 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a6316240cce103e114fe88b985b2cd12e3e968ed;p=dealii.git Need to also distribute system_rhs. git-svn-id: https://svn.dealii.org/trunk@23394 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-46/step-46.cc b/deal.II/examples/step-46/step-46.cc index 1434dd2d78..64aef01aa9 100644 --- a/deal.II/examples/step-46/step-46.cc +++ b/deal.II/examples/step-46/step-46.cc @@ -299,9 +299,13 @@ void StokesProblem::assemble_system () local_matrix(i,j) = local_matrix(j,i); cell->get_dof_indices (local_dof_indices); - constraints.distribute_local_to_global (local_matrix, + + // local_rhs==0, but need to do + // this here because of + // boundary values + constraints.distribute_local_to_global (local_matrix, local_rhs, local_dof_indices, - system_matrix); + system_matrix, system_rhs); } } @@ -390,7 +394,7 @@ void StokesProblem::run () triangulation.refine_global (3-dim); - for (unsigned int refinement_cycle = 0; refinement_cycle<1; + for (unsigned int refinement_cycle = 0; refinement_cycle<6; ++refinement_cycle) { std::cout << "Refinement cycle " << refinement_cycle << std::endl;