SimpleVector.h

Go to the documentation of this file.

#ifndef HeterogeneousCore_CUDAUtilities_interface_SimpleVector_h
#define HeterogeneousCore_CUDAUtilities_interface_SimpleVector_h
 
//  author: Felice Pantaleo, CERN, 2018
 
#include <type_traits>
#include <utility>
 
#include "HeterogeneousCore/CUDAUtilities/interface/cudaCompat.h"
 
namespace cms {
  namespace cuda {
 
    template <class T>
    struct SimpleVector {
      constexpr SimpleVector() = default;
 
      // ownership of m_data stays within the caller
      constexpr void construct(int capacity, T *data) {
        m_size = 0;
        m_capacity = capacity;
        m_data = data;
      }
 
      inline constexpr int push_back_unsafe(const T &element) {
        auto previousSize = m_size;
        m_size++;
        if (previousSize < m_capacity) {
          m_data[previousSize] = element;
          return previousSize;
        } else {
          --m_size;
          return -1;
        }
      }
 
      template <class... Ts>
      constexpr int emplace_back_unsafe(Ts &&...args) {
        auto previousSize = m_size;
        m_size++;
        if (previousSize < m_capacity) {
          (new (&m_data[previousSize]) T(std::forward<Ts>(args)...));
          return previousSize;
        } else {
          --m_size;
          return -1;
        }
      }
 
      __device__ inline T &back() { return m_data[m_size - 1]; }
 
      __device__ inline const T &back() const {
        if (m_size > 0) {
          return m_data[m_size - 1];
        } else
          return T();  //undefined behaviour
      }
 
      // thread-safe version of the vector, when used in a CUDA kernel
      __device__ int push_back(const T &element) {
        auto previousSize = atomicAdd(&m_size, 1);
        if (previousSize < m_capacity) {
          m_data[previousSize] = element;
          return previousSize;
        } else {
          atomicSub(&m_size, 1);
          return -1;
        }
      }
 
      template <class... Ts>
      __device__ int emplace_back(Ts &&...args) {
        auto previousSize = atomicAdd(&m_size, 1);
        if (previousSize < m_capacity) {
          (new (&m_data[previousSize]) T(std::forward<Ts>(args)...));
          return previousSize;
        } else {
          atomicSub(&m_size, 1);
          return -1;
        }
      }
 
      // thread safe version of resize
      __device__ int extend(int size = 1) {
        auto previousSize = atomicAdd(&m_size, size);
        if (previousSize < m_capacity) {
          return previousSize;
        } else {
          atomicSub(&m_size, size);
          return -1;
        }
      }
 
      __device__ int shrink(int size = 1) {
        auto previousSize = atomicSub(&m_size, size);
        if (previousSize >= size) {
          return previousSize - size;
        } else {
          atomicAdd(&m_size, size);
          return -1;
        }
      }
 
      inline constexpr bool empty() const { return m_size <= 0; }
      inline constexpr bool full() const { return m_size >= m_capacity; }
      inline constexpr T &operator[](int i) { return m_data[i]; }
      inline constexpr const T &operator[](int i) const { return m_data[i]; }
      inline constexpr void reset() { m_size = 0; }
      inline constexpr int size() const { return m_size; }
      inline constexpr int capacity() const { return m_capacity; }
      inline constexpr T const *data() const { return m_data; }
      inline constexpr void resize(int size) { m_size = size; }
      inline constexpr void set_data(T *data) { m_data = data; }
 
    private:
      int m_size;
      int m_capacity;
 
      T *m_data;
    };
 
    // ownership of m_data stays within the caller
    template <class T>
    SimpleVector<T> make_SimpleVector(int capacity, T *data) {
      SimpleVector<T> ret;
      ret.construct(capacity, data);
      return ret;
    }
 
    // ownership of m_data stays within the caller
    template <class T>
    SimpleVector<T> *make_SimpleVector(SimpleVector<T> *mem, int capacity, T *data) {
      auto ret = new (mem) SimpleVector<T>();
      ret->construct(capacity, data);
      return ret;
    }
 
  }  // namespace cuda
}  // namespace cms
 
#endif  // HeterogeneousCore_CUDAUtilities_interface_SimpleVector_h