[Eigen中文文档] 常见的陷阱

简介: 本文将介绍一些Eigen常见的陷阱

文档总目录

英文原文(Common pitfalls)

针对模板方法的编译错误

详见下一节 C++中的template和typename关键字

混叠

详见 [Eigen中文文档] 混叠

对齐问题(运行时断言)

Eigen进行了显式向量化,虽然这受到许多用户的赞赏,但在某些特殊情况下,数据对齐会出现问题。实际上,在C++17之前,C++没有足够好的支持显式数据对齐。在这种情况下,程序会出现断言失败(即“受控崩溃”),并显示一条消息,引导你参考此页面。它包含有关如何处理该问题的每个已知原因的详细信息。

如果你不关心矢量化并且不想处理这些对齐问题,可以阅读如何避免这一情况

C++11 和 auto 关键字

简而言之:不要将 auto 关键字与 Eigen 表达式一起使用,除非你 100% 确定自己在做什么。特别是,不要使用 auto 关键字来替代 Matrix<> 类型。这是一个例子:

MatrixXd A, B;
auto C = A*B;
for(...) {
    ... w = C * v;  ...}

在这个例子中,C 的类型不是MatrixXd,而是表示矩阵乘积并存储对 AB 的引用的抽象表达式。因此,A*B 的乘积将在 for 循环的每次迭代中执行多次。此外,如果 AB 的系数在迭代过程中发生改变,那么 C 将评估为不同的值,就像以下示例一样:

MatrixXd A = ..., B = ...;
auto C = A*B;
MatrixXd R1 = C;
A = ...;
MatrixXd R2 = C;

因此我们最终得到 R1R2

如下是导致段错误的另一个示例:

auto C = ((A+B).eval()).transpose();
// do something with C

问题在于 eval() 返回一个临时对象(在本例中为MatrixXd),然后由Transpose<>表达式引用。然而,这个临时对象在第一行之后立即被删除,然后C表达式引用了一个已经不存在的对象。一个可能的解决方法是在整个表达式上应用eval()

auto C = (A+B).transpose().eval();

当 Eigen 自动计算子表达式时,可能会出现相同的问题,如下例所示:

VectorXd u, v;
auto C = u + (A*v).normalized();
// do something with C

这里,normalized() 方法必须评估耗费资源的乘积 A*v 以避免评估两次。同样,一种可能的修复方法是对整个表达式调用 .eval()

auto C = (u + (A*v).normalized()).eval();

在本例中,C 将是常规 VectorXd 对象。请注意,当底层表达式已经是普通 Matrix<> 时,DenseBase::eval() 足够智能,可以避免复制。

头文件问题(编译失败)

对于所有的库,都必须查看文档以确定要包含哪个头文件。Eigen也是如此,但略有不同的是:对于Eigen而言,一个类中的方法可能需要比类本身更多的# include。例如,如果你想在向量上使用 cross() 方法(它计算叉积),你需要:

#include<Eigen/Geometry>

三元运算符

简而言之:避免将三元运算符 (COND ? THEN : ELSE) 与 Eigen 的 THENELSE 语句表达式一起使用。要了解原因,让我们考虑以下示例:

Vector3f A;
A << 1, 2, 3;
Vector3f B = ((1 < 0) ? (A.reverse()) : A);

这个例子将返回B = 3, 2, 1。原因是在c++中,ELSE语句的类型根据THEN表达式的类型推断,以使两者的类型匹配。由于THEN是一个Reverse<Vector3f>,因此ELSE语句A被转换为一个Reverse<Vector3f>,因此编译器生成:

Vector3f B = ((1 < 0) ? (A.reverse()) : Reverse<Vector3f>(A));

在这种非常特殊的情况下,解决方法是为 THEN 语句调用 A.reverse().eval(),但最安全、最快的方法实际上是避免使用 Eigen 表达式的三元运算符并使用 if/else 构造。

按值传递

如果不知道为什么 Eigen 的值传递是错误的,请先阅读此页

虽然可能非常小心并确保所有显式使用 Eigen 类型的代码都是按引用传递的,但必须注意在编译时定义参数类型的模板。

