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