R語言RcppEigen計算點(diǎn)乘與矩陣乘法連乘算法錯誤解決
計算點(diǎn)乘與矩陣乘法連乘計算錯誤
當(dāng)我們想將 R 中的連乘(如下公式所示)修改成 Rcpp 代碼時,
t(X)^2 %*% X
理論上我們只用在 .cpp 代碼中輸入下述語句即可(默認(rèn)使用了 RcppEigen 庫):
X.adjoint().array().square() * X.array().square();
但實際上這樣會會出現(xiàn)問題,原因是 X.adjoint().array().square() 與 X.array().square() 沒有成功轉(zhuǎn)化成 Eigen::MatrixXd 形式。雖然單獨(dú)顯示兩部分都沒問題(可以直接輸出看到),但是連乘的時候就會出錯。這時我們分開進(jìn)行一個默認(rèn)轉(zhuǎn)換,然后再進(jìn)行矩陣乘法即可,如下所示:
Eigen::MatrixXd A, B; A = X.adjoint().array().square(); B = X.array().square(); A * B;
這樣就可以解決之前的計算錯誤問題了。
以上就是RcppEigen計算點(diǎn)乘與矩陣乘法連乘計算錯誤解決的詳細(xì)內(nèi)容,更多關(guān)于RcppEigen計算點(diǎn)乘與矩陣乘法的資料請關(guān)注腳本之家其它相關(guān)文章!
相關(guān)文章
R語言中平均值、中位數(shù)和模式知識點(diǎn)總結(jié)
在本篇文章里小編給大家整理的是一篇關(guān)于R語言中平均值、中位數(shù)和模式知識點(diǎn)總結(jié)內(nèi)容,有興趣的朋友們跟著學(xué)習(xí)下。2021-05-05
R語言 install.packages 無法讀取索引的解決方案
這篇文章主要介紹了R語言 install.packages 無法讀取索引的解決方案,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
R語言修改下載安裝包install.package的默認(rèn)存儲路徑的操作方法
這篇文章主要介紹了R語言修改下載安裝包install.package的默認(rèn)存儲路徑的操作方法,本文給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價值,需要的朋友可以參考下2021-03-03

