Added script to generate API interface and implementation automatically

1 files changed, 267 insertions, 270 deletions

diff --git a/src/clblast.cc b/src/clblast.cc
index 3303085e..0ced9ff7 100644
--- a/src/clblast.cc
+++ b/src/clblast.cc
@@ -42,11 +42,12 @@
 #include "internal/routines/level3/xtrmm.h"
 
 namespace clblast {
+
 // =================================================================================================
 // BLAS level-1 (vector-vector) routines
 // =================================================================================================
 
-// SWAP
+// Swap two vectors: SSWAP/DSWAP/CSWAP/ZSWAP
 template <typename T>
 StatusCode Swap(const size_t n,
                 cl_mem x_buffer, const size_t x_offset, const size_t x_inc,
@@ -64,21 +65,21 @@ StatusCode Swap(const size_t n,
 template StatusCode Swap<float>(const size_t,
                                 cl_mem, const size_t, const size_t,
                                 cl_mem, const size_t, const size_t,
-                                cl_command_queue* queue, cl_event* event);
+                                cl_command_queue*, cl_event*);
 template StatusCode Swap<double>(const size_t,
                                  cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Swap<float2>(const size_t,
                                  cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Swap<double2>(const size_t,
                                   cl_mem, const size_t, const size_t,
                                   cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
-// SCAL
+// Vector scaling: SSCAL/DSCAL/CSCAL/ZSCAL
 template <typename T>
 StatusCode Scal(const size_t n,
                 const T alpha,
@@ -96,21 +97,21 @@ StatusCode Scal(const size_t n,
 template StatusCode Scal<float>(const size_t,
                                 const float,
                                 cl_mem, const size_t, const size_t,
-                                cl_command_queue* queue, cl_event* event);
+                                cl_command_queue*, cl_event*);
 template StatusCode Scal<double>(const size_t,
                                  const double,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Scal<float2>(const size_t,
                                  const float2,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Scal<double2>(const size_t,
                                   const double2,
                                   cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
-// COPY
+// Vector copy: SCOPY/DCOPY/CCOPY/ZCOPY
 template <typename T>
 StatusCode Copy(const size_t n,
                 const cl_mem x_buffer, const size_t x_offset, const size_t x_inc,
@@ -128,21 +129,21 @@ StatusCode Copy(const size_t n,
 template StatusCode Copy<float>(const size_t,
                                 const cl_mem, const size_t, const size_t,
                                 cl_mem, const size_t, const size_t,
-                                cl_command_queue* queue, cl_event* event);
+                                cl_command_queue*, cl_event*);
 template StatusCode Copy<double>(const size_t,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Copy<float2>(const size_t,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Copy<double2>(const size_t,
                                   const cl_mem, const size_t, const size_t,
                                   cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
-// AXPY
+// Vector-times-constant plus vector: SAXPY/DAXPY/CAXPY/ZAXPY
 template <typename T>
 StatusCode Axpy(const size_t n,
                 const T alpha,
@@ -163,24 +164,24 @@ template StatusCode Axpy<float>(const size_t,
                                 const float,
                                 const cl_mem, const size_t, const size_t,
                                 cl_mem, const size_t, const size_t,
-                                cl_command_queue* queue, cl_event* event);
+                                cl_command_queue*, cl_event*);
 template StatusCode Axpy<double>(const size_t,
                                  const double,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Axpy<float2>(const size_t,
                                  const float2,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Axpy<double2>(const size_t,
                                   const double2,
                                   const cl_mem, const size_t, const size_t,
                                   cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
-// DOT
+// Dot product of two vectors: SDOT/DDOT
 template <typename T>
 StatusCode Dot(const size_t n,
                cl_mem dot_buffer, const size_t dot_offset,
@@ -201,14 +202,14 @@ template StatusCode Dot<float>(const size_t,
                                cl_mem, const size_t,
                                const cl_mem, const size_t, const size_t,
                                const cl_mem, const size_t, const size_t,
-                               cl_command_queue* queue, cl_event* event);
+                               cl_command_queue*, cl_event*);
 template StatusCode Dot<double>(const size_t,
                                 cl_mem, const size_t,
                                 const cl_mem, const size_t, const size_t,
                                 const cl_mem, const size_t, const size_t,
-                                cl_command_queue* queue, cl_event* event);
+                                cl_command_queue*, cl_event*);
 
-// DOTU
+// Dot product of two complex vectors: CDOTU/ZDOTU
 template <typename T>
 StatusCode Dotu(const size_t n,
                 cl_mem dot_buffer, const size_t dot_offset,
@@ -229,14 +230,14 @@ template StatusCode Dotu<float2>(const size_t,
                                  cl_mem, const size_t,
                                  const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Dotu<double2>(const size_t,
                                   cl_mem, const size_t,
                                   const cl_mem, const size_t, const size_t,
                                   const cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
-// DOTC
+// Dot product of two complex vectors, one conjugated: CDOTC/ZDOTC
 template <typename T>
 StatusCode Dotc(const size_t n,
                 cl_mem dot_buffer, const size_t dot_offset,
@@ -257,138 +258,150 @@ template StatusCode Dotc<float2>(const size_t,
                                  cl_mem, const size_t,
                                  const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 cl_command_queue* queue, cl_event* event);
+                                 cl_command_queue*, cl_event*);
 template StatusCode Dotc<double2>(const size_t,
                                   cl_mem, const size_t,
                                   const cl_mem, const size_t, const size_t,
                                   const cl_mem, const size_t, const size_t,
-                                  cl_command_queue* queue, cl_event* event);
+                                  cl_command_queue*, cl_event*);
 
 // =================================================================================================
 // BLAS level-2 (matrix-vector) routines
 // =================================================================================================
 
-// GEMV
+// Generalized matrix-vector multiplication: SGEMV/DGEMV/CGEMV/ZGEMV
 template <typename T>
 StatusCode Gemv(const Layout layout, const Transpose a_transpose,
-                const size_t m, const size_t n, const T alpha,
+                const size_t m, 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,
+                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 = Xgemv<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.DoGemv(layout, a_transpose, m, n, alpha,
+  return routine.DoGemv(layout, a_transpose,
+                        m, n,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(x_buffer), x_offset, x_inc, beta,
+                        Buffer<T>(x_buffer), x_offset, x_inc,
+                        beta,
                         Buffer<T>(y_buffer), y_offset, y_inc);
 }
 template StatusCode Gemv<float>(const Layout, const Transpose,
-                                const size_t, const size_t, const float,
+                                const size_t, 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,
+                                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 Gemv<double>(const Layout, const Transpose,
-                                 const size_t, const size_t, const double,
+                                 const size_t, 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,
+                                 const cl_mem, const size_t, const size_t,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Gemv<float2>(const Layout, const Transpose,
-                                 const size_t, const size_t, const float2,
+                                 const size_t, const size_t,
+                                 const float2,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const cl_mem, const size_t, const size_t,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Gemv<double2>(const Layout, const Transpose,
-                                  const size_t, const size_t, const double2,
+                                  const size_t, const size_t,
+                                  const double2,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const cl_mem, const size_t, const size_t,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// HEMV
+// Hermitian matrix-vector multiplication: CHEMV/ZHEMV
 template <typename T>
 StatusCode Hemv(const Layout layout, const Triangle triangle,
-                const size_t n, const T alpha,
+                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,
+                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 = Xhemv<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.DoHemv(layout, triangle, n, alpha,
+  return routine.DoHemv(layout, triangle,
+                        n,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(x_buffer), x_offset, x_inc, beta,
+                        Buffer<T>(x_buffer), x_offset, x_inc,
+                        beta,
                         Buffer<T>(y_buffer), y_offset, y_inc);
 }
 template StatusCode Hemv<float2>(const Layout, const Triangle,
-                                 const size_t, const float2,
+                                 const size_t,
+                                 const float2,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const cl_mem, const size_t, const size_t,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Hemv<double2>(const Layout, const Triangle,
-                                  const size_t, const double2,
+                                  const size_t,
+                                  const double2,
+                                  const cl_mem, const size_t, const size_t,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// SYMV
+// Symmetric matrix-vector multiplication: SSYMV/DSYMV
 template <typename T>
 StatusCode Symv(const Layout layout, const Triangle triangle,
-                const size_t n, const T alpha,
+                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,
+                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,
+  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>(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 size_t,
+                                const float,
                                 const cl_mem, const size_t, const size_t,
-                                const cl_mem, const size_t, const size_t, const float,
+                                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 size_t,
+                                 const double,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const double,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 
@@ -396,307 +409,343 @@ template StatusCode Symv<double>(const Layout, const Triangle,
 // BLAS level-3 (matrix-matrix) routines
 // =================================================================================================
 
-// GEMM
+// Generalized matrix-matrix multiplication: SGEMM/DGEMM/CGEMM/ZGEMM
 template <typename T>
 StatusCode Gemm(const Layout layout, const Transpose a_transpose, const Transpose b_transpose,
-                const size_t m, const size_t n, const size_t k, const T alpha,
+                const size_t m, const size_t n, const size_t k,
+                const T alpha,
                 const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
-                const cl_mem b_buffer, const size_t b_offset, const size_t b_ld, const T beta,
+                const cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
+                const T beta,
                 cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                 cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xgemm<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.DoGemm(layout, a_transpose, b_transpose, m, n, k, alpha,
+  return routine.DoGemm(layout, a_transpose, b_transpose,
+                        m, n, k,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(b_buffer), b_offset, b_ld, beta,
+                        Buffer<T>(b_buffer), b_offset, b_ld,
+                        beta,
                         Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Gemm<float>(const Layout, const Transpose, const Transpose,
-                                const size_t, const size_t, const size_t, const float,
+                                const size_t, const size_t, const size_t,
+                                const float,
+                                const cl_mem, const size_t, const size_t,
                                 const cl_mem, const size_t, const size_t,
-                                const cl_mem, const size_t, const size_t, const float,
+                                const float,
                                 cl_mem, const size_t, const size_t,
                                 cl_command_queue*, cl_event*);
 template StatusCode Gemm<double>(const Layout, const Transpose, const Transpose,
-                                 const size_t, const size_t, const size_t, const double,
+                                 const size_t, const size_t, 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,
+                                 const cl_mem, const size_t, const size_t,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Gemm<float2>(const Layout, const Transpose, const Transpose,
-                                 const size_t, const size_t, const size_t, const float2,
+                                 const size_t, const size_t, const size_t,
+                                 const float2,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Gemm<double2>(const Layout, const Transpose, const Transpose,
-                                  const size_t, const size_t, const size_t, const double2,
+                                  const size_t, const size_t, const size_t,
+                                  const double2,
+                                  const cl_mem, const size_t, const size_t,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// SYMM
+// Symmetric matrix-matrix multiplication: SSYMM/DSYMM/CSYMM/ZSYMM
 template <typename T>
 StatusCode Symm(const Layout layout, const Side side, const Triangle triangle,
-                const size_t m, const size_t n, const T alpha,
+                const size_t m, 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 b_buffer, const size_t b_offset, const size_t b_ld, const T beta,
+                const cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
+                const T beta,
                 cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                 cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xsymm<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.DoSymm(layout, side, triangle, m, n, alpha,
+  return routine.DoSymm(layout, side, triangle,
+                        m, n,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(b_buffer), b_offset, b_ld, beta,
+                        Buffer<T>(b_buffer), b_offset, b_ld,
+                        beta,
                         Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Symm<float>(const Layout, const Side, const Triangle,
-                                const size_t, const size_t, const float,
+                                const size_t, const size_t,
+                                const float,
+                                const cl_mem, const size_t, const size_t,
                                 const cl_mem, const size_t, const size_t,
-                                const cl_mem, const size_t, const size_t, const float,
+                                const float,
                                 cl_mem, const size_t, const size_t,
                                 cl_command_queue*, cl_event*);
 template StatusCode Symm<double>(const Layout, const Side, const Triangle,
-                                 const size_t, const size_t, const double,
+                                 const size_t, const size_t,
+                                 const double,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const double,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Symm<float2>(const Layout, const Side, const Triangle,
-                                 const size_t, const size_t, const float2,
+                                 const size_t, const size_t,
+                                 const float2,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Symm<double2>(const Layout, const Side, const Triangle,
-                                  const size_t, const size_t, const double2,
+                                  const size_t, const size_t,
+                                  const double2,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const cl_mem, const size_t, const size_t,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// HEMM
+// Hermitian matrix-matrix multiplication: CHEMM/ZHEMM
 template <typename T>
 StatusCode Hemm(const Layout layout, const Side side, const Triangle triangle,
-                const size_t m, const size_t n, const T alpha,
+                const size_t m, 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 b_buffer, const size_t b_offset, const size_t b_ld, const T beta,
+                const cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
+                const T beta,
                 cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                 cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xhemm<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.DoHemm(layout, side, triangle, m, n, alpha,
+  return routine.DoHemm(layout, side, triangle,
+                        m, n,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(b_buffer), b_offset, b_ld, beta,
+                        Buffer<T>(b_buffer), b_offset, b_ld,
+                        beta,
                         Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Hemm<float2>(const Layout, const Side, const Triangle,
-                                 const size_t, const size_t, const float2,
+                                 const size_t, const size_t,
+                                 const float2,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Hemm<double2>(const Layout, const Side, const Triangle,
-                                  const size_t, const size_t, const double2,
+                                  const size_t, const size_t,
+                                  const double2,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const cl_mem, const size_t, const size_t,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// SYRK
+// Rank-K update of a symmetric matrix: SSYRK/DSYRK/CSYRK/ZSYRK
 template <typename T>
 StatusCode Syrk(const Layout layout, const Triangle triangle, const Transpose a_transpose,
-                const size_t n, const size_t k, const T alpha,
-                const cl_mem a_buffer, const size_t a_offset, const size_t a_ld, const T beta,
+                const size_t n, const size_t k,
+                const T alpha,
+                const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
+                const T beta,
                 cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                 cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xsyrk<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.DoSyrk(layout, triangle, a_transpose, n, k, alpha,
-                        Buffer<T>(a_buffer), a_offset, a_ld, beta,
+  return routine.DoSyrk(layout, triangle, a_transpose,
+                        n, k,
+                        alpha,
+                        Buffer<T>(a_buffer), a_offset, a_ld,
+                        beta,
                         Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Syrk<float>(const Layout, const Triangle, const Transpose,
-                                const size_t, const size_t, const float,
-                                const cl_mem, const size_t, const size_t, const float,
+                                const size_t, const size_t,
+                                const float,
+                                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 Syrk<double>(const Layout, const Triangle, const Transpose,
-                                 const size_t, const size_t, const double,
-                                 const cl_mem, const size_t, const size_t, const double,
+                                 const size_t, const size_t,
+                                 const double,
+                                 const cl_mem, const size_t, const size_t,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Syrk<float2>(const Layout, const Triangle, const Transpose,
-                                 const size_t, const size_t, const float2,
-                                 const cl_mem, const size_t, const size_t, const float2,
+                                 const size_t, const size_t,
+                                 const float2,
+                                 const cl_mem, const size_t, const size_t,
+                                 const float2,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Syrk<double2>(const Layout, const Triangle, const Transpose,
-                                  const size_t, const size_t, const double2,
-                                  const cl_mem, const size_t, const size_t, const double2,
+                                  const size_t, const size_t,
+                                  const double2,
+                                  const cl_mem, const size_t, const size_t,
+                                  const double2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// HERK
+// Rank-K update of a hermitian matrix: CHERK/ZHERK
 template <typename T>
 StatusCode Herk(const Layout layout, const Triangle triangle, const Transpose a_transpose,
-                const size_t n, const size_t k, const T alpha,
-                const cl_mem a_buffer, const size_t a_offset, const size_t a_ld, const T beta,
+                const size_t n, const size_t k,
+                const T alpha,
+                const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
+                const T beta,
                 cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                 cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xherk<std::complex<T>,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.DoHerk(layout, triangle, a_transpose, n, k, alpha,
-                        Buffer<std::complex<T>>(a_buffer), a_offset, a_ld, beta,
+  return routine.DoHerk(layout, triangle, a_transpose,
+                        n, k,
+                        alpha,
+                        Buffer<std::complex<T>>(a_buffer), a_offset, a_ld,
+                        beta,
                         Buffer<std::complex<T>>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Herk<float>(const Layout, const Triangle, const Transpose,
-                                const size_t, const size_t, const float,
-                                const cl_mem, const size_t, const size_t, const float,
+                                const size_t, const size_t,
+                                const float,
+                                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 Herk<double>(const Layout, const Triangle, const Transpose,
-                                 const size_t, const size_t, const double,
-                                 const cl_mem, const size_t, const size_t, const double,
+                                 const size_t, const size_t,
+                                 const double,
+                                 const cl_mem, const size_t, const size_t,
+                                 const double,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// SYR2K
+// Rank-2K update of a symmetric matrix: SSYR2K/DSYR2K/CSYR2K/ZSYR2K
 template <typename T>
 StatusCode Syr2k(const Layout layout, const Triangle triangle, const Transpose ab_transpose,
-                 const size_t n, const size_t k, const T alpha,
+                 const size_t n, const size_t k,
+                 const T alpha,
                  const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
-                 const cl_mem b_buffer, const size_t b_offset, const size_t b_ld, const T beta,
+                 const cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
+                 const T beta,
                  cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                  cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xsyr2k<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.DoSyr2k(layout, triangle, ab_transpose, n, k, alpha,
+  return routine.DoSyr2k(layout, triangle, ab_transpose,
+                         n, k,
+                         alpha,
                          Buffer<T>(a_buffer), a_offset, a_ld,
-                         Buffer<T>(b_buffer), b_offset, b_ld, beta,
+                         Buffer<T>(b_buffer), b_offset, b_ld,
+                         beta,
                          Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Syr2k<float>(const Layout, const Triangle, const Transpose,
-                                 const size_t, const size_t, const float,
+                                 const size_t, const size_t,
+                                 const float,
+                                 const cl_mem, const size_t, const size_t,
                                  const cl_mem, const size_t, const size_t,
-                                 const cl_mem, const size_t, const size_t, const float,
+                                 const float,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
 template StatusCode Syr2k<double>(const Layout, const Triangle, const Transpose,
-                                  const size_t, const size_t, const double,
+                                  const size_t, 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,
+                                  const cl_mem, const size_t, const size_t,
+                                  const double,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 template StatusCode Syr2k<float2>(const Layout, const Triangle, const Transpose,
-                                  const size_t, const size_t, const float2,
+                                  const size_t, const size_t,
+                                  const float2,
+                                  const cl_mem, const size_t, const size_t,
                                   const cl_mem, const size_t, const size_t,
-                                  const cl_mem, const size_t, const size_t, const float2,
+                                  const float2,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 template StatusCode Syr2k<double2>(const Layout, const Triangle, const Transpose,
-                                   const size_t, const size_t, const double2,
+                                   const size_t, const size_t,
+                                   const double2,
+                                   const cl_mem, const size_t, const size_t,
                                    const cl_mem, const size_t, const size_t,
-                                   const cl_mem, const size_t, const size_t, const double2,
+                                   const double2,
                                    cl_mem, const size_t, const size_t,
                                    cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// SYR2K
+// Rank-2K update of a hermitian matrix: CHER2K/ZHER2K
 template <typename T, typename U>
 StatusCode Her2k(const Layout layout, const Triangle triangle, const Transpose ab_transpose,
-                 const size_t n, const size_t k, const T alpha,
+                 const size_t n, const size_t k,
+                 const T alpha,
                  const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
-                 const cl_mem b_buffer, const size_t b_offset, const size_t b_ld, const U beta,
+                 const cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
+                 const U beta,
                  cl_mem c_buffer, const size_t c_offset, const size_t c_ld,
                  cl_command_queue* queue, cl_event* event) {
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xher2k<T,U>(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.DoHer2k(layout, triangle, ab_transpose, n, k, alpha,
+  return routine.DoHer2k(layout, triangle, ab_transpose,
+                         n, k,
+                         alpha,
                          Buffer<T>(a_buffer), a_offset, a_ld,
-                         Buffer<T>(b_buffer), b_offset, b_ld, beta,
+                         Buffer<T>(b_buffer), b_offset, b_ld,
+                         beta,
                          Buffer<T>(c_buffer), c_offset, c_ld);
 }
 template StatusCode Her2k<float2,float>(const Layout, const Triangle, const Transpose,
-                                        const size_t, const size_t, const float2,
+                                        const size_t, const size_t,
+                                        const float2,
+                                        const cl_mem, const size_t, const size_t,
                                         const cl_mem, const size_t, const size_t,
-                                        const cl_mem, const size_t, const size_t, const float,
+                                        const float,
                                         cl_mem, const size_t, const size_t,
                                         cl_command_queue*, cl_event*);
 template StatusCode Her2k<double2,double>(const Layout, const Triangle, const Transpose,
-                                          const size_t, const size_t, const double2,
+                                          const size_t, const size_t,
+                                          const double2,
+                                          const cl_mem, const size_t, const size_t,
                                           const cl_mem, const size_t, const size_t,
-                                          const cl_mem, const size_t, const size_t, const double,
+                                          const double,
                                           cl_mem, const size_t, const size_t,
                                           cl_command_queue*, cl_event*);
 
-// =================================================================================================
-
-// TRMM
+// Triangular matrix-matrix multiplication: STRMM/DTRMM/CTRMM/ZTRMM
 template <typename T>
-StatusCode Trmm(const Layout layout, const Side side, const Triangle triangle,
-                const Transpose a_transpose, const Diagonal diagonal,
+StatusCode Trmm(const Layout layout, const Side side, const Triangle triangle, const Transpose a_transpose, const Diagonal diagonal,
                 const size_t m, const size_t n,
                 const T alpha,
                 const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
@@ -705,90 +754,38 @@ StatusCode Trmm(const Layout layout, const Side side, const Triangle triangle,
   auto queue_cpp = Queue(*queue);
   auto event_cpp = Event(*event);
   auto routine = Xtrmm<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.DoTrmm(layout, side, triangle, a_transpose, diagonal, m, n, alpha,
+  return routine.DoTrmm(layout, side, triangle, a_transpose, diagonal,
+                        m, n,
+                        alpha,
                         Buffer<T>(a_buffer), a_offset, a_ld,
                         Buffer<T>(b_buffer), b_offset, b_ld);
 }
-template StatusCode Trmm<float>(const Layout, const Side, const Triangle,
-                                const Transpose, const Diagonal,
-                                const size_t, const size_t, const float,
+template StatusCode Trmm<float>(const Layout, const Side, const Triangle, const Transpose, const Diagonal,
+                                const size_t, const size_t,
+                                const float,
                                 const cl_mem, const size_t, const size_t,
                                 cl_mem, const size_t, const size_t,
                                 cl_command_queue*, cl_event*);
-template StatusCode Trmm<double>(const Layout, const Side, const Triangle,
-                                 const Transpose, const Diagonal,
-                                 const size_t, const size_t, const double,
+template StatusCode Trmm<double>(const Layout, const Side, const Triangle, const Transpose, const Diagonal,
+                                 const size_t, const size_t,
+                                 const double,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
-template StatusCode Trmm<float2>(const Layout, const Side, const Triangle,
-                                 const Transpose, const Diagonal,
-                                 const size_t, const size_t, const float2,
+template StatusCode Trmm<float2>(const Layout, const Side, const Triangle, const Transpose, const Diagonal,
+                                 const size_t, const size_t,
+                                 const float2,
                                  const cl_mem, const size_t, const size_t,
                                  cl_mem, const size_t, const size_t,
                                  cl_command_queue*, cl_event*);
-template StatusCode Trmm<double2>(const Layout, const Side, const Triangle,
-                                  const Transpose, const Diagonal,
-                                  const size_t, const size_t, const double2,
+template StatusCode Trmm<double2>(const Layout, const Side, const Triangle, const Transpose, const Diagonal,
+                                  const size_t, const size_t,
+                                  const double2,
                                   const cl_mem, const size_t, const size_t,
                                   cl_mem, const size_t, const size_t,
                                   cl_command_queue*, cl_event*);
 
 // =================================================================================================
-
-// TRSM
-/*
-template <typename T>
-StatusCode Trsm(const Layout layout, const Side side, const Triangle triangle,
-                const Transpose a_transpose, const Diagonal diagonal,
-                const size_t m, const size_t n,
-                const T alpha,
-                const cl_mem a_buffer, const size_t a_offset, const size_t a_ld,
-                cl_mem b_buffer, const size_t b_offset, const size_t b_ld,
-                cl_command_queue* queue, cl_event* event) {
-  auto queue_cpp = Queue(*queue);
-  auto event_cpp = Event(*event);
-  auto routine = Xtrsm<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.DoTrsm(layout, side, triangle, a_transpose, diagonal, m, n, alpha,
-                        Buffer<T>(a_buffer), a_offset, a_ld,
-                        Buffer<T>(b_buffer), b_offset, b_ld);
-}
-template StatusCode Trsm<float>(const Layout, const Side, const Triangle,
-                                const Transpose, const Diagonal,
-                                const size_t, const size_t, const float,
-                                const cl_mem, const size_t, const size_t,
-                                cl_mem, const size_t, const size_t,
-                                cl_command_queue*, cl_event*);
-template StatusCode Trsm<double>(const Layout, const Side, const Triangle,
-                                 const Transpose, const Diagonal,
-                                 const size_t, const size_t, const double,
-                                 const cl_mem, const size_t, const size_t,
-                                 cl_mem, const size_t, const size_t,
-                                 cl_command_queue*, cl_event*);
-template StatusCode Trsm<float2>(const Layout, const Side, const Triangle,
-                                 const Transpose, const Diagonal,
-                                 const size_t, const size_t, const float2,
-                                 const cl_mem, const size_t, const size_t,
-                                 cl_mem, const size_t, const size_t,
-                                 cl_command_queue*, cl_event*);
-template StatusCode Trsm<double2>(const Layout, const Side, const Triangle,
-                                  const Transpose, const Diagonal,
-                                  const size_t, const size_t, const double2,
-                                  const cl_mem, const size_t, const size_t,
-                                  cl_mem, const size_t, const size_t,
-                                  cl_command_queue*, cl_event*);
-*/
-// =================================================================================================
 } // namespace clblast