]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
TODO's removed
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 15 Mar 2004 23:36:34 +0000 (23:36 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 15 Mar 2004 23:36:34 +0000 (23:36 +0000)
git-svn-id: https://svn.dealii.org/trunk@8759 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/numerics/matrices.all_dimensions.cc

index d49cb052a531a08ed91352c7e735252c59307630..7834575749d7b4075cc1f600867fc5a98eacad2a 100644 (file)
@@ -35,6 +35,8 @@ MatrixTools::apply_boundary_values (const std::map<unsigned int,double> &boundar
                                    Vector<number>   &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<unsigned int,double> &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<last; ++j)
        matrix.global_entry(j) = 0.;
@@ -311,7 +312,6 @@ MatrixTools::apply_boundary_values (const std::map<unsigned int,double> &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 ?

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.