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语言 title()函数的参数用法说明
这篇文章主要介绍了R语言 title()函数的参数用法说明，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2021-04-04
R语言检验多重共线性的操作
这篇文章主要介绍了R语言检验多重共线性的操作，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2021-04-04
R语言的数据输入深入讲解
这篇文章主要介绍了R语言的数据输入深入讲解，代码实例讲解的很清晰，有需要的同学可以学习下
2021-03-03
R语言数据重塑知识点总结
在本篇文章里小编给大家整理的是一篇关于R语言数据重塑的相关知识点内容，有兴趣的朋友们可以跟着学习下。
2021-03-03
使用R语言填补缺失值的方法
这篇文章主要介绍了使用R语言填补缺失值的方法,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2021-03-03
R语言的历史介绍
在本篇内容里小编给大家介绍了关于R语言的历史内容，有兴趣学习的朋友可以阅读下。
2021-03-03
用R语言绘制函数曲线图
这篇文章主要介绍了如何用R语言绘制函数曲线图，帮助大家更好的理解和学习使用R语言，感兴趣的朋友可以了解下
2021-03-03
R语言及RStudio下载与安装方法
这篇文章主要介绍了R、RStudio下载与安装方法,本文通过图文图文并茂的方法给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2021-03-03
R语言学习代码格式一键美化
这篇文章主要为大家讲解介绍了R语言学习中代码格式一键美化的实现过程步骤，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步
2021-11-11
R studio 批量注释的快捷方式
这篇文章主要介绍了R studio 批量注释的快捷方式，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2021-04-04