关于Eigen中矩阵混淆的问题(讨论贴)

可以先看一下这两篇

(12条消息) Eigen有哪些需要注意的操作_eigen noalias_Leo-Peng的博客-CSDN博客

数值计算库Eigen:别名(Aliasing) - 知乎 (zhihu.com)

        大概意思就是如果不同矩阵共享内存空间,会导致矩阵计算出错。尤其是矩阵乘法中,等号左右出现同一个变量时,matA  = matA*matB。这时候就需要使用eval(),相当于两个矩阵相乘后再赋值给临时变量,再将临时变量给matA。

        又说方阵相乘不会出现混淆,但是如果出现转置和逆,就有可能出现。

(非方阵相乘一定会出现混淆吗?????什么时候会出现混淆???意思每次矩阵乘法,如果等号两边有同一个变量时就加上eval()函数吗??)

        如果绝对不会出现内存重叠,如matC=matA*matB,则可以matC使用noalias()提高速度,不创建中间变量。(why????什么时候会自动创建中间临时变量???怎么提高速度的??意思是,我以后都使用自己创建的临时变量,都可以用noalias()???)

//如:
state_.noalias() += kalmanGainSubset * innovationSubset;

Eigen::MatrixXd gainResidual = identity_;
gainResidual.noalias() -= kalmanGainSubset * stateToMeasurementSubset;
estimateErrorCovariance_ = gainResidual * estimateErrorCovariance_ gainResidual.transpose();
estimateErrorCovariance_.noalias() += kalmanGainSubset *
                                            measurementCovarianceSubset *
                                            kalmanGainSubset.transpose();

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

XiangrongZ

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值