transpose_op.h 3.4 KB
Newer Older
1
/* Copyright (c) 2016 PaddlePaddle Authors. All Rights Reserved.
X
xzl 已提交
2

L
Luo Tao 已提交
3 4 5
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
X
xzl 已提交
6

L
Luo Tao 已提交
7
    http://www.apache.org/licenses/LICENSE-2.0
X
xzl 已提交
8

L
Luo Tao 已提交
9 10 11 12 13
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. */
X
xzl 已提交
14 15 16

#pragma once

17
#include <vector>
Y
Yi Wang 已提交
18 19
#include "paddle/fluid/framework/op_registry.h"
#include "paddle/fluid/operators/math/math_function.h"
X
xzl 已提交
20 21 22 23

namespace paddle {
namespace operators {

24 25
enum { kTransposeMKLDNNFP32 = 1, kTransposeMKLDNNINT8 = 2 };

Q
QI JUN 已提交
26 27
template <typename DeviceContext, typename T>
inline void TransCompute(const int dim, const DeviceContext& dev_ctx,
28 29 30 31
                         const framework::Tensor& in, framework::Tensor* out,
                         const std::vector<int>& axis) {
  switch (dim) {
    case 1:
Q
QI JUN 已提交
32
      math::Transpose<DeviceContext, T, 1> trans1;
33 34 35
      trans1(dev_ctx, in, out, axis);
      break;
    case 2:
Q
QI JUN 已提交
36
      math::Transpose<DeviceContext, T, 2> trans2;
37 38 39
      trans2(dev_ctx, in, out, axis);
      break;
    case 3:
Q
QI JUN 已提交
40
      math::Transpose<DeviceContext, T, 3> trans3;
41 42 43
      trans3(dev_ctx, in, out, axis);
      break;
    case 4:
Q
QI JUN 已提交
44
      math::Transpose<DeviceContext, T, 4> trans4;
45 46 47
      trans4(dev_ctx, in, out, axis);
      break;
    case 5:
Q
QI JUN 已提交
48
      math::Transpose<DeviceContext, T, 5> trans5;
49 50 51
      trans5(dev_ctx, in, out, axis);
      break;
    case 6:
Q
QI JUN 已提交
52
      math::Transpose<DeviceContext, T, 6> trans6;
53 54 55 56
      trans6(dev_ctx, in, out, axis);
      break;
    default:
      PADDLE_THROW("Tensors with rank at most 6 are supported");
X
xzl 已提交
57 58 59
  }
}

Q
QI JUN 已提交
60
template <typename DeviceContext, typename T>
Y
Yu Yang 已提交
61
class TransposeKernel : public framework::OpKernel<T> {
X
xzl 已提交
62 63
 public:
  void Compute(const framework::ExecutionContext& context) const override {
X
xzl 已提交
64 65 66
    auto* x = context.Input<framework::Tensor>("X");
    auto* out = context.Output<framework::Tensor>("Out");
    out->mutable_data<T>(context.GetPlace());
67 68 69
    if (out->numel() == 0) {
      return;
    }
70
    std::vector<int> axis = context.Attr<std::vector<int>>("axis");
X
xzl 已提交
71
    int ndims = axis.size();
Q
QI JUN 已提交
72 73
    auto& dev_ctx = context.template device_context<DeviceContext>();
    TransCompute<DeviceContext, T>(ndims, dev_ctx, *x, out, axis);
X
xzl 已提交
74 75 76
  }
};

Q
QI JUN 已提交
77
template <typename DeviceContext, typename T>
Y
Yu Yang 已提交
78
class TransposeGradKernel : public framework::OpKernel<T> {
X
xzl 已提交
79 80
 public:
  void Compute(const framework::ExecutionContext& context) const override {
X
xzl 已提交
81 82 83 84
    auto* out_grad =
        context.Input<framework::Tensor>(framework::GradVarName("Out"));
    auto* x_grad =
        context.Output<framework::Tensor>(framework::GradVarName("X"));
85
    if (!x_grad) return;
86

87
    x_grad->mutable_data<T>(context.GetPlace());
88 89 90 91
    if (x_grad->numel() == 0) {
      return;
    }

92 93
    std::vector<int> axis = context.Attr<std::vector<int>>("axis");
    std::vector<int> reversed_axis(axis);
94

95 96
    for (size_t i = 0; i < axis.size(); i++) {
      reversed_axis[axis[i]] = i;
X
xzl 已提交
97
    }
98 99

    int ndims = axis.size();
Q
QI JUN 已提交
100 101 102
    auto& dev_ctx = context.template device_context<DeviceContext>();
    TransCompute<DeviceContext, T>(ndims, dev_ctx, *out_grad, x_grad,
                                   reversed_axis);
X
xzl 已提交
103 104 105 106 107
  }
};

}  // namespace operators
}  // namespace paddle