12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758 |
- // This file is part of Eigen, a lightweight C++ template library
- // for linear algebra.
- //
- // Copyright (C) 2013 Gael Guennebaud <gael.guennebaud@inria.fr>
- //
- // This Source Code Form is subject to the terms of the Mozilla
- // Public License v. 2.0. If a copy of the MPL was not distributed
- // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
- #include "main.h"
- template<typename Scalar> void special_numbers()
- {
- typedef Matrix<Scalar, Dynamic,Dynamic> MatType;
- int rows = internal::random<int>(1,300);
- int cols = internal::random<int>(1,300);
-
- Scalar nan = std::numeric_limits<Scalar>::quiet_NaN();
- Scalar inf = std::numeric_limits<Scalar>::infinity();
- Scalar s1 = internal::random<Scalar>();
-
- MatType m1 = MatType::Random(rows,cols),
- mnan = MatType::Random(rows,cols),
- minf = MatType::Random(rows,cols),
- mboth = MatType::Random(rows,cols);
-
- int n = internal::random<int>(1,10);
- for(int k=0; k<n; ++k)
- {
- mnan(internal::random<int>(0,rows-1), internal::random<int>(0,cols-1)) = nan;
- minf(internal::random<int>(0,rows-1), internal::random<int>(0,cols-1)) = inf;
- }
- mboth = mnan + minf;
-
- VERIFY(!m1.hasNaN());
- VERIFY(m1.allFinite());
-
- VERIFY(mnan.hasNaN());
- VERIFY((s1*mnan).hasNaN());
- VERIFY(!minf.hasNaN());
- VERIFY(!(2*minf).hasNaN());
- VERIFY(mboth.hasNaN());
- VERIFY(mboth.array().hasNaN());
-
- VERIFY(!mnan.allFinite());
- VERIFY(!minf.allFinite());
- VERIFY(!(minf-mboth).allFinite());
- VERIFY(!mboth.allFinite());
- VERIFY(!mboth.array().allFinite());
- }
- EIGEN_DECLARE_TEST(special_numbers)
- {
- for(int i = 0; i < 10*g_repeat; i++) {
- CALL_SUBTEST_1( special_numbers<float>() );
- CALL_SUBTEST_1( special_numbers<double>() );
- }
- }
|