summaryrefslogtreecommitdiff
path: root/src/routines/level2/xtbmv.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/routines/level2/xtbmv.cpp')
-rw-r--r--src/routines/level2/xtbmv.cpp44
1 files changed, 22 insertions, 22 deletions
diff --git a/src/routines/level2/xtbmv.cpp b/src/routines/level2/xtbmv.cpp
index e315c544..f4a58ed2 100644
--- a/src/routines/level2/xtbmv.cpp
+++ b/src/routines/level2/xtbmv.cpp
@@ -29,17 +29,15 @@ Xtbmv<T>::Xtbmv(Queue &queue, EventPointer event, const std::string &name):
// The main routine
template <typename T>
-StatusCode Xtbmv<T>::DoTbmv(const Layout layout, const Triangle triangle,
- const Transpose a_transpose, const Diagonal diagonal,
- const size_t n, const size_t k,
- 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) {
+void Xtbmv<T>::DoTbmv(const Layout layout, const Triangle triangle,
+ const Transpose a_transpose, const Diagonal diagonal,
+ const size_t n, const size_t k,
+ 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) {
// Creates a copy of X: a temporary scratch buffer
auto scratch_buffer = Buffer<T>(context_, n*x_inc + x_offset);
- try {
- x_buffer.CopyTo(queue_, n*x_inc + x_offset, scratch_buffer);
- } catch (...) { } // Continues: error-code is returned in MatVec
+ x_buffer.CopyTo(queue_, n*x_inc + x_offset, scratch_buffer);
// The data is either in the upper or lower triangle
size_t is_upper = ((triangle == Triangle::kUpper && layout != Layout::kRowMajor) ||
@@ -52,20 +50,22 @@ StatusCode Xtbmv<T>::DoTbmv(const Layout layout, const Triangle triangle,
// The specific triangular banded matrix-accesses are implemented in the kernel guarded by the
// ROUTINE_TBMV define.
auto fast_kernels = false;
- auto status = MatVec(layout, a_transpose,
- n, n, static_cast<T>(1),
- a_buffer, a_offset, a_ld,
- scratch_buffer, x_offset, x_inc, static_cast<T>(0),
- x_buffer, x_offset, x_inc,
- fast_kernels, fast_kernels,
- parameter, false, k, 0);
-
- // Returns the proper error code (renames vector Y to X)
- switch(status) {
- case StatusCode::kInvalidVectorY: return StatusCode::kInvalidVectorX;
- case StatusCode::kInvalidIncrementY: return StatusCode::kInvalidIncrementX;
- case StatusCode::kInsufficientMemoryY: return StatusCode::kInsufficientMemoryX;
- default: return status;
+ try {
+ MatVec(layout, a_transpose,
+ n, n, static_cast<T>(1),
+ a_buffer, a_offset, a_ld,
+ scratch_buffer, x_offset, x_inc, static_cast<T>(0),
+ x_buffer, x_offset, x_inc,
+ fast_kernels, fast_kernels,
+ parameter, false, k, 0);
+ } catch (BLASError &e) {
+ // Returns the proper error code (renames vector Y to X)
+ switch (e.status()) {
+ case StatusCode::kInvalidVectorY: throw BLASError(StatusCode::kInvalidVectorX, e.details());
+ case StatusCode::kInvalidIncrementY: throw BLASError(StatusCode::kInvalidIncrementX, e.details());
+ case StatusCode::kInsufficientMemoryY: throw BLASError(StatusCode::kInsufficientMemoryX, e.details());
+ default: throw;
+ }
}
}