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