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