dotu.hpp

//*****************************************************************************
// Copyright (c) 2024, Intel Corporation
// All rights reserved.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are met:
// - Redistributions of source code must retain the above copyright notice,
//   this list of conditions and the following disclaimer.
// - Redistributions in binary form must reproduce the above copyright notice,
//   this list of conditions and the following disclaimer in the documentation
//   and/or other materials provided with the distribution.
// - Neither the name of the copyright holder nor the names of its contributors
//   may be used to endorse or promote products derived from this software
//   without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
// AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
// ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
// LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
// CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
// SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
// INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
// CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
// ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
// THE POSSIBILITY OF SUCH DAMAGE.
//*****************************************************************************
 
#pragma once
 
#include <stdexcept>
 
#include <pybind11/pybind11.h>
 
#include "dot_common.hpp"
 
namespace dpnp::extensions::blas
{
namespace mkl_blas = oneapi::mkl::blas;
namespace type_utils = dpnp::tensor::type_utils;
namespace py = pybind11;
 
template <typename T>
static sycl::event dotu_impl(sycl::queue &exec_q,
                             const std::int64_t n,
                             const char *vectorX,
                             const std::int64_t incx,
                             const char *vectorY,
                             const std::int64_t incy,
                             char *result,
                             const std::vector<sycl::event> &depends)
{
    type_utils::validate_type_for_device<T>(exec_q);
 
    const T *x = reinterpret_cast<const T *>(vectorX);
    const T *y = reinterpret_cast<const T *>(vectorY);
    T *res = reinterpret_cast<T *>(result);
 
    std::stringstream error_msg;
    bool is_exception_caught = false;
 
    sycl::event dotu_event;
    try {
        // Release GIL to avoid serialization of host task submissions
        // to the same queue in OneMKL
        py::gil_scoped_release lock{};
 
        dotu_event =
            mkl_blas::column_major::dotu(exec_q,
                                         n,    // size of the input vectors
                                         x,    // Pointer to vector x.
                                         incx, // Stride of vector x.
                                         y,    // Pointer to vector y.
                                         incy, // Stride of vector y.
                                         res,  // Pointer to result.
                                         depends);
    } catch (oneapi::mkl::exception const &e) {
        error_msg
            << "Unexpected MKL exception caught during dotu() call:\nreason: "
            << e.what();
        is_exception_caught = true;
    } catch (sycl::exception const &e) {
        error_msg << "Unexpected SYCL exception caught during dotu() call:\n"
                  << e.what();
        is_exception_caught = true;
    }
 
    if (is_exception_caught) // an unexpected error occurs
    {
        throw std::runtime_error(error_msg.str());
    }
 
    return dotu_event;
}
 
template <typename fnT, typename varT>
struct DotuContigFactory
{
    fnT get()
    {
        if constexpr (types::DotuTypePairSupportFactory<varT>::is_defined) {
            return dotu_impl<varT>;
        }
        else {
            return nullptr;
        }
    }
};
} // namespace dpnp::extensions::blas