Question

我有一个Eigen::MatrixXd，我希望通过应用组件功能来修改其所有元素。例如：

MatrixXd m = ...;

for each m[i][j]:
  m[i][j] = exp(m[i][j]);

有没有办法实现这个结果？

Answer 1

是的，使用Eigen::MatrixBase<>::unaryExpr()成员函数。例如：

#include <cmath>
#include <iostream>

#include <Eigen/Core>

double Exp(double x) // the functor we want to apply
{
    return std::exp(x);
}

int main()
{
    Eigen::MatrixXd m(2, 2);
    m << 0, 1, 2, 3;
    std::cout << m << std::endl << "becomes: ";
    std::cout << std::endl << m.unaryExpr(&Exp) << std::endl;
}

Answer 2

vsoftco的答案很通用，对自定义函数很有用。但是，对于许多常用功能，有一种更简单的方法。调整他的示例我们可以使用array s，它看起来像这样：

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

int main()
{
    Eigen::MatrixXd m(2, 2);
    m << 0, 1, 2, 3;
    std::cout << m << "\nbecomes:\n";
    std::cout << m.array().exp() << std::endl;
    return 0;
}

Answer 3

FWIW，在C ++ 11及更高版本中，也可用于lambda函数。

#include <cmath>
#include <iostream>

#include <Eigen/Core>

int main()
{
    Eigen::MatrixXd m(2, 2);
    m << 0, 1, 2, 3;
    std::cout << m << std::endl << " ->  " 
    std::cout << m.unaryExpr([](double x){return x + 1}) << std::endl;
}

Answer 4

@vsoftco的答案使我获得了99％的解决方案，但是由于某种原因，将&Exp传递给.unaryExpr()给了我编译错误（g ++，c + 11，Eigen 3.3.5给了与base type ‘double (*)(double)’ fails to be a struct or class type有关的错误。

但是，我发现创建一个std::function对象并传递该对象可以解决此问题。复制@vsoftco的示例：

#include <cmath>
#include <iostream>

#include <Eigen/Core>

double Exp(double x) 
{
    return std::exp(x);
}

int main()
{
    Eigen::MatrixXd m(2, 2);
    m << 0, 1, 2, 3;
    std::function<double(double)> exp_wrap = Exp; //added to @vsoftco's answer
    std::cout << m << std::endl << "becomes: ";
    std::cout << std::endl << m.unaryExpr(exp_wrap) << std::endl; //and used here
}

我不确定使用std::function对象（或std::ptr_fun）与传递&Exp相比会产生多少开销，但是如果没有这些替代方法，我将无法使其正常工作。

欢呼

将函数应用于所有Eigen矩阵元素

4 个答案: