From cb5015cd4a0292790271b0b57e26a24b9b95be23 Mon Sep 17 00:00:00 2001 From: guido Date: Wed, 28 May 2003 14:58:27 +0000 Subject: [PATCH] PSOR checked git-svn-id: https://svn.dealii.org/trunk@7677 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/sparse_matrix.templates.h | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) 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]]; -- 2.39.5