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() 沒有成功轉(zhuǎn)化成 Eigen::MatrixXd 形式。雖然單獨顯示兩部分都沒問題（可以直接輸出看到），但是連乘的時候就會出錯。這時我們分開進行一個默認轉(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語言for循環(huán)的替換方案
這篇文章主要介紹了基于R語言for循環(huán)的替換方案，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧
2021-04-04
R語言中aggregate?函數(shù)詳解
這篇文章主要介紹了R語言中aggregate?函數(shù)的相關(guān)資料,aggregate函數(shù)是數(shù)據(jù)處理中常用到的函數(shù)，具有強大的功能，可以按照要求把數(shù)據(jù)打組聚合，然后對聚合以后的數(shù)據(jù)進行加和、求平均等各種操作，需要的朋友可以參考下
2023-04-04
R語言繪制line?plot線圖示例詳解
這篇文章主要為大家介紹了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語言常用命令集合
這篇文章主要介紹了R語言常用命令集合，文中整理的都是常用的一些命令，有感興趣的可以多學習下
2021-03-03
R語言實現(xiàn)將factor轉(zhuǎn)換成numeric方法
這篇文章主要介紹了R語言實現(xiàn)將factor轉(zhuǎn)換成numeric方法，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧
2021-03-03
R語言 data.frame()命令的使用說明
這篇文章主要介紹了R語言 data.frame()命令的使用說明，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧
2021-04-04
R語言判斷語句的使用詳解
這篇文章主要介紹了R語言判斷語句的使用詳解，幫助大家更好的理解和學習使用R語言，感興趣的朋友可以了解下
2021-03-03
R語言ComplexHeatmap繪制復雜熱圖heatmap
ComplexHeatmap可以繪制很復雜的熱圖，能滿足日常以及文章所需，本文介紹一下文章常見的復雜熱圖的繪制方式，含代碼，文章比較長，建議耐心看一下，也許就有你需要的熱圖
2022-06-06