计算点乘与矩阵乘法连乘计算错误 当我们想将 R 中的连乘(如下公式所示)修改成 Rcpp 代码时, t(X)^2 %*% X 理论上我们只用在 .cpp 代码中输入下述语句即
当我们想将 R 中的连乘(如下公式所示)修改成 Rcpp 代码时,
t(X)^2 %*% X
理论上我们只用在 .cpp 代码中输入下述语句即可(默认使用了 RcppEigen
库):
X.adjoint().array().square() * X.array().square();
但实际上这样会会出现问题,原因是 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计算点乘与矩阵乘法的资料请关注编程网其它相关文章!
--结束END--
本文标题: R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决
本文链接: https://www.lsjlt.com/news/156600.html(转载时请注明来源链接)
有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341
下载Word文档到电脑,方便收藏和打印~
2024-03-01
2024-03-01
2024-02-29
2024-02-29
2024-02-29
2024-02-29
2024-02-29
2024-02-29
2024-02-29
2024-02-29
回答
回答
回答
回答
回答
回答
回答
回答
回答
回答
0