概述
计算点乘与矩阵乘法连乘计算错误
当我们想将 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计算点乘与矩阵乘法的资料请关注靠谱客其它相关文章！
最后
以上就是大方鱼为你收集整理的R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决的全部内容，希望文章能够帮你解决R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复