GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-20 08:39:37 +08:00
Code Issues Packages Projects Releases Wiki Activity

Use absolute error to test sum in which cancellation may occur.

Browse Source
This commit is contained in:
Jitse Niesen 2011-02-25 08:56:37 +00:00
parent ef73265987
commit 58abf0eb98
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
test/redux.cpp
View File

@ -89,7 +89,7 @@ template<typename VectorType> void vectorRedux(const VectorType& w)
minc = std::min(minc, internal::real(v[j])); minc = std::min(minc, internal::real(v[j]));
maxc = std::max(maxc, internal::real(v[j])); maxc = std::max(maxc, internal::real(v[j]));
} }
VERIFY_IS_APPROX(s, v.head(i).sum()); VERIFY_IS_MUCH_SMALLER_THAN(internal::abs(s - v.head(i).sum()), Scalar(1));
VERIFY_IS_APPROX(p, v.head(i).prod()); VERIFY_IS_APPROX(p, v.head(i).prod());
VERIFY_IS_APPROX(minc, v.real().head(i).minCoeff()); VERIFY_IS_APPROX(minc, v.real().head(i).minCoeff());
VERIFY_IS_APPROX(maxc, v.real().head(i).maxCoeff()); VERIFY_IS_APPROX(maxc, v.real().head(i).maxCoeff());
@ -123,7 +123,7 @@ template<typename VectorType> void vectorRedux(const VectorType& w)
minc = std::min(minc, internal::real(v[j])); minc = std::min(minc, internal::real(v[j]));
maxc = std::max(maxc, internal::real(v[j])); maxc = std::max(maxc, internal::real(v[j]));
} }
VERIFY_IS_APPROX(s, v.segment(i, size-2*i).sum()); VERIFY_IS_MUCH_SMALLER_THAN(internal::abs(s - v.segment(i, size-2*i).sum()), Scalar(1));
VERIFY_IS_APPROX(p, v.segment(i, size-2*i).prod()); VERIFY_IS_APPROX(p, v.segment(i, size-2*i).prod());
VERIFY_IS_APPROX(minc, v.real().segment(i, size-2*i).minCoeff()); VERIFY_IS_APPROX(minc, v.real().segment(i, size-2*i).minCoeff());
VERIFY_IS_APPROX(maxc, v.real().segment(i, size-2*i).maxCoeff()); VERIFY_IS_APPROX(maxc, v.real().segment(i, size-2*i).maxCoeff());