溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》

R語言RcppEigen計算點乘與矩陣乘法連乘算法錯誤解決方法是什么

發(fā)布時間:2021-11-06 13:32:22 來源:億速云 閱讀:266 作者:iii 欄目:開發(fā)技術(shù)

本篇內(nèi)容主要講解“R語言RcppEigen計算點乘與矩陣乘法連乘算法錯誤解決方法是什么”,感興趣的朋友不妨來看看。本文介紹的方法操作簡單快捷,實用性強。下面就讓小編來帶大家學(xué)習(xí)“R語言RcppEigen計算點乘與矩陣乘法連乘算法錯誤解決方法是什么”吧!

計算點乘與矩陣乘法連乘計算錯誤

當我們想將 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;

這樣就可以解決之前的計算錯誤問題了。

到此,相信大家對“R語言RcppEigen計算點乘與矩陣乘法連乘算法錯誤解決方法是什么”有了更深的了解,不妨來實際操作一番吧!這里是億速云網(wǎng)站,更多相關(guān)內(nèi)容可以進入相關(guān)頻道進行查詢,關(guān)注我們,繼續(xù)學(xué)習(xí)!

向AI問一下細節(jié)

免責聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI