concat.cc 3.8 KB
Newer Older
C
chengduoZH 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
/* Copyright (c) 2018 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/math/concat.h"
16
#include <vector>
C
chengduoZH 已提交
17 18 19 20 21 22

namespace paddle {
namespace operators {
namespace math {

/*
C
chengduoZH 已提交
23
 * All tensors' dimension should be the same and the values of
24
 * each dimension must be the same, except the axis dimension.
C
chengduoZH 已提交
25 26 27 28 29
 */
template <typename T>
class ConcatFunctor<platform::CPUDeviceContext, T> {
 public:
  void operator()(const platform::CPUDeviceContext& context,
C
chengduoZH 已提交
30
                  const std::vector<framework::Tensor>& input, int axis,
C
chengduoZH 已提交
31
                  framework::Tensor* output) {
C
chengduoZH 已提交
32
    // TODO(zcd): Add input data validity checking
C
chengduoZH 已提交
33 34 35 36 37 38 39 40 41
    int num = input.size();

    int rows = 1;
    auto dim_0 = input[0].dims();
    for (int i = 0; i < axis; ++i) {
      rows *= dim_0[i];
    }
    int out_rows = rows, out_cols = 0;

C
chengduoZH 已提交
42
    std::vector<int64_t> input_cols(input.size());
C
chengduoZH 已提交
43 44 45
    for (int i = 0; i < num; ++i) {
      int t_cols = input[i].numel() / rows;
      out_cols += t_cols;
C
chengduoZH 已提交
46
      input_cols[i] = t_cols;
C
chengduoZH 已提交
47
    }
X
xuwei06 已提交
48
    auto cpu_place = boost::get<platform::CPUPlace>(context.GetPlace());
C
chengduoZH 已提交
49

C
chengduoZH 已提交
50
    // computation
L
luotao1 已提交
51 52 53 54 55 56 57 58
    auto output_data = output->data<T>();
    int col_idx = 0;
    for (int j = 0; j < num; ++j) {
      int col_len = input_cols[j];
      auto input_data = input[j].data<T>();
      for (int k = 0; k < out_rows; ++k) {
        memory::Copy(cpu_place, output_data + k * out_cols + col_idx, cpu_place,
                     input_data + k * col_len, sizeof(T) * col_len);
C
chengduoZH 已提交
59
      }
L
luotao1 已提交
60
      col_idx += col_len;
C
chengduoZH 已提交
61
    }
C
chengduoZH 已提交
62 63 64
  }
};

C
chengduoZH 已提交
65 66
/*
 * All tensors' dimension should be the same and the values of
67
 * each dimension must be the same, except the axis dimension.
C
chengduoZH 已提交
68
 */
C
chengduoZH 已提交
69 70 71 72
template <typename T>
class ConcatGradFunctor<platform::CPUDeviceContext, T> {
 public:
  void operator()(const platform::CPUDeviceContext& context,
Q
qiaolongfei 已提交
73
                  const framework::Tensor& input,
C
chengduoZH 已提交
74
                  const std::vector<const framework::Tensor*>& ref_inputs,
Q
qiaolongfei 已提交
75
                  const int axis, std::vector<framework::Tensor*>* outputs) {
C
chengduoZH 已提交
76
    // TODO(zcd): Add input data validity checking
Q
qiaolongfei 已提交
77
    size_t num = outputs->size();
C
chengduoZH 已提交
78

C
chengduoZH 已提交
79
    int input_rows = 1;
Q
qiaolongfei 已提交
80
    auto dim_0 = ref_inputs[0]->dims();
C
chengduoZH 已提交
81 82 83
    for (int i = 0; i < axis; ++i) {
      input_rows *= dim_0[i];
    }
Q
qiaolongfei 已提交
84

C
chengduoZH 已提交
85 86
    int input_cols = 0;

87
    std::vector<int64_t> output_cols(outputs->size());
Q
qiaolongfei 已提交
88 89
    for (size_t i = 0; i < num; ++i) {
      int t_cols = ref_inputs[i]->numel() / input_rows;
C
chengduoZH 已提交
90 91 92
      input_cols += t_cols;
      output_cols[i] = t_cols;
    }
X
xuwei06 已提交
93
    auto cpu_place = boost::get<platform::CPUPlace>(context.GetPlace());
C
chengduoZH 已提交
94 95

    // computation
C
chengduoZH 已提交
96
    for (int k = 0; k < input_rows; ++k) {
C
chengduoZH 已提交
97 98
      const T* src_ptr = input.data<T>() + k * input_cols;
      int col_idx = 0;
C
chengduoZH 已提交
99
      for (size_t j = 0; j < num; ++j) {
C
chengduoZH 已提交
100
        int col_len = output_cols[j];
Q
qiaolongfei 已提交
101
        auto* out_tensor = outputs->at(j);
Q
qiaolongfei 已提交
102 103 104 105 106
        if (out_tensor != nullptr) {
          T* dst_ptr = out_tensor->data<T>() + k * col_len;
          memory::Copy(cpu_place, dst_ptr, cpu_place, src_ptr + col_idx,
                       sizeof(T) * col_len);
        }
C
chengduoZH 已提交
107 108 109
        col_idx += col_len;
      }
    }
C
chengduoZH 已提交
110 111
  }
};
C
chengduoZH 已提交
112 113 114
#define DEFINE_FUNCTOR(type)                                      \
  template class ConcatFunctor<platform::CPUDeviceContext, type>; \
  template class ConcatGradFunctor<platform::CPUDeviceContext, type>;
C
chengduoZH 已提交
115

C
chengduoZH 已提交
116
FOR_ALL_TYPES(DEFINE_FUNCTOR);
C
chengduoZH 已提交
117

C
chengduoZH 已提交
118 119 120
}  // namespace math
}  // namespace operators
}  // namespace paddle