獲取特徵矩陣元素類型的好方法

Question

我正在爲我的項目編寫一個簡單的界面，介於std::vector和特徵之間。對於簡單的矩陣 - 矩陣乘法碼：

template<typename MatrixElementType1, typename MatrixElementType2> 
    inline auto _matmat(const vector<MatrixElementType1>& Mat1, 
         const vector<MatrixElementType2>& Mat2, 
         const size_t M, const size_t N, const size_t K) 
    { 
     using MatrixType1 = Matrix<MatrixElementType1, Dynamic, Dynamic>; 
     using MatrixType2 = Matrix<MatrixElementType2, Dynamic, Dynamic>; 
     Map<const MatrixType1> m1(Mat1.data(), M, N); 
     Map<const MatrixType2> m2(Mat2.data(), N, K); 
     auto rst = m1 * m2; 
     return vector<???>(rst.data(), rst.data() + rst.size()); 

    } 

問題是我應該用什麼？ ..我知道decltype(rst(0,0))作品，但是有更加優雅的方式嗎？這似乎爲decltype(rst)::value_type本徵矩陣不工作..

Answer 1

首先，因爲m1*m2返回表達式和產品沒有結果，你需要評估它明確如果使用自動：

auto rst = (m1 * m2).eval(); 

然後，您可以使用decltype(rst)::Scalar獲得標量類型。一個更好的策略雖然可能是聲明適當類型和大小的std :: vector，並將其映射：

typedef decltype(m1*m2)::Scalar ResScalar; 
vector<ResScalar> res(M,K); 
Matrix<ResScalar,Dynamic,Dynamic>::Map(res.data(),M,K).noalias() = m1*m2; 
return res;