首页 » 软件编程 » R语言

R语言RcppEigen计算点乘与矩阵乘法连乘算法错误解决

R语言 2021-12-12

计算点乘与矩阵乘法连乘计算错误

当我们想将 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计算点乘与矩阵乘法的资料请关注www.580doc.com其它相关文章!


上一篇:R语言包ggplot实现分面去掉小标题的灰色底色小技巧下一篇:Rcpp入门R代码提速方法过程
程序园_程序员的世界 Copyright © 2020- www.580doc.com. Some Rights Reserved.