Eigen：按行计算矩阵的范数，比在向量上迭代计算它们要慢

Question

我正在尝试使用 Eigen 加速向量范数的计算。

尝试了两种方式：

方式一：将向量分别存储在一个数组中（std::vector）

std::vector<Eigen::Matrix<double, 1, VECTOR_SIZE>> my_vectors(num_vectors);
Eigen::Matrix<double, Eigen::Dynamic, 1> norms(num_vectors, 1);
for (int i = 0 ; i < my_vectors.size() ; i++) {
    norms(i, 0) = my_vectors[i].norm();
}

方式2：将向量存储为同一矩阵的行并使用行范数

Eigen::Matrix<double, Eigen::Dynamic, VECTOR_SIZE> my_vectors(num_vectors, VECTOR_SIZE);
Eigen::Matrix<double, Eigen::Dynamic, 1> norms = my_vectors.rowwise().norm();

我很惊讶地发现 Way1 很慢，而 Way2 更慢。

我做错了吗？有没有办法更快地计算规范？

在我做的另一个测试中，向量减法也是如此。对单独的向量进行迭代，从每个向量中减少相同的向量，比将所有向量存储为矩阵行并使用 .rowwise() - vector_to_subtract

更快

Answer 1

来自https://eigen.tuxfamily.org/dox/group__TopicStorageOrders.html

如果未指定存储顺序，则 Eigen 默认以列为主存储条目。如果使用方便的 typedef（Matrix3f、ArrayXXd 等）之一，情况也是如此。

我最好的猜测是内存访问问题：rowwise() 可能在内存中进行了一堆“跳过”，即不读取连续区域 - 而单独存储每一行​​不会遇到这个问题。

警告：虽然我是出于好意写下我的答案，但我无法自己测试它，而且显然它会进一步减慢速度。感谢SomethingSomething 测试这个想法。