C++/STL：std::transform 给定的步幅？

Question

我有一个包含 Nd 数据的一维数组，我想用 std::transform 或 std::for_each 有效地遍历它。

unigned int nelems;
unsigned int stride=3;// we are going to have 3D points
float *pP;// this will keep xyzxyzxyz...
Load(pP);
std::transform(pP, pP+nelems, strMover<float>(pP, stride));//How to define the strMover??

Answer 1

好吧，我决定使用 for_each 而不是转换，欢迎任何其他决定：

generator<unsigned int> gen(0, 1);
            vector<unsigned int> idx(m_nelem);//make an index
            std::generate(idx.begin(), idx.end(),gen);
            std::for_each(idx.begin(), idx.end(), strMover<float>(&pPOS[0],&m_COM[0],stride));

从

template<class T> T op_sum (T i, T j) { return i+j; }
template<class T> 
class strMover
    {
    T *pP_;
    T *pMove_;
    unsigned int stride_;
    public:
        strMover(T *pP,T *pMove, unsigned int stride):pP_(pP), pMove_(pMove),stride_(stride)
            {}
        void operator() ( const unsigned int ip )
            {
            std::transform(&pP_[ip*stride_], &pP_[ip*stride_]+stride_, 
                pMove_, &pP_[ip*stride_], op_sum<T>);
            }
    };

第一眼看，这是一个线程安全的解决方案。

Answer 2

使用升压适配器。你可以从中得到迭代器。唯一的缺点是编译时间。

vector<float> pp = vector_load(pP);
boost::for_each(pp|stride(3)|transformed(dosmtn()));

Answer 3

答案不是改变 strMover，而是改变你的迭代器。定义一个新的迭代器类，它包装一个 float * 但在调用 operator++ 时向前移动 3 个位置。

您可以使用 boost 的 排列迭代器 并使用一个非严格排列，仅包括您感兴趣的范围。

如果您尝试推出自己的迭代器，则会遇到一些问题：为了严格遵守标准，您需要仔细考虑这样的正确的“结束”迭代器stride 迭代器是这样的，因为天真的实现会愉快地跨过允许的“one-past-the-end”，到达远远超过数组末尾的模糊区域，指针永远不应该进入该区域，因为担心 鼻恶魔。

但我不得不问：为什么首先要将 3d 点数组存储为 float 数组？只需定义一个 Point3D 数据类型并创建一个数组即可。简单多了。

Answer 4

这太糟糕了，人们告诉你改用跨步迭代器。除了无法通过这种方法使用标准库中的函数对象之外，还使编译器通过使用这样的拐杖来生成多核或 sse 优化变得非常非常复杂。
寻找“跨步迭代器”以获得正确的解决方案，例如在 C++ 食谱中。

回到原来的问题...使用 valarray 和 stride 来模拟多维数组。