From: guido Date: Mon, 15 Mar 2004 23:36:34 +0000 (+0000) Subject: TODO's removed X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e84c1fc92dbc038ca633cea040fa79a737c5a271;p=dealii-svn.git TODO's removed git-svn-id: https://svn.dealii.org/trunk@8759 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 d49cb052a5..7834575749 100644 --- a/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc +++ b/deal.II/deal.II/source/numerics/matrices.all_dimensions.cc @@ -35,6 +35,8 @@ MatrixTools::apply_boundary_values (const std::map &boundar Vector &right_hand_side, const bool preserve_symmetry) { + // Require that diagonals are first + // in each row Assert (matrix.get_sparsity_pattern().optimize_diagonal(), typename SparsityPattern::ExcDiagonalNotOptimized()); Assert (matrix.n() == right_hand_side.size(), @@ -84,7 +86,6 @@ MatrixTools::apply_boundary_values (const std::map &boundar // we shall not set // matrix.global_entry( // sparsity_rowstart[dof.first]) -//TODO[GK]: This is going to break when the diagonal element of the matrix is not stored const unsigned int last = sparsity_rowstart[dof_number+1]; for (unsigned int j=sparsity_rowstart[dof_number]+1; j &boundar // matrix so the diagonal // element is the first of // this row. -//TODO[GK]: This is going to break when the diagonal element of the matrix is not stored const unsigned int last = local_sparsity.get_rowstart_indices()[block_index.second+1], first = (block_col == block_index.first ?