从multi_array中提取子阵列

我在尝试从multi_array中提取子数组。对于这个演示，我们假设没有折叠维度（即子数组的维度与原始数组相同）。我认为我正在用正确的请求扩展构建一个视图（虽然很笨拙......），但是现在如何将数据从请求的索引（又称为视图的所有索引）复制到子数组中？这里是一个轮廓：从multi_array中提取子阵列

#include <boost/multi_array.hpp> 

const unsigned int Dimension = 3; 
using ArrayType = boost::multi_array<double, Dimension>; 
using IndexType = boost::array<ArrayType::index, Dimension>; 

ArrayType ExtractSubGrid(const ArrayType& array, const typename boost::detail::multi_array::index_gen<Dimension, Dimension>& indices) 
{ 
    typename ArrayType::template const_array_view<Dimension>::type view = array[indices]; 

    IndexType subArraySize; 
    for(size_t dim = 0 ; dim < Dimension; ++dim) { 
     subArraySize[dim] = indices.ranges_[dim].finish_ - indices.ranges_[dim].start_; 
    } 

    ArrayType subArray; 
    subArray.resize(subArraySize); 

    // How to do this? 
    //subArray.data() = view.data(); 

    return subArray; 
} 

int main() 
{ 
    ArrayType myArray(IndexType({{3,4,2}})); 

    boost::detail::multi_array::index_gen<3,3> indices = boost::indices[ArrayType::index_range(0,2)][ArrayType::index_range(1,3)][ArrayType::index_range(0,4)]; 

    ArrayType subGrid = ExtractSubGrid(myArray, indices); 

    return 0; 
}

来源

2016-02-29 David Doria

我预计'子阵=视图;'更逻辑接近预期的接口。没时间现在通过 – sehe

@seow哇，测试东西，编译，而且正是我想要的。现在只需测试一下正确的数据就可以了。 –

会是这样的帮助：

std::copy(view.begin(), view.end(), subArray.begin());

来源

2016-02-29 14:56:38

只要我们确定按照正确的顺序复制索引，那会很好。我试图通过填充一些数据来测试，但是看起来很像这样简单的错误：'unsigned int i = 0; （auto itr = myArray.begin（）; itr！= myArray.end（）; ++ itr）{itr = data [i]; i ++; }'？它似乎没有将'itr'解引用到一个可变的double？ –

这是因为'const' ness？ –

我在循环中尝试了'ArrayType :: iterator'而不是'auto'，以确保它没有获得const迭代器，并且它仍然具有相同的编译器错误。它们深入subarray.hpp（'other'中成员num_dimensions的请求，它是非类类型'const double'等）。 –

从multi_array中提取子阵列

回答

相关问题