unpooling.cc 4.3 KB
Newer Older
1
/* Copyright (c) 2016 PaddlePaddle Authors. All Rights Reserved.
S
sweetsky0901 已提交
2 3 4 5 6 7 8 9 10 11 12 13 14

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. */

Y
Yi Wang 已提交
15
#include "paddle/fluid/operators/math/unpooling.h"
16

S
sweetsky0901 已提交
17 18 19
namespace paddle {
namespace operators {
namespace math {
S
sweetsky0901 已提交
20
template <typename T>
Q
QI JUN 已提交
21
class Unpool2dMaxFunctor<platform::CPUDeviceContext, T> {
S
sweetsky0901 已提交
22
 public:
Q
QI JUN 已提交
23
  void operator()(const platform::CPUDeviceContext& context,
S
sweetsky0901 已提交
24 25
                  const framework::Tensor& input,
                  const framework::Tensor& indices, framework::Tensor* output) {
S
sweetsky0901 已提交
26 27 28 29 30 31 32 33 34
    const int batch_size = input.dims()[0];
    const int input_height = input.dims()[2];
    const int input_width = input.dims()[3];
    const int output_channels = output->dims()[1];
    const int output_height = output->dims()[2];
    const int output_width = output->dims()[3];
    int input_feasize = input_height * input_width;
    int output_feasize = output_height * output_width;
    const T* input_data = input.data<T>();
S
sweetsky0901 已提交
35
    const int* indices_data = indices.data<int>();
S
sweetsky0901 已提交
36 37 38 39
    T* output_data = output->mutable_data<T>(context.GetPlace());
    for (int b = 0; b < batch_size; ++b) {
      for (int c = 0; c < output_channels; ++c) {
        for (int i = 0; i < input_feasize; ++i) {
S
sweetsky0901 已提交
40
          int index = indices_data[i];
41 42 43 44 45 46 47
          PADDLE_ENFORCE_LT(
              index, output_feasize,
              platform::errors::InvalidArgument(
                  "index should less than output tensor height * output tensor "
                  "width. Expected %ld < %ld, but got "
                  "%ld >= %ld. Please check input value.",
                  index, output_feasize, index, output_feasize));
S
sweetsky0901 已提交
48 49 50 51 52 53 54 55 56
          output_data[index] = input_data[i];
        }
        input_data += input_feasize;
        indices_data += input_feasize;
        output_data += output_feasize;
      }
    }
  }
};
S
sweetsky0901 已提交
57
template <class T>
Q
QI JUN 已提交
58
class Unpool2dMaxGradFunctor<platform::CPUDeviceContext, T> {
S
sweetsky0901 已提交
59
 public:
Q
QI JUN 已提交
60
  void operator()(const platform::CPUDeviceContext& context,
S
sweetsky0901 已提交
61 62 63 64 65
                  const framework::Tensor& input,
                  const framework::Tensor& indices,
                  const framework::Tensor& output,
                  const framework::Tensor& output_grad,
                  framework::Tensor* input_grad) {
S
sweetsky0901 已提交
66 67 68
    const int batch_size = input.dims()[0];
    const int input_height = input.dims()[2];
    const int input_width = input.dims()[3];
S
sweetsky0901 已提交
69 70 71
    const int output_channels = output.dims()[1];
    const int output_height = output.dims()[2];
    const int output_width = output.dims()[3];
S
sweetsky0901 已提交
72 73
    int input_feasize = input_height * input_width;
    int output_feasize = output_height * output_width;
S
sweetsky0901 已提交
74
    const int* indices_data = indices.data<int>();
S
sweetsky0901 已提交
75 76 77 78 79
    const T* output_grad_data = output_grad.data<T>();
    T* input_grad_data = input_grad->mutable_data<T>(context.GetPlace());

    for (int b = 0; b < batch_size; ++b) {
      for (int c = 0; c < output_channels; ++c) {
S
sweetsky0901 已提交
80
        for (int i = 0; i < input_feasize; ++i) {
S
sweetsky0901 已提交
81
          int index = indices_data[i];
82 83 84 85 86 87 88
          PADDLE_ENFORCE_LT(
              index, output_feasize,
              platform::errors::InvalidArgument(
                  "index should less than output tensor height * output tensor "
                  "width. Expected %ld < %ld, but got "
                  "%ld >= %ld. Please check input value.",
                  index, output_feasize, index, output_feasize));
S
sweetsky0901 已提交
89 90 91 92 93 94 95 96 97
          input_grad_data[i] = output_grad_data[index];
        }
        input_grad_data += input_feasize;
        indices_data += input_feasize;
        output_grad_data += output_feasize;
      }
    }
  }
};
Q
QI JUN 已提交
98 99 100 101
template class Unpool2dMaxGradFunctor<platform::CPUDeviceContext, float>;
template class Unpool2dMaxGradFunctor<platform::CPUDeviceContext, double>;
template class Unpool2dMaxFunctor<platform::CPUDeviceContext, float>;
template class Unpool2dMaxFunctor<platform::CPUDeviceContext, double>;
S
sweetsky0901 已提交
102 103 104
}  // namespace math
}  // namespace operators
}  // namespace paddle