犰狳子视图问题

Question

以下片段产生编译错误：

arma::Mat<double> a(10,10,arma::fill::zeros);
arma::ucolvec w = whatever1;
whatever2 = a.rows(w).each_col() + another-col-vector;

错误是 arma::subview_elem2 没有名为 each_col 的成员。

在 Armadillo 的许多情况下，标准数组函数并不总是可用于其他函数调用的表达式或结果。显然 rows() 函数不返回一个 Mat 对象，而是一个 subview_elem2 对象，大概是为了优化。另一种方法是在 Mat 和其他内部类（例如子视图）实现的接口/纯抽象类中声明所有数组函数。除了只生成 r 值的表达式的写操作之外，似乎应该可以使所有 Armadillo 数组表达式与数组对象互换。

所以...我希望得到以下结果

a) 说明哪些方法不适用于哪些结果。 b) 最好启用所有有意义的数组方法组合。

如果没有上述，如何才能达到想要的结果，也就是对表达式求值：

a.rows(w).each_col()

??

Answer 1

一些关于犰狳的先验信息

犰狳库大量使用模板，大多数操作返回expression templates。只有当您将结果分配给变量时，才会执行实际计算。这就是为什么你不应该使用 auto 存储犰狳计算的结果。

例如，给定一些矩阵A、B 和C，类似于

auto D = A * B + C;

不会执行计算，只有表达式模板存储在D 中。另一方面，使用

arma::mat D = A * B + C;

将强制进行计算并将结果存储在D。

解决您的问题

特别是对于您的问题，a.rows(w) 之类的内容会返回 subview_elem2 类型的表达式模板（此文件在源代码 armadillo_bits/subview_elem2_bones.hpp 中定义）。这种“临时类型”没有.each_col 方法，这会导致您得到错误。一种解决方法是将a.rows(w) 的结果存储在一个变量中，但由于您对该变量不感兴趣，您可以使用.eval() 方法。 .eval() 方法强制表达式模板执行到该点的实际计算，因此随后对 .each_col 的调用将起作用。也就是替换

a.rows(w).each_col() + another-col-vector;

与

a.rows(w).eval().each_col() + another-col-vector;