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语言版本升级完全攻略基于Ubuntu
这篇文章主要为大家介绍了在Ubuntu上针对不同的版本对R语言版本升级的完全攻略，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步
2021-11-11
ComplexHeatmap绘制单个热图
这篇文章主要为大家介绍了ComplexHeatmap绘制单个热图示例详解，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2022-06-06
R语言学习VennDiagram包绘制韦恩图示例
这篇文章主要为大家介绍了R语言学习VennDiagram包绘制韦恩图示例，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2022-06-06
R语言中的vector(向量)，array(数组)使用总结
这篇文章主要介绍了R语言中的vector(向量)，array(数组)使用总结，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
2021-03-03
R语言决策基础知识点详解
在本篇文章里小编给大家整理了一篇关于R语言决策基础知识点内容，有兴趣的朋友们可以学习参考下。
2021-03-03
使用R语言绘制棒棒糖图火柴杆图教程
本篇文章为大家介绍几种利用R语言绘制棒棒糖图（火柴杆图）的方法，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2021-11-11
R语言数据结构之矩阵、数组与数据框详解
进行数据分析的第一步是先拿到数据,下面这篇文章主要给大家介绍了关于R语言数据结构之矩阵、数组与数据框的相关资料,文中通过实例代码介绍的非常详细,需要的朋友可以参考下
2022-07-07
R语言作图之density plot(密度图)的制作步骤
这篇文章主要介绍了R语言作图之density plot(密度图)的制作步骤,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2021-03-03
R语言中corrplot标题居中及eps格式输出
由于论文的需要，很多R语言绘图需要eps格式矢量图的方式进行输出，同时最近也需要绘制相关性图，因此决定写个博客作为记录，有需要的朋友可以借鉴参考下
2021-11-11
详解R语言实现前向逐步回归(前向选择模型)
本文主要介绍了详解R语言实现前向逐步回归，从实现原理开始，文中通过示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2021-08-08