我正在寻找干净的语法糖来从显式值初始化 boost::multi_array。我能想到的最好的办法是

double g[5][5] = {
{-0.0009    ,  0.003799 ,    0.00666   ,   0.00374   ,   0.00186 },
{-0.0008    ,  0.0176   ,    0.0619    ,   0.0159    ,   0.00324 },
{0.00099    ,  0.0475   ,    0.666     ,   0.0376    ,   0.00758 },
{0.00242    ,  0.02189  ,    0.0624    ,   0.0192    ,   0.0008 },
{0.00182    ,  0.00404  ,    0.00479   ,   0.00924   ,   0.00189 }};

boost::multi_array_ref<double,2> mg((double*)g,boost::extents[5][5]);

我不喜欢这样，因为它需要 2 个变量而不是 1 个，三重冗余 [5][5] 维度(大小可以通过嵌套的大括号列表推断)，以及从 double[][] 到 double* 的转换。

我确实喜欢这样一个事实，即没有从 g 执行复制到 mg，并且 g 的初始化方式是赏心悦目(即嵌套的、结构化的初始化列表，具有最少的多余字符)。

最佳答案

有几个选项可用。它们都需要一些模板魔法；它们在句法表达能力和效率上有所不同。如果boost::multi_array，生活会更轻松和 friend 们实际上提供了一些更有用的构造函数，但遗憾的是，在撰写本文时情况并非如此。

1) 使用平面 initializer_list

这个基本选项删除了一些冗余，并提供了相当不错的语法糖。它使用一个带有 initializer_list<t> 的辅助函数。 , 将其转储到 std::vector 中，并使用它首先创建一个 const_multi_array_ref然后将其深复制到 multi-array 中.

#include <boost/multi_array.hpp>
#include <cassert>
#include <initializer_list>
#include <vector>

// Helper class to determine the full number of elements in the
// multi-dimensional array
template <std::size_t... vs> struct ArraySize;
template <std::size_t v, std::size_t... vs> struct ArraySize<v, vs...>
{ static constexpr std::size_t size = v * ArraySize<vs...>::size; };
template <> struct ArraySize<>
{ static constexpr std::size_t size = 1; };

// Creates your multi_array
template <typename T, int... dims>
boost::multi_array<T, sizeof...(dims)>
makeMultiArray(std::initializer_list<T> l)
{
  constexpr std::size_t asize = ArraySize<dims...>::size;
  assert(l.size() == asize); // could be a static assert in C++14

  // Dump data into a vector (because it has the right kind of ctor)
  const std::vector<T> a(l);
  // This can be used in a multi_array_ref ctor.
  boost::const_multi_array_ref<T, sizeof...(dims)> mar(
    &a[0],
    std::array<int, sizeof...(dims)>{dims...});
  // Finally, deep-copy it into the structure we can return.
  return boost::multi_array<T, sizeof...(dims)>(mar);
}

// Usage example
auto mg = makeMultiArray<double, 5, 5>({
  -0.0009, 0.003799, 0.00666, 0.00374, 0.00186,
  -0.0008, 0.0176,   0.0619,  0.0159,  0.00324,
  0.00099, 0.0475,   0.666,   0.0376,  0.00758,
  0.00242, 0.02189,  0.0624,  0.0192,  0.0008,
  0.00182, 0.00404,  0.00479, 0.00924, 0.00189});

在此版本中，initializer_list 的适当大小只在运行时检查，但我认为在 C++14 std::initializer_list::size()将是 constexpr ，这应该允许您使用 static_assert .

• 优点:声明和维度中的冗余已消失。
• 缺点:至少创建一个拷贝；使用不易读的平面列表。

2) 从 C 数组初始化

这个更接近您的原始数组，但您需要单独定义数组——我认为您不能直接将其作为具有冗余转换的函数参数提供。优点是因为您首先构建了一个标准的 C 数组，所以您可以在 multi_array_ref 中重复使用它。而且您不需要拷贝。与第一个选项相比，您需要一些额外的构造:CArray用于从模板参数构造 C 数组类型。

