特徵矩陣庫係數模運算

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