]> https://gitweb.dealii.org/ - dealii.git/commitdiff
bug removed
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Wed, 31 May 2000 04:42:59 +0000 (04:42 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Wed, 31 May 2000 04:42:59 +0000 (04:42 +0000)
git-svn-id: https://svn.dealii.org/trunk@2981 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/sparse_matrix.h

index 88faaad694ddb45876921251e5430af8580f33b1..035d296032990a9de66b8ef232970a94632eba06 100644 (file)
@@ -822,7 +822,8 @@ void SparseMatrix<number>::set (const unsigned int i,
 template <typename number>
 inline
 void SparseMatrix<number>::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<number>::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];
 };
 
 

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.