您好,登錄后才能下訂單哦!
本篇內(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í)!
免責聲明:本站發(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)容。