Another baby step towards getting rid of global_entry().

git-svn-id: https://svn.dealii.org/trunk@27973 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/source/numerics/matrix_tools.cc b/deal.II/source/numerics/matrix_tools.cc
index 9e33193338ef6f1151348632ece10c1ff11ba648..61717252ba7320a88412b834b526e73576c48947 100644 (file)
--- a/deal.II/source/numerics/matrix_tools.cc
+++ b/deal.II/source/numerics/matrix_tools.cc
@@ -2299,13 +2299,12 @@ namespace MatrixTools
                 // don't use the
                 // diagonal element of
                 // the diagonal block
-                const unsigned int
-                first = (block_index.first == block_row ?
-                         transpose_sparsity.get_rowstart_indices()[block_index.second]+1 :
-                         transpose_sparsity.get_rowstart_indices()[block_index.second]),
-                        last  = transpose_sparsity.get_rowstart_indices()[block_index.second+1];
-
-                for (unsigned int j=first; j<last; ++j)
+                for (typename SparseMatrix<number>::iterator
+                     q = (block_index.first == block_row ?
+                         transpose_matrix.begin(block_index.second)+1 :
+                         transpose_matrix.begin(block_index.second));
+                    q != transpose_matrix.end(block_index.second);
+                    ++q)
                   {
                     // get the number
                     // of the column in
@@ -2317,7 +2316,7 @@ namespace MatrixTools
                     // block which has
                     // an entry in the
                     // interesting row
-                    const unsigned int row = transpose_sparsity.get_column_numbers()[j];
+                    const unsigned int row = q->column();
 
                     // find the
                     // position of