diff options
author | 2012-12-20 18:08:26 +0100 | |
---|---|---|
committer | 2012-12-20 18:08:26 +0100 | |
commit | 85005ffbd12625bfe24d6e486a4a5c20633ca9db (patch) | |
tree | abf15cb27810271d818011fd6224667bcfeaa842 /test/cwiseop.cpp | |
parent | b7ea59556dbf9c73b05638a041cff1c37b664d90 (diff) |
Make sure sqrt and the likes are not compiled for integer type in cwiseop unit test.
Diffstat (limited to 'test/cwiseop.cpp')
-rw-r--r-- | test/cwiseop.cpp | 58 |
1 files changed, 38 insertions, 20 deletions
diff --git a/test/cwiseop.cpp b/test/cwiseop.cpp index ca6e4211c..24fc26d52 100644 --- a/test/cwiseop.cpp +++ b/test/cwiseop.cpp @@ -28,6 +28,42 @@ template<typename Scalar> struct AddIfNull { enum { Cost = NumTraits<Scalar>::AddCost }; }; +template<typename MatrixType> +typename Eigen::internal::enable_if<!NumTraits<typename MatrixType::Scalar>::IsInteger,typename MatrixType::Scalar>::type +cwiseops_real_only(MatrixType& m1, MatrixType& m2, MatrixType& m3, MatrixType& mones) +{ + typedef typename MatrixType::Scalar Scalar; + typedef typename NumTraits<Scalar>::Real RealScalar; + + VERIFY_IS_APPROX(m1.cwise() / m2, m1.cwise() * (m2.cwise().inverse())); + m3 = m1.cwise().abs().cwise().sqrt(); + VERIFY_IS_APPROX(m3.cwise().square(), m1.cwise().abs()); + VERIFY_IS_APPROX(m1.cwise().square().cwise().sqrt(), m1.cwise().abs()); + VERIFY_IS_APPROX(m1.cwise().abs().cwise().log().cwise().exp() , m1.cwise().abs()); + + VERIFY_IS_APPROX(m1.cwise().pow(2), m1.cwise().square()); + m3 = (m1.cwise().abs().cwise()<=RealScalar(0.01)).select(mones,m1); + VERIFY_IS_APPROX(m3.cwise().pow(-1), m3.cwise().inverse()); + m3 = m1.cwise().abs(); + VERIFY_IS_APPROX(m3.cwise().pow(RealScalar(0.5)), m3.cwise().sqrt()); + +// VERIFY_IS_APPROX(m1.cwise().tan(), m1.cwise().sin().cwise() / m1.cwise().cos()); + VERIFY_IS_APPROX(mones, m1.cwise().sin().cwise().square() + m1.cwise().cos().cwise().square()); + m3 = m1; + m3.cwise() /= m2; + VERIFY_IS_APPROX(m3, m1.cwise() / m2); + + return Scalar(0); +} + +template<typename MatrixType> +typename Eigen::internal::enable_if<NumTraits<typename MatrixType::Scalar>::IsInteger,typename MatrixType::Scalar>::type +cwiseops_real_only(MatrixType& , MatrixType& , MatrixType& , MatrixType& ) +{ + typedef typename MatrixType::Scalar Scalar; + return 0; +} + template<typename MatrixType> void cwiseops(const MatrixType& m) { typedef typename MatrixType::Index Index; @@ -101,26 +137,6 @@ template<typename MatrixType> void cwiseops(const MatrixType& m) VERIFY_IS_APPROX(m3, m1.cwise() * m2); VERIFY_IS_APPROX(mones, m2.cwise()/m2); - if(!NumTraits<Scalar>::IsInteger) - { - VERIFY_IS_APPROX(m1.cwise() / m2, m1.cwise() * (m2.cwise().inverse())); - m3 = m1.cwise().abs().cwise().sqrt(); - VERIFY_IS_APPROX(m3.cwise().square(), m1.cwise().abs()); - VERIFY_IS_APPROX(m1.cwise().square().cwise().sqrt(), m1.cwise().abs()); - VERIFY_IS_APPROX(m1.cwise().abs().cwise().log().cwise().exp() , m1.cwise().abs()); - - VERIFY_IS_APPROX(m1.cwise().pow(2), m1.cwise().square()); - m3 = (m1.cwise().abs().cwise()<=RealScalar(0.01)).select(mones,m1); - VERIFY_IS_APPROX(m3.cwise().pow(-1), m3.cwise().inverse()); - m3 = m1.cwise().abs(); - VERIFY_IS_APPROX(m3.cwise().pow(RealScalar(0.5)), m3.cwise().sqrt()); - -// VERIFY_IS_APPROX(m1.cwise().tan(), m1.cwise().sin().cwise() / m1.cwise().cos()); - VERIFY_IS_APPROX(mones, m1.cwise().sin().cwise().square() + m1.cwise().cos().cwise().square()); - m3 = m1; - m3.cwise() /= m2; - VERIFY_IS_APPROX(m3, m1.cwise() / m2); - } // check min VERIFY_IS_APPROX( m1.cwise().min(m2), m2.cwise().min(m1) ); @@ -150,6 +166,8 @@ template<typename MatrixType> void cwiseops(const MatrixType& m) VERIFY( (m1.cwise()<m1.unaryExpr(bind2nd(plus<Scalar>(), Scalar(1)))).all() ); VERIFY( !(m1.cwise()<m1.unaryExpr(bind2nd(minus<Scalar>(), Scalar(1)))).all() ); VERIFY( !(m1.cwise()>m1.unaryExpr(bind2nd(plus<Scalar>(), Scalar(1)))).any() ); + + cwiseops_real_only(m1, m2, m3, mones); } void test_cwiseop() |