diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh index baa2c52e7a13916e4db7a5ec42782901ff13ba72..61cffac213b1dfff6dab740eb7e37d952ce27e0b 100644 --- a/dune/matrix-vector/triangularsolve.hh +++ b/dune/matrix-vector/triangularsolve.hh @@ -3,6 +3,11 @@ namespace Dune { namespace MatrixVector { + + /** + * \brief Solves L*x=b where L is a lower triangular matrix. + * \param transpose Set true if L is passed as an upper triangular matrix. + */ template <class Matrix, class Vector, class BitVector> static void lowerTriangularSolve(Matrix const& L, Vector b, Vector& x, BitVector const* ignore, @@ -40,6 +45,10 @@ namespace MatrixVector { } } + /** + * \brief Solves U*x=b where U is an upper triangular matrix. + * \param transpose Set true if U is passed as a lower triangular matrix. + */ template <class Matrix, class Vector, class BitVector> static void upperTriangularSolve(Matrix const& U, Vector b, Vector& x, BitVector const* ignore,