diff options
Diffstat (limited to 'src/clblast.cc')
| -rw-r--r-- | src/clblast.cc | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/src/clblast.cc b/src/clblast.cc index eddb8022..c04462f3 100644 --- a/src/clblast.cc +++ b/src/clblast.cc @@ -22,6 +22,7 @@ // BLAS level-2 includes #include "internal/routines/level2/xgemv.h" +#include "internal/routines/level2/xsymv.h" // BLAS level-3 includes #include "internal/routines/level3/xgemm.h" @@ -36,6 +37,7 @@ namespace clblast { // ================================================================================================= // BLAS level-1 (vector-vector) routines +// ================================================================================================= // AXPY template <typename T> @@ -75,6 +77,7 @@ template StatusCode Axpy<double2>(const size_t, const double2, // ================================================================================================= // BLAS level-2 (matrix-vector) routines +// ================================================================================================= // GEMV template <typename T> @@ -125,7 +128,46 @@ template StatusCode Gemv<double2>(const Layout, const Transpose, cl_command_queue*, cl_event*); // ================================================================================================= + +// SYMV +template <typename T> +StatusCode Symv(const Layout layout, const Triangle triangle, + const size_t n, const T alpha, + const cl_mem a_buffer, const size_t a_offset, const size_t a_ld, + const cl_mem x_buffer, const size_t x_offset, const size_t x_inc, const T beta, + cl_mem y_buffer, const size_t y_offset, const size_t y_inc, + cl_command_queue* queue, cl_event* event) { + + auto queue_cpp = Queue(*queue); + auto event_cpp = Event(*event); + auto routine = Xsymv<T>(queue_cpp, event_cpp); + + // Compiles the routine's device kernels + auto status = routine.SetUp(); + if (status != StatusCode::kSuccess) { return status; } + + // Runs the routine + return routine.DoSymv(layout, triangle, n, alpha, + Buffer<T>(a_buffer), a_offset, a_ld, + Buffer<T>(x_buffer), x_offset, x_inc, beta, + Buffer<T>(y_buffer), y_offset, y_inc); +} +template StatusCode Symv<float>(const Layout, const Triangle, + const size_t, const float, + const cl_mem, const size_t, const size_t, + const cl_mem, const size_t, const size_t, const float, + cl_mem, const size_t, const size_t, + cl_command_queue*, cl_event*); +template StatusCode Symv<double>(const Layout, const Triangle, + const size_t, const double, + const cl_mem, const size_t, const size_t, + const cl_mem, const size_t, const size_t, const double, + cl_mem, const size_t, const size_t, + cl_command_queue*, cl_event*); + +// ================================================================================================= // BLAS level-3 (matrix-matrix) routines +// ================================================================================================= // GEMM template <typename T> |