消除中间特征数组

Question

Eigen 是否制作任何中间数组来计算 x 或 Eigen 只是将值放入 simd 寄存器并进行计算？

一般来说，如何知道 Eigen 制作了多少个中间体？

Eigen 会在循环的每个循环中为中间体分配新的内存吗？

无论如何确保本征不会产生任何中间体？它有像“EIGEN_NO_INTERMEDIATE”这样的宏吗？

#include <Eigen/Eigen>
#include <iostream>

using namespace Eigen;

template<typename T>
void fill(T& x) {
    for (int i = 0; i < x.size(); ++i) x.data()[i] = i + 1;
}

int main() {
    int n = 10;  // n is actually about 400
    ArrayXXf x(n, n);
    ArrayXf y(n);
    fill(x);
    fill(y);
    for (int i = 0; i < 10; ++i) {  // many cycles
        x = x * ((x.colwise() / y).rowwise() / y.transpose()).exp();
    }
    std::cout << x << "\n";
}

Answer 1

您可以DenseStorage像这样在构造函数中添加一个钩子：

#include <iostream>

static long int nb_temporaries;
inline void on_temporary_creation(long int size) {
  if(size!=0) nb_temporaries++;
}

// must be defined before including any Eigen header!
#define EIGEN_DENSE_STORAGE_CTOR_PLUGIN { on_temporary_creation(size); }

#define VERIFY_EVALUATION_COUNT(XPR,N) {\
    nb_temporaries = 0; \
    XPR; \
    if(nb_temporaries!=N) { std::cerr << "nb_temporaries == " << nb_temporaries << "\n"; }\
  }

#include <Eigen/Core>
using namespace Eigen;

template<typename T>
void fill(T& x) { for(int i=0; i<x.size(); ++i) x(i)= i+1; }

int main() {
    int n=10;
    ArrayXXf x(n,n); fill(x);
    ArrayXf y(n); fill(y);

    for(int i=0; i<10; ++i)
    {
        VERIFY_EVALUATION_COUNT( x = x * ((x.colwise()/y).rowwise()/y.transpose()).exp(), 0);
    }
    std::cout << x << '\n';
}

本质上，这就是 Eigen 在其测试套件中在某些时候所做的：请参阅此处了解测试套件中的原始定义，以及此处了解测试套件中的示例用法。

或者，如果您只关心中间内存分配，您可以尝试使用宏EIGEN_RUNTIME_NO_MALLOC——这将允许固定大小的表达式计算为临时值，因为它们只会在堆栈上分配。