C++ 将函数应用于所有特征矩阵元素

我有一个

Eigen:：MatrixXd

，我想通过应用函数组件来修改它的所有元素。例如：

MatrixXd m = ...;

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

有没有办法实现此结果？

是的，使用成员函数。例如：

#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;
}

#包括
#包括
#包括
double Exp（double x）//我们要应用的函子
{
返回std：：exp（x）；
}
int main（）
{
本征：：矩阵m（2,2）；
m的答案非常笼统，适用于自定义函数。但是，对于许多常用函数，有一种更简单的方法。根据他的示例，我们可以使用
数组
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;
}

#包括
#包括
int main（）
{
本征：：矩阵m（2,2）；
m@vsoftco的答案在这个问题上给了我99%的答案，但由于某种原因，将
&Exp
传递给
。unaryExpr（）
给了我编译错误（g++，c+11，Eigen 3.3.5给出了与以下相关的错误：
基类型'double（*）（double）'不能成为结构或类类型
）

但是，我发现创建一个
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
}

#包括
#包括
#包括
双经验（双x）
{
返回std：：exp（x）；
}
int main（）
{
本征：：矩阵m（2,2）；
mFWIW，在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;
}

#包括
#包括
#包括
int main（）
{
本征：：矩阵m（2,2）；
m@XingShi No，因为有重载，并且
unaryExpr
无法单独从
std:：exp
推断出函子签名（即，有多个候选项）。如果您通过指定函子类型来“帮助”编译器，如
m.unaryExpr（&std:：exp）
，它将起作用。根据API引用，必须使用std:：ptr_func来转换函数中函数的指针，该指针作为参数传递给unaryExpr.m.unaryExp（std:：ptr_func（&Exp））“亚历克斯感谢您的评论。但是，这不是必须的，API说：<代码> STD:：pTrfFooY STD:：pTrfFooS >在C++ 11中被禁止，参见。这种方法似乎可以产生更高效（读取矢量化）的机器代码