From: kronbichler Date: Fri, 17 Oct 2008 07:55:53 +0000 (+0000) Subject: Fixed an issue with boundary conditions. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0a4fb9cad265080f71a334115b2b43f847b3faa4;p=dealii-svn.git Fixed an issue with boundary conditions. git-svn-id: https://svn.dealii.org/trunk@17256 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc b/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc index 869ad7104f..4366991f8e 100644 --- a/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc +++ b/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc @@ -745,28 +745,33 @@ namespace TrilinosWrappers (dof->first < local_range.second)) constrained_rows.push_back (dof->first); - // then eliminate these rows and set - // their diagonal entry to what we have - // determined above. if the value already - // is nonzero, it will be preserved, - // in accordance with the basic + // then eliminate these rows and + // set their diagonal entry to + // what we have determined + // above. if the value already is + // nonzero, it will be preserved, + // in accordance with the basic // matrix classes in deal.II. matrix.clear_rows (constrained_rows, average_nonzero_diagonal_entry); - // the next thing is to set right hand - // side to the wanted value. there's one - // drawback: if we write to individual - // vector elements, then we have to do - // that on all processors. however, some + // the next thing is to set right + // hand side to the wanted + // value. there's one drawback: + // if we write to individual + // vector elements, then we have + // to do that on all + // processors. however, some // processors may not need to set - // anything because their chunk of - // matrix/rhs do not contain any boundary - // nodes. therefore, rather than using - // individual calls, we use one call for - // all elements, thereby making sure that - // all processors call this function, - // even if some only have an empty set of - // elements to set + // anything because their chunk + // of matrix/rhs do not contain + // any boundary nodes. therefore, + // rather than using individual + // calls, we use one call for all + // elements, thereby making sure + // that all processors call this + // function, even if some only + // have an empty set of elements + // to set right_hand_side.compress (); solution.compress (); @@ -784,10 +789,10 @@ namespace TrilinosWrappers } solution.set (indices, solution_values); - // now also set appropriate values for - // the rhs + // now also set appropriate + // values for the rhs for (unsigned int i=0; i