// CArray<double,1,2,3>::type is double[1][2][3]
template <typename T, std::size_t... vs> struct CArray;
template <typename T, std::size_t v, std::size_t... vs> struct CArray<T, v, vs...>
{ typedef typename CArray<T, vs...>::type type[v]; };
template <typename T> struct CArray<T> { typedef T type; };

// Creates a multi_array_ref
template <typename T, int... dims>
boost::multi_array_ref<T, sizeof...(dims)>
makeMultiArray(typename CArray<T, dims...>::type l)
{
  constexpr std::size_t asize = ArraySize<dims...>::size;
  return boost::multi_array_ref<T, sizeof...(dims)>(
           reinterpret_cast<double*>(l),
           std::array<int, sizeof...(dims)>{dims...});
}

// Usage example
double g[5][5] =
  { { -0.0009, 0.003799, 0.00666, 0.00374, 0.00186 },
    { -0.0008, 0.0176,   0.0619,  0.0159,  0.00324 },
    { 0.00099, 0.0475,   0.666,   0.0376,  0.00758 },
    { 0.00242, 0.02189,  0.0624,  0.0192,  0.0008  },
    { 0.00182, 0.00404,  0.00479, 0.00924, 0.00189 } };
auto mg = makeMultiArray<double, 5, 5>(g);

• 优点:保持初始化程序的层次结构，以 boost 可读性；避免复制。
• 缺点:仍然有一些您想要摆脱的冗余。

3) 使用嵌套 initializer_list s

这个似乎更合适，但效率较低。除了上面的代码，我们还需要一个方法来构造嵌套的initializer_list s，并将它们复制到一个数组中。

// Nested initializer lists
template <typename T, std::size_t level> struct NestedInitializerList
{
  typedef std::initializer_list<typename NestedInitializerList<T, level-1>::type> type;
};
template <typename T> struct NestedInitializerList<T, 1>
{
  typedef std::initializer_list<T> type;
};

// Helpers which fill the array from a nested initializer_list
template <typename T>
void fillArray(const T& l, typename CArray<T>::type& a)
{
  a = l;
}
template <typename T, int dim, int... dims>
void fillArray(typename NestedInitializerList<T, sizeof...(dims)+1>::type l,
               typename CArray<T, dim, dims...>::type& a)
{
  assert(l.size() == dim); // could be a static assert in C++14
  int i=0;
  for (auto it = l.begin(); it != l.end(); ++it, ++i)
  {
    fillArray<T, dims...>(*it, a[i]);
  }
}

// Creates your multi_array
template <typename T, int... dims>
boost::multi_array<T, sizeof...(dims)>
makeMultiArray(typename NestedInitializerList<T, sizeof...(dims)>::type l)
{
  typename CArray<T, dims...>::type a; // Multidimensional C array
  fillArray<T, dims...>(l, a);         // Fill from l
  // Turn into multi_array_ref.
  boost::const_multi_array_ref<T, sizeof...(dims)> mar(
    reinterpret_cast<const double*>(a),
    std::array<int, sizeof...(dims)>{dims...});
  // Finally, deep-copy it into the structure we can return.
  return boost::multi_array<T, sizeof...(dims)>(mar);
}

// Usage example
auto mg = makeMultiArray<double, 5, 5>(
  { { -0.0009, 0.003799, 0.00666, 0.00374, 0.00186 },
    { -0.0008, 0.0176,   0.0619,  0.0159,  0.00324 },
    { 0.00099, 0.0475,   0.666,   0.0376,  0.00758 },
    { 0.00242, 0.02189,  0.0624,  0.0192,  0.0008  },
    { 0.00182, 0.00404,  0.00479, 0.00924, 0.00189 } });

最后一个来自 this article 的一点灵感。 .

• Pro:用法完全如您所愿；无冗余，清晰的层次初始化结构
• 缺点:fillArray()例程是递归的，因此效率较低(我预计编译器也无法对其进行优化)。

关于c++ - 用实际元素初始化 boost::multi_array 的干净方法，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/35209121/