R語言RcppEigen計算點乘與矩陣乘法連乘算法錯誤解決
更新時間：2021年11月06日 11:29:49 作者：Kanny廣小隸
這篇文章主要為大家介紹了RcppEigen計算點乘與矩陣乘法時發(fā)生連乘計算錯誤的解決方式，有需要的朋友可以借鑒參考下，希望能夠有所幫助
計算點乘與矩陣乘法連乘計算錯誤
當我們想將 R 中的連乘（如下公式所示）修改成 Rcpp 代碼時,
t(X)^2 %*% X
理論上我們只用在 .cpp 代碼中輸入下述語句即可（默認使用了 RcppEigen 庫）：
X.adjoint().array().square() * X.array().square();
但實際上這樣會會出現(xiàn)問題，原因是 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計算點乘與矩陣乘法的資料請關注腳本之家其它相關文章！
相關文章
這篇文章主要介紹了基于R語言for循環(huán)的替換方案，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
這篇文章主要介紹了R語言中aggregate?函數(shù)的相關資料,aggregate函數(shù)是數(shù)據(jù)處理中常用到的函數(shù)，具有強大的功能，可以按照要求把數(shù)據(jù)打組聚合，然后對聚合以后的數(shù)據(jù)進行加和、求平均等各種操作，需要的朋友可以參考下2023-04-04
這篇文章主要為大家介紹了R語言繪制line?plot線圖的示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進步，早日升職加薪2022-02-02
R語言繪圖數(shù)據(jù)可視化Ridgeline plot山脊圖畫法
這篇文章主要為大家介紹了R語言繪圖數(shù)據(jù)可視化Ridgeline plot山脊圖畫法的示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助2022-02-02
R語言-如何循環(huán)讀取excel并保存為RData
這篇文章主要介紹了R語言循環(huán)讀取excel并保存為RData的操作，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
這篇文章主要介紹了R語言常用命令集合，文中整理的都是常用的一些命令，有感興趣的可以多學習下2021-03-03
R語言 實現(xiàn)將factor轉換成numeric方法
這篇文章主要介紹了R語言 實現(xiàn)將factor轉換成numeric方法，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-03-03
這篇文章主要介紹了R語言 data.frame()命令的使用說明，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
這篇文章主要介紹了R語言判斷語句的使用詳解，幫助大家更好的理解和學習使用R語言，感興趣的朋友可以了解下2021-03-03
R語言ComplexHeatmap繪制復雜熱圖heatmap
ComplexHeatmap可以繪制很復雜的熱圖，能滿足日常以及文章所需，本文介紹一下文章常見的復雜熱圖的繪制方式，含代碼，文章比較長，建議耐心看一下，也許就有你需要的熱圖2022-06-06