Eigen库：向模板函数传递Block引用时出现编译错误

我太懂这种两难的感觉了——用Eigen写模板函数时，要么临时对象传不进去报编译错，要么把参数改成const后又没法修改矩阵，简直卡得死死的！先帮你把问题场景还原出来，方便咱们对着分析：

比如你可能写了类似这样的代码：

#include <Eigen/Core>

template<typename T>
void fcn(T& mat) {
    mat(0,0) = 1.0; // 想修改矩阵元素
}

int main() {
    Eigen::MatrixXd mat = Eigen::MatrixXd::Zero(2,2);
    fcn(mat); // 这里调用没问题
    // 但传临时对象就炸了
    fcn(Eigen::MatrixXd::Zero(2,2)); // 编译错误：无法将临时值绑定到非const左值引用
}

问题根源 #

编译器报错的核心是C++的基础规则：临时对象是右值，非const左值引用（T&）只能绑定到可修改的左值。但如果把参数改成const T& mat，函数里的mat(0,0) = 1.0又会因为const限制报错，确实陷入了两难。

最优解决方案：用Eigen::Ref搞定一切 #

Eigen专门为这种场景提供了Eigen::Ref，它是一个智能包装类，既能灵活接受左值、临时对象（右值），又允许你修改参数（只要不是const版本）。

方案1：直接用Eigen::Ref作为函数参数（非模板场景）

如果你的函数只针对特定类型的Eigen矩阵（比如MatrixXd），直接这么写就行：

#include <Eigen/Core>

void fcn(Eigen::Ref<Eigen::MatrixXd> mat) {
    mat(0,0) = 1.0; // 正常修改完全没问题
}

int main() {
    Eigen::MatrixXd mat = Eigen::MatrixXd::Zero(2,2);
    fcn(mat); // 左值调用正常
    fcn(Eigen::MatrixXd::Zero(2,2)); // 临时对象也能正常编译执行
}

方案2：结合模板保持通用性（支持多种Eigen类型）

如果想让函数支持任意Eigen矩阵/向量类型（比如Matrix<float,2,2>、Vector3d等），可以把Eigen::Ref和模板参数结合：

#include <Eigen/Core>

template<typename Derived>
void fcn(Eigen::Ref<Eigen::MatrixBase<Derived>> mat) {
    mat(0,0) = 1.0; // 支持修改任意兼容的Eigen表达式
}

int main() {
    Eigen::MatrixXd mat_d = Eigen::MatrixXd::Zero(2,2);
    fcn(mat_d); // double类型矩阵正常

    Eigen::Matrix<float,2,2> mat_f = Eigen::Matrix<float,2,2>::Zero();
    fcn(mat_f); // float类型矩阵也支持

    fcn(Eigen::MatrixXd::Zero(2,2)); // 临时对象照样能调用
}

为什么Eigen::Ref能解决问题？ #

Eigen::Ref的本质是一个适配器：

• 当传入左值矩阵时，它直接引用原对象，没有额外开销；
• 当传入临时对象（右值）时，它会自动创建一个可修改的临时拷贝，让你能正常修改，同时遵循C++的规则；
• 它还支持Eigen的块操作（比如mat.block(0,0,2,2)），适用性极强。

如果你的函数不需要修改参数，那直接用const Eigen::Ref<const Eigen::MatrixBase<Derived>>或者const Eigen::MatrixBase<Derived>&就好，既安全又高效。

内容的提问来源于stack exchange，提问作者T Devel