Skip to content
Snippets Groups Projects
Commit dad27284 authored by Max Kahnt's avatar Max Kahnt
Browse files

Add docu.

parent f4f36927
Branches
Tags
No related merge requests found
...@@ -3,6 +3,11 @@ ...@@ -3,6 +3,11 @@
namespace Dune { namespace Dune {
namespace MatrixVector { 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> template <class Matrix, class Vector, class BitVector>
static void lowerTriangularSolve(Matrix const& L, Vector b, Vector& x, static void lowerTriangularSolve(Matrix const& L, Vector b, Vector& x,
BitVector const* ignore, BitVector const* ignore,
...@@ -40,6 +45,10 @@ namespace MatrixVector { ...@@ -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> template <class Matrix, class Vector, class BitVector>
static void upperTriangularSolve(Matrix const& U, Vector b, Vector& x, static void upperTriangularSolve(Matrix const& U, Vector b, Vector& x,
BitVector const* ignore, BitVector const* ignore,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment