如何比較向量大約在特徵？

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

如何比較向量大約在特徵？

回答

相關問題