R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决
计算点乘与矩阵乘法连乘计算错误
当我们想将 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;
这样就可以解决之前的计算错误问题了。
【R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决】以上就是RcppEigen计算点乘与矩阵乘法连乘计算错误解决的详细内容,更多关于RcppEigen计算点乘与矩阵乘法的资料请关注脚本之家其它相关文章!
推荐阅读
- 【生信技能树】R语言练习题|【生信技能树】R语言练习题 - 中级
- 一起来学习C语言的字符串转换函数
- C语言字符函数中的isalnum()和iscntrl()你都知道吗
- C语言浮点函数中的modf和fmod详解
- C语言中的时间函数clock()和time()你都了解吗
- 使用协程爬取网页,计算网页数据大小
- C语言学习|第十一届蓝桥杯省赛 大学B组 C/C++ 第一场
- ACSL|ACSL 美国计算机科学联赛 2016-2017 R4 摩天大楼-Skyscraper 题解
- C语言解方程的根和判断是否是闰年
- C语言的版本比较