0

我有 Eigen::vector3d 类型的 std::vector ，当我使用 Microsoft Visual Studio 2012 编译此代码时，打开/Qvec-report:2标志以报告矢量化详细信息。它显示由于 msdn 页面上指定的原因 1304（循环包含不同类型的分配）而未矢量化循环 - https://msdn.microsoft.com/en-us/library/jj658585.aspx

我的代码如下：

#include <iostream>
#include <vector>
#include <time.h>
#include<Eigen/StdVector>
int main(char *argv[], int argc)
{
    int tempSize=100;
/** I am aligning these vectors as specfied on http://eigen.tuxfamily.org/dox/group__TopicStlContainers.html */
    std::vector<Eigen::Vector3d,Eigen::aligned_allocator<Eigen::Vector3d>> eiVec(tempSize);
    std::vector<Eigen::Vector3d,Eigen::aligned_allocator<Eigen::Vector3d>> eiVec1(tempSize);
    std::vector<Eigen::Vector3d,Eigen::aligned_allocator<Eigen::Vector3d>> eiVec2(tempSize);

    for(int i=0;i<100;i++)
    {
        eiVec1[i] = Eigen::Vector3d::Zero();
        eiVec2[i] = Eigen::Vector3d::Zero();
    }

    Eigen::Vector3d *eV = &eiVec.front();
    const Eigen::Vector3d *eV1 = &eiVec1.front();
    const Eigen::Vector3d *eV2 = &eiVec2.front();

/** Below loop is not vectorized due to code 1304  */
    for(int i=0;i<100;i++)
    {
        eV[i] = eV1[i] - eV2[i];
    }
    return 0;
}

所以我无法理解如何告诉编译器，这是固定大小的数据，固定大小的分配将在这里发生。

4

2 回答 2

1

来自本征文档

在 Eigen 中，诸如算术运算符operator+本身不执行任何计算，它们只是返回一个描述要执行的计算的“表达式对象”。实际计算发生在稍后，当整个表达式被评估时，通常在 operator=

所以eV1[i] - eV2[i]不返回Eigen::Vector3d- 所以原因 1304 适用。你也许可以施放它——但我不会打赌。

于 2016-10-17T11:59:08.963 回答
0

考虑将eiVec1andeiVec12更改为 amatrix<double, 3, 100>并进行减法eiVec = eiVec1 - eiVec2

于 2017-06-26T00:08:36.900 回答