transform.h 2.7 KB
Newer Older
Y
Yu Yang 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
/* Copyright (c) 2016 PaddlePaddle Authors. All Rights Reserve.

   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 "paddle/platform/enforce.h"
#include "paddle/platform/hostdevice.h"
#include "paddle/platform/place.h"

#include <algorithm>
Y
Yu Yang 已提交
22
#include <type_traits>
Y
Yu Yang 已提交
23
#ifdef __NVCC__
Y
Yu Yang 已提交
24
#include <thrust/device_ptr.h>
Y
Yu Yang 已提交
25 26 27 28 29 30
#include <thrust/transform.h>
#endif

namespace paddle {
namespace platform {

Y
Yu Yang 已提交
31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58
#ifdef __NVCC__
template <typename T, bool is_ptr>
struct DevicePtrCast;

template <typename T>
struct DevicePtrCast<T, true> {
  using ELEM = typename std::remove_pointer<T>::type;
  using RTYPE = thrust::device_ptr<ELEM>;

  inline thrust::device_ptr<ELEM> operator()(ELEM* ele) const {
    return thrust::device_pointer_cast(ele);
  }
};

template <typename T>
struct DevicePtrCast<T, false> {
  using RTYPE = T;
  inline RTYPE operator()(RTYPE it) const { return it; }
};

template <typename T>
auto DevCast(T t) ->
    typename DevicePtrCast<T, std::is_pointer<T>::value>::RTYPE {
  DevicePtrCast<T, std::is_pointer<T>::value> cast;
  return cast(t);
}
#endif

Y
Yu Yang 已提交
59 60 61 62 63 64 65 66 67
// Transform on host or device. It provides the same API in std library.
template <typename Place, typename InputIter, typename OutputIter,
          typename UnaryOperation>
void Transform(Place place, InputIter first, InputIter last, OutputIter result,
               UnaryOperation op) {
  if (is_cpu_place(place)) {
    std::transform(first, last, result, op);
  } else {
#ifdef __NVCC__
Y
Yu Yang 已提交
68
    thrust::transform(DevCast(first), DevCast(last), DevCast(result), op);
Y
Yu Yang 已提交
69 70 71 72 73 74 75 76 77 78 79 80 81 82
#else
    PADDLE_THROW("Do not invoke `Transform<GPUPlace>` in .cc file");
#endif
  }
}

template <typename Place, typename InputIter1, typename InputIter2,
          typename OutputIter, typename BinaryOperation>
void Transform(Place place, InputIter1 first1, InputIter1 last1,
               InputIter2 first2, OutputIter result, BinaryOperation op) {
  if (is_cpu_place(place)) {
    std::transform(first1, last1, first2, result, op);
  } else {
#ifdef __NVCC__
Y
Yu Yang 已提交
83 84
    thrust::transform(DevCast(first1), DevCast(last1), DevCast(first2),
                      DevCast(result), op);
Y
Yu Yang 已提交
85 86 87 88 89 90 91 92
#else
    PADDLE_THROW("Do not invoke `Transform<GPUPlace>` in .cc file");
#endif
  }
};

}  // namespace platform
}  // namespace paddle