計算點乘與矩陣乘法連乘計算錯誤
當我們想將 R 中的連乘（如下公式所示）修改成 Rcpp 代碼時,
|
1
|
t(X)^2 %*% X |
理論上我們只用在 .cpp 代碼中輸入下述語句即可（默認使用了 RcppEigen 庫）：
|
1
|
X.adjoint().array().square() * X.array().square(); |
但實際上這樣會會出現問題，原因是 X.adjoint().array().square() 與 X.array().square() 沒有成功轉化成 Eigen::MatrixXd 形式。雖然單獨顯示兩部分都沒問題（可以直接輸出看到），但是連乘的時候就會出錯。這時我們分開進行一個默認轉換，然后再進行矩陣乘法即可，如下所示：
|
1
2
3
4
|
Eigen::MatrixXd A, B;A = X.adjoint().array().square();B = X.array().square();A * B; |
這樣就可以解決之前的計算錯誤問題了。
以上就是RcppEigen計算點乘與矩陣乘法連乘計算錯誤解決的詳細內容，更多關于RcppEigen計算點乘與矩陣乘法的資料請關注服務器之家其它相關文章！
原文鏈接：https://kanny.blog.csdn.net/article/details/102795674
