From: guido Date: Wed, 28 May 2003 14:58:27 +0000 (+0000) Subject: PSOR checked X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cb5015cd4a0292790271b0b57e26a24b9b95be23;p=dealii-svn.git PSOR checked git-svn-id: https://svn.dealii.org/trunk@7677 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/sparse_matrix.templates.h b/deal.II/lac/include/lac/sparse_matrix.templates.h index 3a9659c2db..67c396cb43 100644 --- a/deal.II/lac/include/lac/sparse_matrix.templates.h +++ b/deal.II/lac/include/lac/sparse_matrix.templates.h @@ -1018,11 +1018,17 @@ SparseMatrix::PSOR ( { const unsigned int row = permutation[urow]; somenumber s = dst(row); + for (unsigned int j=cols->rowstart[row]; jrowstart[row+1]; ++j) { const unsigned int col = cols->colnums[j]; if (inverse_permutation[col] < urow) - s -= val[j] * dst(col); + { + s -= val[j] * dst(col); + } + + std::cerr << std::endl; + } dst(row) = s * om / val[cols->rowstart[row]];