strided_memcpy.h

/* Copyright (c) 2023 PaddlePaddle Authors. All Rights Reserved.
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
    http://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License. */

#pragma once
#include <vector>

#include "paddle/phi/kernels/funcs/detail/strided_memcpy.h"

#include "paddle/fluid/platform/device_context.h"
#include "paddle/phi/core/dense_tensor.h"

namespace phi {
namespace funcs {

// Strided memory copy from src to dst.
//
// The src and dst should be both on dev_ctx.GetPlace(), otherwise, there will
// be a segment fault.
//
// The stride of an array (also referred to as increment, pitch or step size) is
// the number of locations in memory between beginnings of successive array
// elements
//
// For example, for tensor like [1, 3, 300, 300]. If there is no padding, the
// stride is [270000, 90000, 300, 1].
//
// NOTE: When use GPU, the memcpy is async. To sync memcpy, please invoke
// `dev_ctx.Wait()`.
template <typename T>
inline void StridedMemcpy(const phi::DeviceContext& dev_ctx,
                          const T* src,
                          const phi::DDim& src_stride,
                          const phi::DDim& dst_dim,
                          const phi::DDim& dst_stride,
                          T* dst) {
  detail::StridedCopyDimVisitor<T> func(
      dev_ctx, src, src_stride, dst_stride, dst);
  dst_dim.apply_visitor(func);
}

// Strided numel memory copy from src to dst by the specified axis
//
// For example, for a tensor dims [4, 20, 100], the strieded numel is
// [8000, 2000, 100]
//
// NOTE: The src and dst tensor should have the same elements
// except the specified axis.
template <typename T>
inline void StridedNumelCopyWithAxis(const phi::DeviceContext& ctx,
                                     int64_t axis,
                                     T* dst,
                                     const phi::DDim& dst_stride_numel,
                                     const T* src,
                                     const phi::DDim& src_stride_numel,
                                     int64_t size) {
  int64_t before = dst_stride_numel[0] / dst_stride_numel[axis];
  int64_t src_after = src_stride_numel[axis];
  int64_t dst_after = dst_stride_numel[axis];
  auto place = ctx.GetPlace();

  PADDLE_ENFORCE_EQ(src_stride_numel.size(),
                    dst_stride_numel.size(),
                    phi::errors::InvalidArgument(
                        "Source and destination tensor should have the same "
                        "dimension size, but source tensor dimension size is "
                        "%u, destination tensor size is %u.",
                        src_stride_numel.size(),
                        dst_stride_numel.size()));

  for (int64_t i = 0; i < axis; ++i) {
    if (i < axis) {
      PADDLE_ENFORCE_EQ(
          src_stride_numel[i] / src_stride_numel[axis],
          dst_stride_numel[i] / dst_stride_numel[axis],
          phi::errors::InvalidArgument(
              "Source and destination tensor should have the same number of "
              "elements except the specified axis, but the source elements "
              "number is %d, destination elements number is %d.",
              src_stride_numel[i] / src_stride_numel[axis],
              dst_stride_numel[i] / dst_stride_numel[axis]));
    } else if (i == axis) {
      continue;
    } else {
      PADDLE_ENFORCE_EQ(
          src_stride_numel[i],
          dst_stride_numel[i],
          phi::errors::InvalidArgument(
              "Source and destination tensor should have the same number of "
              "elements except the specified axis, but the source elements "
              "number is %d, destination elements number is %d.",
              src_stride_numel[i],
              dst_stride_numel[i]));
    }
  }

  for (int64_t i = 0; i < before; ++i) {
    if (place.GetType() == phi::AllocationType::CPU) {
      auto& cpu_place = place;
      paddle::memory::Copy(cpu_place,
                           dst + i * dst_after,
                           cpu_place,
                           src + i * src_after,
                           sizeof(T) * size);
    } else {
#if defined(PADDLE_WITH_CUDA) || defined(PADDLE_WITH_HIP)
      auto& gpu_place = place;
      auto& cuda_ctx = reinterpret_cast<const phi::GPUContext&>(ctx);
      paddle::memory::Copy(gpu_place,
                           dst + i * dst_after,
                           gpu_place,
                           src + i * src_after,
                           sizeof(T) * size,
                           cuda_ctx.stream());
#elif defined(PADDLE_WITH_ASCEND_CL)
      auto& npu_place = place;
      auto& npu_ctx = reinterpret_cast<const platform::NPUDeviceContext&>(ctx);
      paddle::memory::Copy(npu_place,
                           dst + i * dst_after,
                           npu_place,
                           src + i * src_after,
                           sizeof(T) * size,
                           npu_ctx.stream());
#elif defined(PADDLE_WITH_MLU)
      auto& mlu_place = place;
      auto& mlu_ctx = reinterpret_cast<const platform::MLUDeviceContext&>(ctx);
      paddle::memory::Copy(mlu_place,
                           dst + i * dst_after,
                           mlu_place,
                           src + i * src_after,
                           sizeof(T) * size,
                           mlu_ctx.stream());
#else
      PADDLE_THROW(
          phi::errors::PreconditionNotMet("Paddle is not compiled with GPU."));
#endif
    }
  }
}

template <typename T>
inline void StridedMemcpyWithAxis0(
    const phi::DeviceContext& dev_ctx,
    const phi::DenseTensor& input,
    const std::vector<const phi::DenseTensor*>& shape_refer,
    std::vector<phi::DenseTensor*>* outputs) {
  const phi::DDim in_stride = stride_numel(input.dims());
  const int axis = 0;
  size_t input_offset = 0;

  for (size_t i = 0; i < outputs->size(); ++i) {
    auto out_stride = stride_numel(shape_refer[i]->dims());
    auto out = outputs->at(i);
    if (out != nullptr && out->initialized() && out->numel() > 0) {
      StridedNumelCopyWithAxis<T>(dev_ctx,
                                  axis,
                                  out->data<T>(),
                                  out_stride,
                                  input.data<T>() + input_offset,
                                  in_stride,
                                  out_stride[axis]);
    }
    input_offset += out_stride[axis];
  }
}

}  // namespace funcs
}  // namespace phi