如何比较向量大约在特征？

Eigen中是否有函数用于比较使用相对和绝对容差的矢量（矩阵）（又名numpy.allclose？如果其中一个矢量非常接近零，则标准isApprox失败。如何比较向量大约在特征？

2013-02-24 DikobrAz

请参阅http://eigen.tuxfamily.org/dox-2.0/TutorialCore.html – user1929959 2013-02-24 11:58:39

没有内置的函数实现numpy.allclose，但如果这真的是你需要的，你可以轻松地自己写一个函数。但是，我宁愿建议使用的isMuchSmallerThan参考值：

(a-b).isMuchSmallerThan(ref)

哪里裁判是你的问题有代表性的非为零。

编辑：参考这里是一个可能的实现allclose的：

template<typename DerivedA, typename DerivedB> 
bool allclose(const Eigen::DenseBase<DerivedA>& a, 
       const Eigen::DenseBase<DerivedB>& b, 
       const typename DerivedA::RealScalar& rtol 
        = Eigen::NumTraits<typename DerivedA::RealScalar>::dummy_precision(), 
       const typename DerivedA::RealScalar& atol 
        = Eigen::NumTraits<typename DerivedA::RealScalar>::epsilon()) 
{ 
    return ((a.derived() - b.derived()).array().abs() 
      <= (atol + rtol * b.derived().array().abs())).all(); 
}

来源

2013-02-24 13:20:27 ggael

isMuchSmallerThan可用于绝对比较，即（ab）.isMuchSmallerThan（1.0，atol）等价于np.allclose（a， b，0.0，atol），所以模仿np.allclose我们必须做这样的事情：（ab）.isMuchSmallerThan（1.0，atol）|| a.isApprox（b，rtol）。我对么？ – DikobrAz 2013-02-24 13:38:38

不完全是因为isMuchStylerThan和isApprox是基于L2矩阵范数而不是基于元素的比较（无限范数）。 – ggael 2013-02-24 16:46:39

还有这不是为我工作isApprox功能。我只是使用（expect-res）.norm（）<一些小数字。

来源

2017-09-01 21:22:35 nnrales

如何比较向量大约在特征？

回答

相关问题