From 63df05d33ab5702383e701867e65b6a5e9b99b5e Mon Sep 17 00:00:00 2001 From: fwenguang <95677191+fwenguang@users.noreply.github.com> Date: Wed, 3 Aug 2022 14:11:12 +0800 Subject: [PATCH] [MLU] add strided_slice kernel (#44460) --- paddle/fluid/memory/memcpy.cc | 2 +- .../fluid/operators/strided_slice_op_mlu.cc | 411 ++++++++++++++++++ paddle/fluid/operators/utils.h | 2 +- 3 files changed, 413 insertions(+), 2 deletions(-) create mode 100644 paddle/fluid/operators/strided_slice_op_mlu.cc diff --git a/paddle/fluid/memory/memcpy.cc b/paddle/fluid/memory/memcpy.cc index 05f46dd396..70583e127d 100644 --- a/paddle/fluid/memory/memcpy.cc +++ b/paddle/fluid/memory/memcpy.cc @@ -1195,7 +1195,7 @@ void Copy(platform::MLUPlace dst_place, dst, src, num, reinterpret_cast(stream)); } else { platform::DeviceContextPool& pool = platform::DeviceContextPool::Instance(); - static_cast(pool.Get(src_place))->Wait(); + static_cast(pool.Get(dst_place))->Wait(); VLOG(4) << "Sync memory::Copy " << num << " Bytes from " << src_place << " to " << dst_place; diff --git a/paddle/fluid/operators/strided_slice_op_mlu.cc b/paddle/fluid/operators/strided_slice_op_mlu.cc new file mode 100644 index 0000000000..95972d8159 --- /dev/null +++ b/paddle/fluid/operators/strided_slice_op_mlu.cc @@ -0,0 +1,411 @@ +/* Copyright (c) 2022 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. */ + +#include "paddle/fluid/operators/mlu/mlu_baseop.h" +#include "paddle/fluid/operators/slice_op.h" +#include "paddle/phi/kernels/funcs/strided_slice.h" + +namespace paddle { +namespace operators { + +static void ProcessStridedSliceParams( + const std::vector& axes, + const DDim& input_dims, + const std::vector& starts, + const std::vector& ends, + const std::vector& strides, + const std::vector& infer_flags, + const std::vector& decrease_axis, + std::vector* starts_indices_vector, + std::vector* ends_indices_vector, + std::vector* strides_indices_vector) { + for (size_t axis = 0; axis < axes.size(); axis++) { + int64_t start = starts[axis]; + int64_t end = ends[axis]; + int64_t stride = strides[axis]; + + int axis_index = axes[axis]; + int64_t dim_size = input_dims[axis_index]; + + bool decrease_axis_affect = false; + if (start == -1 && end == 0 && infer_flags[axis] == -1) { + auto ret = + std::find(decrease_axis.begin(), decrease_axis.end(), axis_index); + if (ret != decrease_axis.end()) { + decrease_axis_affect = true; + } + } + + if (stride < 0) { + if (start < 0) { + start = std::max(start, -dim_size); + } else { + start = std::min(start, dim_size - 1) - dim_size; + } + if (end < 0) { + end = std::max(end, -dim_size - 1); + } else { + end = end - dim_size; + } + } else { + if (start < 0) { + start = std::max(start, -dim_size) + dim_size; + } else { + start = std::min(start, dim_size - 1); + } + if (end < 0) { + end = end + dim_size; + } else { + end = std::min(end, dim_size); + } + } + + if (decrease_axis_affect) { + if (stride < 0) { + end = start - 1; + } else { + end = start + 1; + } + } + + (*starts_indices_vector)[axis_index] = static_cast(start); + (*ends_indices_vector)[axis_index] = static_cast(end); + (*strides_indices_vector)[axis_index] = static_cast(stride); + } +} + +template +class StridedSliceMLUKernel : public framework::OpKernel { + public: + void Compute(const framework::ExecutionContext& ctx) const override { + const Variable* input_var = ctx.InputVar("Input"); + bool is_tensor_array = input_var->IsType(); + PADDLE_ENFORCE_EQ(is_tensor_array, + false, + platform::errors::InvalidArgument( + "Tensor array as input is not supported.")); + int rank = ctx.Input("Input")->dims().size(); + switch (rank) { + case 1: + StridedSliceCompute<1>(ctx); + break; + case 2: + StridedSliceCompute<2>(ctx); + break; + case 3: + StridedSliceCompute<3>(ctx); + break; + case 4: + StridedSliceCompute<4>(ctx); + break; + case 5: + StridedSliceCompute<5>(ctx); + break; + case 6: + StridedSliceCompute<6>(ctx); + break; + case 7: + StridedSliceCompute<7>(ctx); + break; + case 8: + StridedSliceCompute<8>(ctx); + break; + default: + PADDLE_THROW(platform::errors::InvalidArgument( + "The rank of input is supported up to 8.")); + break; + } + } + + private: + template + void StridedSliceCompute(const framework::ExecutionContext& ctx) const { + auto place = ctx.GetPlace(); + + auto in = ctx.Input("Input"); + auto out = ctx.Output("Out"); + auto in_dims = in->dims(); + + // list + auto starts_int = ctx.Attr>("starts"); + auto ends_int = ctx.Attr>("ends"); + auto strides_int = ctx.Attr>("strides"); + + std::vector starts(starts_int.begin(), starts_int.end()); + std::vector ends(ends_int.begin(), ends_int.end()); + std::vector strides(strides_int.begin(), strides_int.end()); + + auto axes = ctx.Attr>("axes"); + auto infer_flags = ctx.Attr>("infer_flags"); + auto decrease_axis = ctx.Attr>("decrease_axis"); + + // vector> + auto list_new_starts_tensor = + ctx.MultiInput("StartsTensorList"); + auto list_new_ends_tensor = + ctx.MultiInput("EndsTensorList"); + auto list_new_strides_tensor = + ctx.MultiInput("StridesTensorList"); + + // Tensor + if (list_new_starts_tensor.size() > 0) { + starts = GetDataFromTensorList(list_new_starts_tensor); + } else if (ctx.HasInput("StartsTensor")) { + auto* starts_tensor = ctx.Input("StartsTensor"); + starts = GetDataFromTensor(starts_tensor); + } + + if (list_new_ends_tensor.size() > 0) { + ends = GetDataFromTensorList(list_new_ends_tensor); + } else if (ctx.HasInput("EndsTensor")) { + auto* ends_tensor = ctx.Input("EndsTensor"); + ends = GetDataFromTensor(ends_tensor); + } + + if (list_new_strides_tensor.size() > 0) { + strides = GetDataFromTensorList(list_new_strides_tensor); + } else if (ctx.HasInput("StridesTensor")) { + auto* strides_tensor = ctx.Input("StridesTensor"); + strides = GetDataFromTensor(strides_tensor); + } + + // out dims calculation + std::vector out_dims_vector(in_dims.size(), -1); + phi::funcs::StridedSliceOutDims(starts, + ends, + strides, + axes, + infer_flags, + in_dims, + decrease_axis, + out_dims_vector.data(), + axes.size(), + false); + framework::DDim out_dims(phi::make_ddim(out_dims_vector)); + + // construct the starts_indices, ends_indices and strides_indices tensor for + // calling StridedSlice op + std::vector starts_indices_vector(D, 0); + std::vector ends_indices_vector(out_dims_vector.begin(), + out_dims_vector.end()); + std::vector strides_indices_vector(D, 1); + + ProcessStridedSliceParams(axes, + in_dims, + starts, + ends, + strides, + infer_flags, + decrease_axis, + &starts_indices_vector, + &ends_indices_vector, + &strides_indices_vector); + + auto out_dims_origin = out_dims; + if (decrease_axis.size() > 0) { + std::vector new_out_shape; + for (size_t i = 0; i < decrease_axis.size(); ++i) { + PADDLE_ENFORCE_EQ( + out_dims[decrease_axis[i]], + 1, + platform::errors::InvalidArgument( + "the size of decrease dimension should be 1, but received %d.", + out_dims[decrease_axis[i]])); + out_dims_origin[decrease_axis[i]] = 0; + } + + for (int i = 0; i < out_dims_origin.size(); ++i) { + if (out_dims_origin[i] != 0) { + new_out_shape.push_back(out_dims_origin[i]); + } + } + if (new_out_shape.size() == 0) { + new_out_shape.push_back(1); + } + out_dims_origin = phi::make_ddim(new_out_shape); + } + + out->Resize(out_dims_origin); + out->mutable_data(place); + + MLUCnnlTensorDesc in_desc(*in); + MLUCnnlTensorDesc out_desc( + out_dims_vector.size(), out_dims_vector.data(), ToCnnlDataType()); + MLUCnnl::StridedSlice(ctx, + starts_indices_vector.data(), + ends_indices_vector.data(), + strides_indices_vector.data(), + in_desc.get(), + GetBasePtr(in), + out_desc.get(), + GetBasePtr(out)); + } +}; + +template +class StridedSliceGradMLUKernel : public framework::OpKernel { + public: + void Compute(const framework::ExecutionContext& ctx) const override { + const Variable* input_var = ctx.InputVar("Input"); + bool is_tensor_array = input_var->IsType(); + PADDLE_ENFORCE_EQ(is_tensor_array, + false, + platform::errors::InvalidArgument( + "Tensor array as input is not supported.")); + int rank = ctx.Input("Input")->dims().size(); + + switch (rank) { + case 1: + StridedSliceGradCompute<1>(ctx); + break; + case 2: + StridedSliceGradCompute<2>(ctx); + break; + case 3: + StridedSliceGradCompute<3>(ctx); + break; + case 4: + StridedSliceGradCompute<4>(ctx); + break; + case 5: + StridedSliceGradCompute<5>(ctx); + break; + case 6: + StridedSliceGradCompute<6>(ctx); + break; + case 7: + StridedSliceGradCompute<7>(ctx); + break; + case 8: + StridedSliceGradCompute<8>(ctx); + break; + default: + PADDLE_THROW(platform::errors::InvalidArgument( + "The rank of input is supported up to 8.")); + break; + } + } + + private: + template + void StridedSliceGradCompute(const framework::ExecutionContext& ctx) const { + auto place = ctx.GetPlace(); + + auto* input = ctx.Input("Input"); + auto input_dims = input->dims(); + auto* dout = ctx.Input(framework::GradVarName("Out")); + auto* dx = ctx.Output(framework::GradVarName("Input")); + dx->mutable_data(input_dims, place); + + auto starts_int = ctx.Attr>("starts"); + auto ends_int = ctx.Attr>("ends"); + auto strides_int = ctx.Attr>("strides"); + + std::vector starts(starts_int.begin(), starts_int.end()); + std::vector ends(ends_int.begin(), ends_int.end()); + std::vector strides(strides_int.begin(), strides_int.end()); + + auto axes = ctx.Attr>("axes"); + auto infer_flags = ctx.Attr>("infer_flags"); + auto decrease_axis = ctx.Attr>("decrease_axis"); + + auto list_new_ends_tensor = + ctx.MultiInput("EndsTensorList"); + auto list_new_starts_tensor = + ctx.MultiInput("StartsTensorList"); + auto list_new_strides_tensor = + ctx.MultiInput("StridesTensorList"); + + if (list_new_starts_tensor.size() > 0) { + starts = GetDataFromTensorList(list_new_starts_tensor); + } else if (ctx.HasInput("StartsTensor")) { + auto* starts_tensor = ctx.Input("StartsTensor"); + starts = GetDataFromTensor(starts_tensor); + } + + if (list_new_ends_tensor.size() > 0) { + ends = GetDataFromTensorList(list_new_ends_tensor); + } else if (ctx.HasInput("EndsTensor")) { + auto* ends_tensor = ctx.Input("EndsTensor"); + ends = GetDataFromTensor(ends_tensor); + } + + if (list_new_strides_tensor.size() > 0) { + strides = GetDataFromTensorList(list_new_strides_tensor); + } else if (ctx.HasInput("StridesTensor")) { + auto* strides_tensor = ctx.Input("StridesTensor"); + strides = GetDataFromTensor(strides_tensor); + } + + std::vector out_dims_vector(input_dims.size(), -1); + phi::funcs::StridedSliceOutDims(starts, + ends, + strides, + axes, + infer_flags, + input_dims, + decrease_axis, + out_dims_vector.data(), + axes.size(), + false); + + std::vector starts_indices_vector(D, 0); + std::vector ends_indices_vector(out_dims_vector.begin(), + out_dims_vector.end()); + std::vector strides_indices_vector(D, 1); + + ProcessStridedSliceParams(axes, + input_dims, + starts, + ends, + strides, + infer_flags, + decrease_axis, + &starts_indices_vector, + &ends_indices_vector, + &strides_indices_vector); + + MLUCnnlTensorDesc dout_desc( + out_dims_vector.size(), out_dims_vector.data(), ToCnnlDataType()); + MLUCnnlTensorDesc dx_desc(*input); + MLUCnnl::StridedSliceGrad(ctx, + starts_indices_vector.data(), + ends_indices_vector.data(), + strides_indices_vector.data(), + dout_desc.get(), + GetBasePtr(dout), + dx_desc.get(), + GetBasePtr(dx)); + } +}; + +} // namespace operators +} // namespace paddle + +namespace ops = paddle::operators; +namespace plat = paddle::platform; + +REGISTER_OP_MLU_KERNEL(strided_slice, + ops::StridedSliceMLUKernel, + ops::StridedSliceMLUKernel, + ops::StridedSliceMLUKernel, + ops::StridedSliceMLUKernel, + ops::StridedSliceMLUKernel); + +REGISTER_OP_MLU_KERNEL(strided_slice_grad, + ops::StridedSliceGradMLUKernel, + ops::StridedSliceGradMLUKernel, + ops::StridedSliceGradMLUKernel, + ops::StridedSliceGradMLUKernel, + ops::StridedSliceGradMLUKernel); diff --git a/paddle/fluid/operators/utils.h b/paddle/fluid/operators/utils.h index 51a5e2d52d..7315f3a287 100644 --- a/paddle/fluid/operators/utils.h +++ b/paddle/fluid/operators/utils.h @@ -122,7 +122,7 @@ inline T GetValue(const framework::Tensor* x) { if (!platform::is_cpu_place(x->place())) { framework::Tensor cpu_x; framework::TensorCopy(*x, platform::CPUPlace(), &cpu_x); -#ifdef PADDLE_WITH_ASCEND_CL +#if defined(PADDLE_WITH_ASCEND_CL) || defined(PADDLE_WITH_MLU) platform::DeviceContextPool& pool = platform::DeviceContextPool::Instance(); const platform::DeviceContext* dev_ctx = pool.Get(x->place()); dev_ctx->Wait(); -- GitLab