R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决
更新时间：2021年11月06日 11:29:49 作者：Kanny广小隶
这篇文章主要为大家介绍了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;
这样就可以解决之前的计算错误问题了。
以上就是RcppEigen计算点乘与矩阵乘法连乘计算错误解决的详细内容，更多关于RcppEigen计算点乘与矩阵乘法的资料请关注脚本之家其它相关文章！
相关文章
这篇文章主要介绍了基于R语言for循环的替换方案，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2021-04-04
这篇文章主要介绍了R语言中aggregate 函数的相关资料,aggregate函数是数据处理中常用到的函数，具有强大的功能，可以按照要求把数据打组聚合，然后对聚合以后的数据进行加和、求平均等各种操作，需要的朋友可以参考下2023-04-04
这篇文章主要为大家介绍了R语言绘制line plot线图的示例详解，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪2022-02-02
这篇文章主要为大家介绍了R语言绘图数据可视化Ridgeline plot山脊图画法的示例详解，有需要的朋友可以借鉴参考下，希望能够有所帮助2022-02-02
这篇文章主要介绍了R语言循环读取excel并保存为RData的操作，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2021-04-04
这篇文章主要介绍了R语言常用命令集合，文中整理的都是常用的一些命令，有感兴趣的可以多学习下2021-03-03
这篇文章主要介绍了R语言 实现将factor转换成numeric方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2021-03-03
这篇文章主要介绍了R语言 data.frame()命令的使用说明，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2021-04-04
这篇文章主要介绍了R语言判断语句的使用详解，帮助大家更好的理解和学习使用R语言，感兴趣的朋友可以了解下2021-03-03
R语言ComplexHeatmap绘制复杂热图heatmap
ComplexHeatmap可以绘制很复杂的热图，能满足日常以及文章所需，本文介绍一下文章常见的复杂热图的绘制方式，含代码，文章比较长，建议耐心看一下，也许就有你需要的热图2022-06-06
最新评论