From: guido Date: Wed, 31 May 2000 04:42:59 +0000 (+0000) Subject: bug removed X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8a8cedb2fcdbf6a2c13390e4d6abc56f221e27e8;p=dealii-svn.git bug removed git-svn-id: https://svn.dealii.org/trunk@2981 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/sparse_matrix.h b/deal.II/lac/include/lac/sparse_matrix.h index 88faaad694..035d296032 100644 --- a/deal.II/lac/include/lac/sparse_matrix.h +++ b/deal.II/lac/include/lac/sparse_matrix.h @@ -822,7 +822,8 @@ void SparseMatrix::set (const unsigned int i, template inline void SparseMatrix::add (const unsigned int i, const unsigned int j, - const number value) { + const number value) +{ Assert (cols != 0, ExcMatrixNotInitialized()); const unsigned int index = cols->operator()(i,j); @@ -842,9 +843,11 @@ number SparseMatrix::operator () (const unsigned int i, const unsigned int j) const { Assert (cols != 0, ExcMatrixNotInitialized()); - Assert (cols->operator()(i,j) != SparsityPattern::invalid_entry, - ExcInvalidIndex(i,j)); - return val[cols->operator()(i,j)]; + + const unsigned int index = cols->operator()(i,j); + if (index == SparsityPattern::invalid_entry) + return 0.; + return val[index]; };