首页> R语言教程

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

withpy 2021-11-12

简介这篇文章主要为大家介绍了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计算点乘与矩阵乘法的资料请关注本站其它相关文章!

相关文章

  • docker镜像完全卸载的操作步骤

    这篇文章主要介绍了docker镜像完全卸载的操作步骤,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧

  • C语言实现扫雷游戏详解

    这篇文章主要为大家详细介绍了C语言实现扫雷游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下

  • idea2020.2卡死在reading maven projects

    这篇文章主要介绍了idea2020.2卡死在reading maven projects,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

  • Python filter过滤器原理及实例应用

    这篇文章主要介绍了Python filter过滤器原理及实例应用,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下

  • Python实现的服务器示例小结【单进程、多进程、多线程、非阻塞式】

    这篇文章主要介绍了Python实现的服务器,结合实例形式分析了Python实现单进程、多进程、多线程、非阻塞式服务器的相关操作技巧,需要的朋友可以参考下

  • java实现读取txt文件并以在每行以空格取数据

    今天小编就为大家分享一篇java实现读取txt文件并以在每行以空格取数据,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