diff --git a/dune/istl/multitypeblockmatrix.hh b/dune/istl/multitypeblockmatrix.hh index cf3f5e5ccc1ffed678739098add02c5283062d38..c2d3e79a231f905a1aa04f079c0f86d4b01751f2 100644 --- a/dune/istl/multitypeblockmatrix.hh +++ b/dune/istl/multitypeblockmatrix.hh @@ -164,7 +164,7 @@ namespace Dune { public: /** - * do y -= A x in next row + * do y += A x in next row */ static void umv(TVecY& y, const TMatrix& A, const TVecX& x) { static const int rowlen = mpl::size< typename mpl::at_c<TMatrix,crow>::type >::value; diff --git a/dune/istl/tutorial/example.cc b/dune/istl/tutorial/example.cc index cc52ba88b06ae948fdfe7c401c48dffa32bfbe8e..53d7722a05a273550105e5a82eb75e2df229da68 100644 --- a/dune/istl/tutorial/example.cc +++ b/dune/istl/tutorial/example.cc @@ -584,7 +584,7 @@ void test_MultiTypeBlockVector_MultiTypeBlockMatrix() { Dune::MatrixAdapter<CM_BCRS,TestVector,TestVector> op(A); // make linear operator from A Dune::SeqJac<CM_BCRS,TestVector,TestVector,2> jac(A,1,1); // Jacobi preconditioner Dune::SeqGS<CM_BCRS,TestVector,TestVector,2> gs(A,1,1); // GS preconditioner - Dune::SeqSOR<CM_BCRS,TestVector,TestVector,2> sor(A,1,1.9520932); // SSOR preconditioner + Dune::SeqSOR<CM_BCRS,TestVector,TestVector,2> sor(A,1,1.9520932); // SOR preconditioner Dune::SeqSSOR<CM_BCRS,TestVector,TestVector,2> ssor(A,1,1.0); // SSOR preconditioner Dune::LoopSolver<TestVector> loop(op,gs,1E-4,18000,2); // an inverse operator