diff options
Diffstat (limited to 'src/routines/level2/xtrsv.hpp')
-rw-r--r-- | src/routines/level2/xtrsv.hpp | 17 |
1 files changed, 15 insertions, 2 deletions
diff --git a/src/routines/level2/xtrsv.hpp b/src/routines/level2/xtrsv.hpp index 4a73b5eb..dc3f32f0 100644 --- a/src/routines/level2/xtrsv.hpp +++ b/src/routines/level2/xtrsv.hpp @@ -7,7 +7,9 @@ // Author(s): // Cedric Nugteren <www.cedricnugteren.nl> // -// This file implements the Xtrsv routine. +// This file implements the Xtrsv routine. It uses a block-algorithm and performs small triangular +// forward and backward substitutions on the diagonal parts of the matrix in combination with larger +// GEMV computation on the remainder of the matrix. // // ================================================================================================= @@ -25,9 +27,12 @@ class Xtrsv: public Xgemv<T> { public: // Uses the generic matrix-vector routine + using Xgemv<T>::routine_name_; using Xgemv<T>::queue_; using Xgemv<T>::context_; - using Xgemv<T>::MatVec; + using Xgemv<T>::device_; + using Xgemv<T>::db_; + using Xgemv<T>::DoGemv; // Constructor Xtrsv(Queue &queue, EventPointer event, const std::string &name = "TRSV"); @@ -38,6 +43,14 @@ class Xtrsv: public Xgemv<T> { const size_t n, const Buffer<T> &a_buffer, const size_t a_offset, const size_t a_ld, const Buffer<T> &x_buffer, const size_t x_offset, const size_t x_inc); + + // Performs forward or backward substitution on a small triangular matrix + void Substitution(const Layout layout, const Triangle triangle, + const Transpose a_transpose, const Diagonal diagonal, + const size_t n, + const Buffer<T> &a_buffer, const size_t a_offset, const size_t a_ld, + const Buffer<T> &b_buffer, const size_t b_offset, const size_t b_inc, + const Buffer<T> &x_buffer, const size_t offset_x, const size_t x_inc); }; // ================================================================================================= |