特征向量对数错误使用不完整类型无效

Question

我正在尝试使用 Eigen 库计算向量的元素自然对数，这是我的代码：

#include <Eigen/Core>
#include <Eigen/Dense>

void function(VectorXd p, VectorXd q) {
    VectorXd kld = p.cwiseQuotient(q);
    kld = kld.log();
    std::cout << kld << std::endl;
}

但是在编译时

g++ -I eigen_lib -std=c++11 -march=native test_eigen.cpp -o test_eigen

我得到

test_eigen.cpp:15:23: error: invalid use of incomplete type ‘const class Eigen::MatrixLogarithmReturnValue<Eigen::Matrix<double, -1, 1> >’ kld = kld.log();

我错过了什么？

Answer 1

VectorXd::log()是MatrixBase<...>::log()计算方阵的矩阵对数。如果你想要元素对数，你需要使用数组功能：

kld = kld.array().log();
// or:
kld = log(kld.array());

如果您的所有操作都是按元素进行的，请考虑使用ArrayXd代替VectorXd：

void function(const Eigen::ArrayXd& p, const Eigen::ArrayXd& q) {
    Eigen::ArrayXd kld = log(p/q);
    std::cout << kld << std::endl;
}