如果一个模板有一个以按值方式传递参数的函数,并且相关的模板参数最终成为Eigen类型,那么你当然会遇到与明确定义函数以引用方式传递Eigen类型相同的对齐问题。

使用Eigen类型与其他第三方库或甚至STL一起使用可能会遇到同样的问题。例如,boost::bind 使用按值传递来存储返回的函数对象中的参数。

至少有两种方法可以解决这个问题:

  • 如果你传递的值保证在函数对象的生命周期内存在,你可以使用boost::ref()将其包装,然后传递给boost::bind。通常,对于栈上的值,这不是一个解决方案,因为如果该函数对象传递到更低的作用域或独立作用域,该对象可能已经被删除,因此在尝试使用它时会出现问题。
  • 另一种选择是让函数采用引用计数指针(如 boost::shared_ptr)作为参数。这避免了需要管理所传递对象的生命周期。

具有布尔系数的矩阵

目前使用具有布尔系数的 Matrix 的行为不一致,并且可能在 Eigen 的未来版本中发生变化,因此请谨慎使用!

这种不一致的一个简单例子是:

template<int Size>
void foo() {
   
    Eigen::Matrix<bool, Size, Size> A, B, C;
    A.setOnes();
    B.setOnes();

    C = A * B - A * B;
    std::cout << C << "\n";
}

因为调用 foo<3>() 会打印零矩阵,而调用 foo<10>() 会打印单位矩阵。

相关文章
|
存储 编译器
[Eigen中文文档] 深入了解 Eigen - 类层次结构
本页面介绍了Eigen类层次结构中 Core 类的设计及其相互关系。一般用户可能不需要关注这些细节,但对于高级用户和Eigen开发人员可能会有用。
285 0
|
存储 编译器
|
存储 C语言 C++
|
存储 缓存
[Eigen中文文档] 深入了解 Eigen - 惰性求值与混叠(Aliasing)
Eigen具有智能的编译时机制,可以实现惰性求值并在适当的情况下删除临时变量。它会自动处理大多数情况下的混叠问题,例如矩阵乘积。自动行为可以通过使用MatrixBase::eval()和MatrixBase::noalias()方法手动覆盖。
296 0
|
存储 算法 NoSQL
[Eigen中文文档] 稀疏矩阵操作
在许多应用中(例如,有限元方法),通常要处理非常大的矩阵,其中只有少数系数不为零。在这种情况下,可以通过使用仅存储非零系数的特殊表示来减少内存消耗并提高性能。这样的矩阵称为稀疏矩阵。
461 0
|
存储 索引
[Eigen中文文档] 扩展/自定义Eigen(三)
本页面针对非常高级的用户,他们不害怕处理一些Eigen的内部细节。在大多数情况下,可以通过使用自定义一元或二元函数避免使用自定义表达式,而极其复杂的矩阵操作可以通过零元函数(nullary-expressions)来实现,如前一页所述。 本页面通过示例介绍了如何在Eigen中实现新的轻量级表达式类型。它由三个部分组成:表达式类型本身、包含有关表达式编译时信息的特性类和评估器类,用于将表达式评估为矩阵。
134 1
[Eigen中文文档] 在 CMake 项目中使用 Eigen
Eigen提供了CMake(CMake 3.0或更高版本)支持,使得该库可以轻松地在CMake项目中使用。
692 1
|
编译器 C语言
[Eigen中文文档] 断言
宏 eigen_assert默认定义为 eigen_plain_assert。我们使用 eigen_plain_assert而不是 assert来解决 GCC <= 4.3 的已知错误。基本上,eigen_plain_assert就是断言。
126 0
|
并行计算 算法 安全
[Eigen中文文档] Eigen 和多线程
某些 Eigen 算法可以利用硬件中存在的多个内核。
436 0
|
编译器 Linux C语言
[Eigen中文文档] 从入门开始...
这是一个非常简短的Eigen入门文章。该文章有两层目的。对于想要尽快开始编码的人来说,该文章是对Eigen库的最简单介绍。你可以把该文章作为教程的第一部分,这更加详细的解释了Eigen库。看完这个教程后可以继续阅读 The Matrix class教程。
529 0