CLBlast/test/wrapper_clblas.h

// =================================================================================================
// This file is part of the CLBlast project. The project is licensed under the MIT license. This
// project loosely follows the Google C++ styleguide and uses a tab-size of two spaces and a max-
// width of 100 characters per line.
//
// Author(s):
//   Cedric Nugteren <www.cedricnugteren.nl>
//
// This file implements a wrapper around the clBLAS library, such that its routines can be called
// in a similar way as the CLBlast routines: using alpha and beta to determine the precision.
//
// =================================================================================================

#ifndef CLBLAST_TEST_WRAPPER_CLBLAS_H_
#define CLBLAST_TEST_WRAPPER_CLBLAS_H_

#include <clBLAS.h>

#include "internal/utilities.h"

namespace clblast {
// =================================================================================================
// BLAS level-1 (vector-vector) routines

// Calls {clblasSaxpy, clblasDaxpy, clblasCaxpy, clblasZaxpy} with the arguments forwarded.
clblasStatus clblasXaxpy(
  size_t n, float alpha,
  const cl_mem x_vec, size_t x_offset, size_t x_inc,
  const cl_mem y_vec, size_t y_offset, size_t y_inc,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasSaxpy(n, alpha,
                       x_vec, x_offset, static_cast<int>(x_inc),
                       y_vec, y_offset, static_cast<int>(y_inc),
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXaxpy(
  size_t n, double alpha,
  const cl_mem x_vec, size_t x_offset, size_t x_inc,
  const cl_mem y_vec, size_t y_offset, size_t y_inc,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasDaxpy(n, alpha,
                       x_vec, x_offset, static_cast<int>(x_inc),
                       y_vec, y_offset, static_cast<int>(y_inc),
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXaxpy(
  size_t n, float2 alpha,
  const cl_mem x_vec, size_t x_offset, size_t x_inc,
  const cl_mem y_vec, size_t y_offset, size_t y_inc,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_float2{{alpha.real(), alpha.imag()}};
    return clblasCaxpy(n, cl_alpha,
                       x_vec, x_offset, static_cast<int>(x_inc),
                       y_vec, y_offset, static_cast<int>(y_inc),
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXaxpy(
  size_t n, double2 alpha,
  const cl_mem x_vec, size_t x_offset, size_t x_inc,
  const cl_mem y_vec, size_t y_offset, size_t y_inc,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_double2{{alpha.real(), alpha.imag()}};
    return clblasZaxpy(n, cl_alpha,
                       x_vec, x_offset, static_cast<int>(x_inc),
                       y_vec, y_offset, static_cast<int>(y_inc),
                       num_queues, queues, num_wait_events, wait_events, events);
}

// =================================================================================================
// BLAS level-2 (matrix-vector) routines

// =================================================================================================
// BLAS level-3 (matrix-matrix) routines

// This calls {clblasSgemm, clblasDgemm, clblasCgemm, clblasZgemm} with the arguments forwarded.
clblasStatus clblasXgemm(
  clblasOrder layout, clblasTranspose tran_a, clblasTranspose tran_b,
  size_t m, size_t n, size_t k, float alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, float beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasSgemm(layout, tran_a, tran_b,
                       m, n, k, alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXgemm(
  clblasOrder layout, clblasTranspose tran_a, clblasTranspose tran_b,
  size_t m, size_t n, size_t k, double alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, double beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasDgemm(layout, tran_a, tran_b,
                       m, n, k, alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXgemm(
  clblasOrder layout, clblasTranspose tran_a, clblasTranspose tran_b,
  size_t m, size_t n, size_t k, float2 alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, float2 beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_float2{{alpha.real(), alpha.imag()}};
    auto cl_beta = cl_float2{{beta.real(), beta.imag()}};
    return clblasCgemm(layout, tran_a, tran_b,
                       m, n, k, cl_alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, cl_beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXgemm(
  clblasOrder layout, clblasTranspose tran_a, clblasTranspose tran_b,
  size_t m, size_t n, size_t k, double2 alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, double2 beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_double2{{alpha.real(), alpha.imag()}};
    auto cl_beta = cl_double2{{beta.real(), beta.imag()}};
    return clblasZgemm(layout, tran_a, tran_b,
                       m, n, k, cl_alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, cl_beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}

// This calls {clblasSsymm, clblasDsymm} with the arguments forwarded.
clblasStatus clblasXsymm(
  clblasOrder layout, clblasSide side, clblasUplo triangle,
  size_t m, size_t n, float alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, float beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasSsymm(layout, side, triangle,
                       m, n, alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXsymm(
  clblasOrder layout, clblasSide side, clblasUplo triangle,
  size_t m, size_t n, double alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, double beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    return clblasDsymm(layout, side, triangle,
                       m, n, alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXsymm(
  clblasOrder layout, clblasSide side, clblasUplo triangle,
  size_t m, size_t n, float2 alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, float2 beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_float2{{alpha.real(), alpha.imag()}};
    auto cl_beta = cl_float2{{beta.real(), beta.imag()}};
    return clblasCsymm(layout, side, triangle,
                       m, n, cl_alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, cl_beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}
clblasStatus clblasXsymm(
  clblasOrder layout, clblasSide side, clblasUplo triangle,
  size_t m, size_t n, double2 alpha,
  const cl_mem a_mat, size_t a_offset, size_t a_ld,
  const cl_mem b_mat, size_t b_offset, size_t b_ld, double2 beta,
  cl_mem c_mat, size_t c_offset, size_t c_ld,
  cl_uint num_queues, cl_command_queue *queues,
  cl_uint num_wait_events, const cl_event *wait_events, cl_event *events) {
    auto cl_alpha = cl_double2{{alpha.real(), alpha.imag()}};
    auto cl_beta = cl_double2{{beta.real(), beta.imag()}};
    return clblasZsymm(layout, side, triangle,
                       m, n, cl_alpha,
                       a_mat, a_offset, a_ld,
                       b_mat, b_offset, b_ld, cl_beta,
                       c_mat, c_offset, c_ld,
                       num_queues, queues, num_wait_events, wait_events, events);
}

// =================================================================================================
} // namespace clblast

// CLBLAST_TEST_WRAPPER_CLBLAS_H_
#endif