对齐数据类型 Eigen::Matrix 的数组或向量声明答案

【问题标题】：Array or vector declaration for aligned datatype Eigen::Matrix对齐数据类型 Eigen::Matrix 的数组或向量声明
【发布时间】：2014-06-24 14:46:14
【问题描述】：

我有一个类ClusterNode，其中包含一个Eigen::Matrix4d 作为类变量。我还有一个函数numNodes()，它告诉我ClusterNodes 的总数，以便我可以将它们收集到某种数组、列表或向量中。

但是，Eigen::Matrix4d 是对齐的，这意味着我不能将这种类型的对象存储在 std::vector<ClusterNode> 中，根据这个问题的答案 error C2719: '_Val': formal parameter with __declspec(align('16')) won't be aligned?

作为替代方案，我尝试使用数组。

但是，我做不到

const int n = numNodes();
ClusterNode array [n];

因为函数的返回值不被视为常量。

我还有什么其他选择？

【问题讨论】：

【解决方案1】：

您可以使用 Eigen 的对齐分配器。详情见this page，但基本上你可以做到：

std::vector<ClusterNode, Eigen::aligned_allocator<ClusterNode> >

也不要忘记为你的类重载 operator new。同样，所有细节都在 Eigen 的 documentation 中。

【讨论】：

【解决方案2】：

你可以尝试使用 memalign 或 posix_memalign 来分配和对齐一个可变的内存空间

【讨论】：

【解决方案3】：

我最终使用了 Alex Faber 的解决方案

ClusterNode* array = new ClusterNode[n];

但我想补充一点，在这种情况下，类还需要返回对齐的指针。它只会在部分时间以一种非常烦人且无法追踪的方式引起问题，但经过多次寻找错误后，我终于收到了一条有用的错误消息，将我指向Eigen's documentation。

文档摘要：Eigen 需要 16 字节对齐的内存。以下情况会导致悲伤。

如果我使用 Alex 提出的动态数组解决方案（ClusterNode 作为变量），我还需要使用 Eigen 提供的宏来确保对象对齐。例如，

class ClusterNode
{
    Eigen::Matrix4d M;

public:
    EIGEN_MAKE_ALIGNED_OPERATOR_NEW
}

如果不包含宏，数组声明有时会静默失败。

【讨论】：