Eigen：具有小整数指数的系数式 pow

Question

在本征中，与

ArrayXXf a;
a = ArrayXXf::Random(1000, 10000);

在做

a = a.pow(4);

在我的电脑上大约需要 500 毫秒，而在做

a = a.square().square();

只需要大约 5 毫秒。我正在使用最近发布的 GCC 进行编译。

这是预期的行为还是我做错了什么？我希望，至少对于小整数（比如

Answer 1

使用 C++17 if constexpr 这可能是可能的，但不是。所以目前，a.pow(x) 相当于为每个i 调用std::pow(a[i],x)。