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