container.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <initializer_list>
#include <vector>
 
#include <iganet/core/options.hpp>
 
namespace iganet::utils {
 
template <std::size_t N, typename T>
inline std::array<T, N> to_array(std::vector<T> &&vector) {
  if (vector.size() != N)
    throw std::invalid_argument("Cannot convert std::vector to std::array: size mismatch");
  
  std::array<T, N> array;
  std::move(vector.begin(), vector.end(), array.begin());
  return array;
}
 
template <typename T, std::size_t N>
inline std::vector<T> to_vector(std::array<T, N> &&array) {
  std::vector<T> vector(N);
  std::move(array.begin(), array.end(), vector.begin());
  return vector;
}
 
template <typename... Args> inline auto to_array(Args &&...args) {
  return std::array<std::common_type_t<Args...>, sizeof...(Args)>{
      std::move(args)...};
}
 
template <typename... Args> inline auto to_vector(Args &&...args) {
  return std::vector<std::common_type_t<Args...>>{std::move(args)...};
}
 
template <typename T, std::size_t N>
inline auto
to_tensor(const std::array<T, N> &array,
          torch::IntArrayRef sizes = torch::IntArrayRef{-1},
          const iganet::Options<T> &options = iganet::Options<T>{}) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(const_cast<T *>(std::data(array)),
                            (sizes == torch::IntArrayRef{-1}) ? array.size()
                                                              : sizes,
                            options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(array)),
                            (sizes == torch::IntArrayRef{-1}) ? array.size()
                                                              : sizes,
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T, std::size_t N>
inline auto to_tensor(const std::array<T, N> &array,
                      const iganet::Options<T> &options) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(const_cast<T *>(std::data(array)), array.size(),
                            options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(array)), array.size(),
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T>
inline auto
to_tensor(std::initializer_list<T> list,
          torch::IntArrayRef sizes = torch::IntArrayRef{-1},
          const iganet::Options<T> &options = iganet::Options<T>{}) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(
               const_cast<T *>(std::data(list)),
               (sizes == torch::IntArrayRef{-1}) ? list.size() : sizes, options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(list)),
                            (sizes == torch::IntArrayRef{-1}) ? list.size()
                                                              : sizes,
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T>
inline auto to_tensor(std::initializer_list<T> &list,
                      const iganet::Options<T> &options) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(const_cast<T *>(std::data(list)), list.size(),
                            options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(list)), list.size(),
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T>
inline auto
to_tensor(const std::vector<T> &vector,
          torch::IntArrayRef sizes = torch::IntArrayRef{-1},
          const iganet::Options<T> &options = iganet::Options<T>{}) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(const_cast<T *>(std::data(vector)),
                            (sizes == torch::IntArrayRef{-1}) ? vector.size()
                                                              : sizes,
                            options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(vector)),
                            (sizes == torch::IntArrayRef{-1}) ? vector.size()
                                                              : sizes,
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T>
inline auto to_tensor(const std::vector<T> &vector,
                      const iganet::Options<T> &options) {
  if (options.device() == torch::kCPU)
    return torch::from_blob(const_cast<T *>(std::data(vector)), vector.size(),
                            options)
        .detach()
        .clone()
        .requires_grad_(options.requires_grad());
  else
    return torch::from_blob(const_cast<T *>(std::data(vector)), vector.size(),
                            options.device(torch::kCPU))
        .detach()
        .clone()
        .to(options.device())
        .requires_grad_(options.requires_grad());
}
 
template <typename T, std::size_t N>
inline auto to_ArrayRef(const std::array<T, N> &array) {
  return at::ArrayRef<T>{array};
}
 
template <typename T, std::size_t... N>
inline auto concat(const std::array<T, N> &...arrays) {
  std::array<T, (N + ...)> result;
  std::size_t index{};
 
  ((std::copy_n(arrays.begin(), N, result.begin() + index), index += N), ...);
 
  return result;
}
 
template <typename T, std::size_t... N>
inline auto concat(std::array<T, N> &&...arrays) {
  std::array<T, (N + ...)> result;
  std::size_t index{};
 
  ((std::copy_n(std::make_move_iterator(arrays.begin()), N,
                result.begin() + index),
    index += N),
   ...);
 
  return result;
}
 
template <typename... Ts>
inline auto concat(const std::vector<Ts> &...vectors) {
  std::vector<std::common_type_t<Ts...>> result;
 
  (result.insert(result.end(), vectors.begin(), vectors.end()), ...);
 
  return result;
}
 
template <typename... Ts> inline auto concat(std::vector<Ts> &&...vectors) {
  std::vector<std::common_type_t<Ts...>> result;
 
  (result.insert(result.end(), std::make_move_iterator(vectors.begin()),
                 std::make_move_iterator(vectors.end())),
   ...);
 
  return result;
}
 
template <typename T>
inline constexpr auto operator+(torch::ArrayRef<T> array, T data) {
  std::vector<T> result{array.vec()};
  result.push_back(data);
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr auto operator+(std::array<T, N> array, T data) {
  std::array<T, N + 1> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = array[i];
  result[N] = data;
  return result;
}
 
template <typename T>
inline constexpr auto operator+(std::vector<T> vector, T data) {
  std::vector<T> result{vector};
  result.push_back(data);
  return result;
}
 
template <typename T>
inline constexpr auto operator+(T data, torch::ArrayRef<T> array) {
  std::vector<T> result{array.vec()};
  result.insert(result.begin(), data);
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr auto operator+(T data, std::array<T, N> array) {
  std::array<T, N + 1> result;
  result[0] = data;
  for (std::size_t i = 0; i < N; ++i)
    result[i + 1] = array[i];
  return result;
}
 
template <typename T>
inline constexpr auto operator+(T data, std::vector<T> vector) {
  std::vector<T> result{vector};
  result.insert(result.begin(), data);
  return result;
}
 
template <typename T, std::size_t N> inline constexpr auto make_array(T value) {
  std::array<T, N> result;
  result.fill(value);
  return result;
}
 
template <typename T, typename U, std::size_t N>
inline constexpr std::array<T, N> make_array(std::array<U, N> array) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = static_cast<T>(array[i]);
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr std::array<T, N> operator-(std::array<T, N> array) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = -array[i];
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr std::array<T, N> operator+(std::array<T, N> lhs,
                                            std::array<T, N> rhs) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = lhs[i] + rhs[i];
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr std::array<T, N> operator-(std::array<T, N> lhs,
                                            std::array<T, N> rhs) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = lhs[i] - rhs[i];
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr std::array<T, N> operator*(std::array<T, N> lhs,
                                            std::array<T, N> rhs) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = lhs[i] * rhs[i];
  return result;
}
 
template <typename T, std::size_t N>
inline constexpr std::array<T, N> operator/(std::array<T, N> lhs,
                                            std::array<T, N> rhs) {
  std::array<T, N> result;
  for (std::size_t i = 0; i < N; ++i)
    result[i] = lhs[i] / rhs[i];
  return result;
}
 
template <typename T, std::size_t N, std::size_t M = 1>
inline constexpr std::array<T, N - M>
remove_from_front(std::array<T, N> array) {
 
  std::array<T, N - M> result;
  for (std::size_t i = 0; i < N - M; ++i)
    result[i] = array[i + M];
  return result;
}
 
template <typename T, std::size_t N, std::size_t M = 1>
inline constexpr std::array<T, N - M> remove_from_back(std::array<T, N> array) {
 
  std::array<T, N - M> result;
  for (std::size_t i = 0; i < N - M; ++i)
    result[i] = array[i];
  return result;
}
 
} // namespace iganet::utils