特征矩阵库系数模运算

Question

在我正在研究的项目中的一个函数中，我需要在除以给定数字时查找我的特征库矩阵的每个元素的其余部分。这里是Matlab相当于我想要做的：

mod(X,num) 

其中X是除数矩阵和num是除数。

实现此目的的最简单方法是什么？

Answer 1

可以使用C++ 11拉姆达与unaryExpr：

MatrixXi A(4,4), B; 
A.setRandom(); 
B = A.unaryExpr([](const int x) { return x%2; }); 

或：

int l = 2; 
B = A.unaryExpr([&](const int x) { return x%l; }); 

Answer 2

为了完整起见，另一种解决方案将是：

1. 变换X到一个特征数组（用于coeffwise操作），
2. 应用模数fo rmula a%b = a - (b * int(a/b))

C++代码，返回一个本征数组：

auto mod_array = X.array() - (num * (X.array()/num)); 

C++代码来获得的矩阵：

auto mod_matrix = (X.array() - (num * (X.array()/num))).matrix(); 

注意，括号是重要，特别是在(X.array()/num)作为本征将优化(num * X.array()/num)到X.array()这不是我们所期望的。

特征数组的第一个版本比使用unaryExpr的版本更快。 矩阵的第二个版本与unaryExpr的版本大致相同。

如果X包含浮点数，则需要在(X.array()/num)投X.array()为int