diff --git a/dune/solvers/iterationsteps/linegsstep.cc b/dune/solvers/iterationsteps/linegsstep.cc index 26dcc79d3f28e89eecfff660a9190f6518c4fe3d..a8640dbf5d6f229a6cf2df53e1074ff9842f25c1 100755 --- a/dune/solvers/iterationsteps/linegsstep.cc +++ b/dune/solvers/iterationsteps/linegsstep.cc @@ -38,7 +38,7 @@ void LineGSStep<OperatorType, DiscFuncType, BitVectorType >::iterate() if ( (*this->ignoreNodes_)[blockStructure_[b_num][k]][0]) permuted_r_i[k] = 0.0; else - residual( blockStructure_[b_num][k], permuted_r_i[k]); // get r[p(i+k)] + this->residual( blockStructure_[b_num][k], permuted_r_i[k]); // get r[p(i+k)] } // permuted_r_i[k] is the residual for r[p(i+k)], which means the residual \tilde{r}[i+k] // note: calling the method 'residual' with the permuted index implies that there is no additional permutation required within 'residual' diff --git a/dune/solvers/iterationsteps/projectedlinegsstep.cc b/dune/solvers/iterationsteps/projectedlinegsstep.cc index 7560df3801cdaf41abdfb88caf6ee2017f86c487..5784106df1db048f7b1eb5147cabe8f9d20f6139 100755 --- a/dune/solvers/iterationsteps/projectedlinegsstep.cc +++ b/dune/solvers/iterationsteps/projectedlinegsstep.cc @@ -191,7 +191,7 @@ void ProjectedLineGSStep<MatrixType, VectorType, BitVectorType>::iterate() if ( (*this->ignoreNodes_)[this->blockStructure_[b_num][k]][0]) permuted_r_i[k] = 0.0; else - residual( this->blockStructure_[b_num][k], permuted_r_i[k]); // get r[p(i+k)] + this->residual( this->blockStructure_[b_num][k], permuted_r_i[k]); // get r[p(i+k)] } // permuted_r_i[k] is the residual for r[p(i+k)], which means the residual \tilde{r}[i+k] // note: calling the method 'residual' with the permuted index implies that there is no additional permutation required within 'residual'