/* ************************************************************************
 * Copyright (C) 2016-2023 Advanced Micro Devices, Inc. All rights reserved.
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to deal
 * in the Software without restriction, including without limitation the rights
 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell cop-
 * ies of the Software, and to permit persons to whom the Software is furnished
 * to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in all
 * copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IM-
 * PLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS
 * FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
 * IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNE-
 * CTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
 *
 * ************************************************************************ */

#include "check_numerics_vector.hpp"
#include "handle.hpp"
#include "rocblas_copy.hpp"

template <typename T, typename U>
ROCBLAS_KERNEL_NO_BOUNDS rocblas_copy_kernel(rocblas_int    n,
                                             const T        xa,
                                             rocblas_stride shiftx,
                                             rocblas_int    incx,
                                             rocblas_stride stridex,
                                             U              ya,
                                             rocblas_stride shifty,
                                             rocblas_int    incy,
                                             rocblas_stride stridey)
{
    int64_t     tid = blockIdx.x * blockDim.x + threadIdx.x;
    const auto* x   = load_ptr_batch(xa, blockIdx.y, shiftx, stridex);
    auto*       y   = load_ptr_batch(ya, blockIdx.y, shifty, stridey);
    if(tid < n)
    {

        y[tid * incy] = x[tid * incx];
    }
}

//! @brief Optimized kernel for the floating points.
//!
template <rocblas_int NB, typename T, typename U>
ROCBLAS_KERNEL(NB)
rocblas_scopy_2_kernel(rocblas_int n,
                       const T __restrict xa,
                       rocblas_stride shiftx,
                       rocblas_stride stridex,
                       U __restrict ya,
                       rocblas_stride shifty,
                       rocblas_stride stridey)
{
    int64_t     tid = (blockIdx.x * blockDim.x + threadIdx.x) * 2;
    const auto* x   = load_ptr_batch(xa, blockIdx.y, shiftx, stridex);
    auto*       y   = load_ptr_batch(ya, blockIdx.y, shifty, stridey);
    if(tid < n - 1)
    {
        for(int j = 0; j < 2; ++j)
        {
            y[tid + j] = x[tid + j];
        }
    }
    if(n % 2 != 0 && tid == n - 1)
        y[tid] = x[tid];
}

template <rocblas_int NB, typename T, typename U>
rocblas_status rocblas_copy_template(rocblas_handle handle,
                                     rocblas_int    n,
                                     T              x,
                                     rocblas_stride offsetx,
                                     rocblas_int    incx,
                                     rocblas_stride stridex,
                                     U              y,
                                     rocblas_stride offsety,
                                     rocblas_int    incy,
                                     rocblas_stride stridey,
                                     rocblas_int    batch_count)
{
    // Quick return if possible.
    if(n <= 0 || batch_count <= 0)
        return rocblas_status_success;

    if(!x || !y)
        return rocblas_status_invalid_pointer;

    static constexpr bool using_rocblas_float
        = std::is_same_v<U, rocblas_float*> || std::is_same_v<U, rocblas_float* const*>;

    if(!using_rocblas_float || incx != 1 || incy != 1)
    {
        // In case of negative inc shift pointer to end of data for negative indexing tid*inc
        int64_t shiftx = offsetx - ((incx < 0) ? int64_t(incx) * (n - 1) : 0);
        int64_t shifty = offsety - ((incy < 0) ? int64_t(incy) * (n - 1) : 0);

        int  blocks = (n - 1) / NB + 1;
        dim3 grid(blocks, batch_count);
        dim3 threads(NB);

        hipLaunchKernelGGL(rocblas_copy_kernel,
                           grid,
                           threads,
                           0,
                           handle->get_stream(),
                           n,
                           x,
                           shiftx,
                           incx,
                           stridex,
                           y,
                           shifty,
                           incy,
                           stridey);
    }
    else
    {
        // Kernel function for improving the performance of SCOPY when incx==1 and incy==1

        // In case of negative inc shift pointer to end of data for negative indexing tid*inc
        int64_t shiftx = offsetx - 0;
        int64_t shifty = offsety - 0;

        int         blocks = 1 + ((n - 1) / (NB * 2));
        dim3        grid(blocks, batch_count);
        dim3        threads(NB);
        hipStream_t scopy_stream = handle->get_stream();

        hipLaunchKernelGGL(rocblas_scopy_2_kernel<NB>,
                           grid,
                           threads,
                           0,
                           scopy_stream,
                           n,
                           x,
                           shiftx,
                           stridex,
                           y,
                           shifty,
                           stridey);
    }
    return rocblas_status_success;
}

template <typename T, typename U>
rocblas_status rocblas_copy_check_numerics(const char*    function_name,
                                           rocblas_handle handle,
                                           rocblas_int    n,
                                           T              x,
                                           rocblas_stride offset_x,
                                           rocblas_int    inc_x,
                                           rocblas_stride stride_x,
                                           U              y,
                                           rocblas_stride offset_y,
                                           rocblas_int    inc_y,
                                           rocblas_stride stride_y,
                                           rocblas_int    batch_count,
                                           const int      check_numerics,
                                           bool           is_input)
{
    rocblas_status check_numerics_status
        = rocblas_internal_check_numerics_vector_template(function_name,
                                                          handle,
                                                          n,
                                                          x,
                                                          offset_x,
                                                          inc_x,
                                                          stride_x,
                                                          batch_count,
                                                          check_numerics,
                                                          is_input);
    if(check_numerics_status != rocblas_status_success)
        return check_numerics_status;

    check_numerics_status = rocblas_internal_check_numerics_vector_template(function_name,
                                                                            handle,
                                                                            n,
                                                                            y,
                                                                            offset_y,
                                                                            inc_y,
                                                                            stride_y,
                                                                            batch_count,
                                                                            check_numerics,
                                                                            is_input);

    return check_numerics_status;
}

// Instantiations below will need to be manually updated to match any change in
// template parameters in the files copy*.cpp

#ifdef INSTANTIATE_COPY_TEMPLATE
#error INSTANTIATE_COPY_TEMPLATE already defined
#endif

#define INSTANTIATE_COPY_TEMPLATE(NB_, T_, U_)                                         \
    template rocblas_status rocblas_copy_template<NB_, T_, U_>(rocblas_handle handle,  \
                                                               rocblas_int    n,       \
                                                               T_             x,       \
                                                               rocblas_stride offsetx, \
                                                               rocblas_int    incx,    \
                                                               rocblas_stride stridex, \
                                                               U_             y,       \
                                                               rocblas_stride offsety, \
                                                               rocblas_int    incy,    \
                                                               rocblas_stride stridey, \
                                                               rocblas_int    batch_count);

INSTANTIATE_COPY_TEMPLATE(512, const float*, float*)
INSTANTIATE_COPY_TEMPLATE(256, const float*, float*)

INSTANTIATE_COPY_TEMPLATE(256, const double*, double*)

INSTANTIATE_COPY_TEMPLATE(256, const rocblas_half*, rocblas_half*)

INSTANTIATE_COPY_TEMPLATE(256, const rocblas_float_complex*, rocblas_float_complex*)

INSTANTIATE_COPY_TEMPLATE(256, const rocblas_double_complex*, rocblas_double_complex*)

INSTANTIATE_COPY_TEMPLATE(512, float*, float*)
INSTANTIATE_COPY_TEMPLATE(512, float*, float* const*)
INSTANTIATE_COPY_TEMPLATE(256, float*, float*)
INSTANTIATE_COPY_TEMPLATE(256, float* const*, float* const*)
INSTANTIATE_COPY_TEMPLATE(256, float const* const*, float* const*)

INSTANTIATE_COPY_TEMPLATE(512, double*, double*)
INSTANTIATE_COPY_TEMPLATE(512, double*, double* const*)
INSTANTIATE_COPY_TEMPLATE(256, double*, double*)
INSTANTIATE_COPY_TEMPLATE(256, double* const*, double* const*)
INSTANTIATE_COPY_TEMPLATE(256, double const* const*, double* const*)

INSTANTIATE_COPY_TEMPLATE(256, const rocblas_half* const*, rocblas_half* const*)

INSTANTIATE_COPY_TEMPLATE(512, rocblas_float_complex*, rocblas_float_complex*)
INSTANTIATE_COPY_TEMPLATE(512, rocblas_float_complex*, rocblas_float_complex* const*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_float_complex*, rocblas_float_complex*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_float_complex* const*, rocblas_float_complex* const*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_float_complex const* const*, rocblas_float_complex* const*)

INSTANTIATE_COPY_TEMPLATE(512, rocblas_double_complex*, rocblas_double_complex*)
INSTANTIATE_COPY_TEMPLATE(512, rocblas_double_complex*, rocblas_double_complex* const*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_double_complex*, rocblas_double_complex*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_double_complex* const*, rocblas_double_complex* const*)
INSTANTIATE_COPY_TEMPLATE(256, rocblas_double_complex const* const*, rocblas_double_complex* const*)

#undef INSTANTIATE_COPY_TEMPLATE

#ifdef INSTANTIATE_COPY_CHECK_NUMERICS
#error INSTANTIATE_COPY_CHECK_NUMERICS already defined
#endif

#define INSTANTIATE_COPY_CHECK_NUMERICS(T_, U_)                                                \
    template rocblas_status rocblas_copy_check_numerics<T_, U_>(const char*    function_name,  \
                                                                rocblas_handle handle,         \
                                                                rocblas_int    n,              \
                                                                T_             x,              \
                                                                rocblas_stride offset_x,       \
                                                                rocblas_int    inc_x,          \
                                                                rocblas_stride stride_x,       \
                                                                U_             y,              \
                                                                rocblas_stride offset_y,       \
                                                                rocblas_int    inc_y,          \
                                                                rocblas_stride stride_y,       \
                                                                rocblas_int    batch_count,    \
                                                                const int      check_numerics, \
                                                                bool           is_input);

INSTANTIATE_COPY_CHECK_NUMERICS(const float*, float*)
INSTANTIATE_COPY_CHECK_NUMERICS(const double*, double*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_half*, rocblas_half*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_float_complex*, rocblas_float_complex*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_double_complex*, rocblas_double_complex*)

INSTANTIATE_COPY_CHECK_NUMERICS(const float* const*, float* const*)
INSTANTIATE_COPY_CHECK_NUMERICS(const double* const*, double* const*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_half* const*, rocblas_half* const*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_float_complex* const*, rocblas_float_complex* const*)
INSTANTIATE_COPY_CHECK_NUMERICS(const rocblas_double_complex* const*,
                                rocblas_double_complex* const*)

#undef INSTANTIATE_COPY_CHECK_NUMERICS
