| // This file is part of Eigen, a lightweight C++ template library |
| // for linear algebra. |
| // |
| // Copyright (C) 2019 David Tellenbach <david.tellenbach@tellnotes.org> |
| // |
| // 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/. |
| |
| #if defined(__GNUC__) && __GNUC__ >= 10 |
| // GCC 10+ has a bug for unsigned char that thinks we're writing past the |
| // end of an array when compiled with -O3. This warning is not triggered for |
| // any other types, nor for other compilers, nor for other optimization levels. |
| #pragma GCC diagnostic ignored "-Wstringop-overflow" |
| #endif |
| |
| #include "main.h" |
| |
| template <typename Scalar, bool is_integer = NumTraits<Scalar>::IsInteger> |
| struct TestMethodDispatching { |
| static void run() {} |
| }; |
| |
| template <typename Scalar> |
| struct TestMethodDispatching<Scalar, 1> { |
| static void run() { |
| { |
| Matrix<Scalar, Dynamic, Dynamic> m{3, 4}; |
| Array<Scalar, Dynamic, Dynamic> a{3, 4}; |
| VERIFY(m.rows() == 3); |
| VERIFY(m.cols() == 4); |
| VERIFY(a.rows() == 3); |
| VERIFY(a.cols() == 4); |
| } |
| { |
| Matrix<Scalar, 1, 2> m{3, 4}; |
| Array<Scalar, 1, 2> a{3, 4}; |
| VERIFY(m(0) == 3); |
| VERIFY(m(1) == 4); |
| VERIFY(a(0) == 3); |
| VERIFY(a(1) == 4); |
| } |
| { |
| Matrix<Scalar, 2, 1> m{3, 4}; |
| Array<Scalar, 2, 1> a{3, 4}; |
| VERIFY(m(0) == 3); |
| VERIFY(m(1) == 4); |
| VERIFY(a(0) == 3); |
| VERIFY(a(1) == 4); |
| } |
| } |
| }; |
| |
| template <typename Vec4, typename Vec5> |
| void fixedsizeVariadicVectorConstruction2() { |
| { |
| Vec4 ref = Vec4::Random(); |
| Vec4 v{ref[0], ref[1], ref[2], ref[3]}; |
| VERIFY_IS_APPROX(v, ref); |
| VERIFY_IS_APPROX(v, (Vec4(ref[0], ref[1], ref[2], ref[3]))); |
| VERIFY_IS_APPROX(v, (Vec4({ref[0], ref[1], ref[2], ref[3]}))); |
| |
| Vec4 v2 = {ref[0], ref[1], ref[2], ref[3]}; |
| VERIFY_IS_APPROX(v2, ref); |
| } |
| { |
| Vec5 ref = Vec5::Random(); |
| Vec5 v{ref[0], ref[1], ref[2], ref[3], ref[4]}; |
| VERIFY_IS_APPROX(v, ref); |
| VERIFY_IS_APPROX(v, (Vec5(ref[0], ref[1], ref[2], ref[3], ref[4]))); |
| VERIFY_IS_APPROX(v, (Vec5({ref[0], ref[1], ref[2], ref[3], ref[4]}))); |
| |
| Vec5 v2 = {ref[0], ref[1], ref[2], ref[3], ref[4]}; |
| VERIFY_IS_APPROX(v2, ref); |
| } |
| } |
| |
| #define CHECK_MIXSCALAR_V5_APPROX(V, A0, A1, A2, A3, A4) \ |
| { \ |
| VERIFY_IS_APPROX(V[0], Scalar(A0)); \ |
| VERIFY_IS_APPROX(V[1], Scalar(A1)); \ |
| VERIFY_IS_APPROX(V[2], Scalar(A2)); \ |
| VERIFY_IS_APPROX(V[3], Scalar(A3)); \ |
| VERIFY_IS_APPROX(V[4], Scalar(A4)); \ |
| } |
| |
| #define CHECK_MIXSCALAR_V5(VEC5, A0, A1, A2, A3, A4) \ |
| { \ |
| typedef VEC5::Scalar Scalar; \ |
| VEC5 v = {A0, A1, A2, A3, A4}; \ |
| CHECK_MIXSCALAR_V5_APPROX(v, A0, A1, A2, A3, A4); \ |
| } |
| |
| template <int> |
| void fixedsizeVariadicVectorConstruction3() { |
| typedef Matrix<double, 5, 1> Vec5; |
| typedef Array<float, 5, 1> Arr5; |
| CHECK_MIXSCALAR_V5(Vec5, 1, 2., -3, 4.121, 5.53252); |
| CHECK_MIXSCALAR_V5(Arr5, 1, 2., 3.12f, 4.121, 5.53252); |
| } |
| |
| template <typename Scalar> |
| void fixedsizeVariadicVectorConstruction() { |
| CALL_SUBTEST((fixedsizeVariadicVectorConstruction2<Matrix<Scalar, 4, 1>, Matrix<Scalar, 5, 1>>())); |
| CALL_SUBTEST((fixedsizeVariadicVectorConstruction2<Matrix<Scalar, 1, 4>, Matrix<Scalar, 1, 5>>())); |
| CALL_SUBTEST((fixedsizeVariadicVectorConstruction2<Array<Scalar, 4, 1>, Array<Scalar, 5, 1>>())); |
| CALL_SUBTEST((fixedsizeVariadicVectorConstruction2<Array<Scalar, 1, 4>, Array<Scalar, 1, 5>>())); |
| } |
| |
| template <typename Scalar> |
| void initializerListVectorConstruction() { |
| Scalar raw[4]; |
| for (int k = 0; k < 4; ++k) { |
| raw[k] = internal::random<Scalar>(); |
| } |
| { |
| Matrix<Scalar, 4, 1> m{{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}; |
| Array<Scalar, 4, 1> a{{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}; |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(m(k) == raw[k]); |
| } |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(a(k) == raw[k]); |
| } |
| VERIFY_IS_EQUAL(m, (Matrix<Scalar, 4, 1>({{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}))); |
| VERIFY((a == (Array<Scalar, 4, 1>({{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}))).all()); |
| } |
| { |
| Matrix<Scalar, 1, 4> m{{raw[0], raw[1], raw[2], raw[3]}}; |
| Array<Scalar, 1, 4> a{{raw[0], raw[1], raw[2], raw[3]}}; |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(m(k) == raw[k]); |
| } |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(a(k) == raw[k]); |
| } |
| VERIFY_IS_EQUAL(m, (Matrix<Scalar, 1, 4>({{raw[0], raw[1], raw[2], raw[3]}}))); |
| VERIFY((a == (Array<Scalar, 1, 4>({{raw[0], raw[1], raw[2], raw[3]}}))).all()); |
| } |
| { |
| Matrix<Scalar, 4, Dynamic> m{{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}; |
| Array<Scalar, 4, Dynamic> a{{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}; |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(m(k) == raw[k]); |
| } |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(a(k) == raw[k]); |
| } |
| VERIFY_IS_EQUAL(m, (Matrix<Scalar, 4, Dynamic>({{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}))); |
| VERIFY((a == (Array<Scalar, 4, Dynamic>({{raw[0]}, {raw[1]}, {raw[2]}, {raw[3]}}))).all()); |
| } |
| { |
| Matrix<Scalar, Dynamic, 4> m{{raw[0], raw[1], raw[2], raw[3]}}; |
| Array<Scalar, Dynamic, 4> a{{raw[0], raw[1], raw[2], raw[3]}}; |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(m(k) == raw[k]); |
| } |
| for (int k = 0; k < 4; ++k) { |
| VERIFY(a(k) == raw[k]); |
| } |
| VERIFY_IS_EQUAL(m, (Matrix<Scalar, Dynamic, 4>({{raw[0], raw[1], raw[2], raw[3]}}))); |
| VERIFY((a == (Array<Scalar, Dynamic, 4>({{raw[0], raw[1], raw[2], raw[3]}}))).all()); |
| } |
| } |
| |
| template <typename Scalar> |
| void initializerListMatrixConstruction() { |
| const Index RowsAtCompileTime = 5; |
| const Index ColsAtCompileTime = 4; |
| const Index SizeAtCompileTime = RowsAtCompileTime * ColsAtCompileTime; |
| |
| Scalar raw[SizeAtCompileTime]; |
| for (int i = 0; i < SizeAtCompileTime; ++i) { |
| raw[i] = internal::random<Scalar>(); |
| } |
| { |
| Matrix<Scalar, Dynamic, Dynamic> m{}; |
| VERIFY(m.cols() == 0); |
| VERIFY(m.rows() == 0); |
| VERIFY_IS_EQUAL(m, (Matrix<Scalar, Dynamic, Dynamic>())); |
| } |
| { |
| Matrix<Scalar, 5, 4> m{{raw[0], raw[1], raw[2], raw[3]}, |
| {raw[4], raw[5], raw[6], raw[7]}, |
| {raw[8], raw[9], raw[10], raw[11]}, |
| {raw[12], raw[13], raw[14], raw[15]}, |
| {raw[16], raw[17], raw[18], raw[19]}}; |
| |
| Matrix<Scalar, 5, 4> m2; |
| m2 << raw[0], raw[1], raw[2], raw[3], raw[4], raw[5], raw[6], raw[7], raw[8], raw[9], raw[10], raw[11], raw[12], |
| raw[13], raw[14], raw[15], raw[16], raw[17], raw[18], raw[19]; |
| |
| int k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| VERIFY(m(i, j) == raw[k]); |
| ++k; |
| } |
| } |
| VERIFY_IS_EQUAL(m, m2); |
| } |
| { |
| Matrix<Scalar, Dynamic, Dynamic> m{{raw[0], raw[1], raw[2], raw[3]}, |
| {raw[4], raw[5], raw[6], raw[7]}, |
| {raw[8], raw[9], raw[10], raw[11]}, |
| {raw[12], raw[13], raw[14], raw[15]}, |
| {raw[16], raw[17], raw[18], raw[19]}}; |
| |
| VERIFY(m.cols() == 4); |
| VERIFY(m.rows() == 5); |
| int k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| VERIFY(m(i, j) == raw[k]); |
| ++k; |
| } |
| } |
| |
| Matrix<Scalar, Dynamic, Dynamic> m2(RowsAtCompileTime, ColsAtCompileTime); |
| k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| m2(i, j) = raw[k]; |
| ++k; |
| } |
| } |
| VERIFY_IS_EQUAL(m, m2); |
| } |
| } |
| |
| template <typename Scalar> |
| void initializerListArrayConstruction() { |
| const Index RowsAtCompileTime = 5; |
| const Index ColsAtCompileTime = 4; |
| const Index SizeAtCompileTime = RowsAtCompileTime * ColsAtCompileTime; |
| |
| Scalar raw[SizeAtCompileTime]; |
| for (int i = 0; i < SizeAtCompileTime; ++i) { |
| raw[i] = internal::random<Scalar>(); |
| } |
| { |
| Array<Scalar, Dynamic, Dynamic> a{}; |
| VERIFY(a.cols() == 0); |
| VERIFY(a.rows() == 0); |
| } |
| { |
| Array<Scalar, 5, 4> m{{raw[0], raw[1], raw[2], raw[3]}, |
| {raw[4], raw[5], raw[6], raw[7]}, |
| {raw[8], raw[9], raw[10], raw[11]}, |
| {raw[12], raw[13], raw[14], raw[15]}, |
| {raw[16], raw[17], raw[18], raw[19]}}; |
| |
| Array<Scalar, 5, 4> m2; |
| m2 << raw[0], raw[1], raw[2], raw[3], raw[4], raw[5], raw[6], raw[7], raw[8], raw[9], raw[10], raw[11], raw[12], |
| raw[13], raw[14], raw[15], raw[16], raw[17], raw[18], raw[19]; |
| |
| int k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| VERIFY(m(i, j) == raw[k]); |
| ++k; |
| } |
| } |
| VERIFY_IS_APPROX(m, m2); |
| } |
| { |
| Array<Scalar, Dynamic, Dynamic> m{{raw[0], raw[1], raw[2], raw[3]}, |
| {raw[4], raw[5], raw[6], raw[7]}, |
| {raw[8], raw[9], raw[10], raw[11]}, |
| {raw[12], raw[13], raw[14], raw[15]}, |
| {raw[16], raw[17], raw[18], raw[19]}}; |
| |
| VERIFY(m.cols() == 4); |
| VERIFY(m.rows() == 5); |
| int k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| VERIFY(m(i, j) == raw[k]); |
| ++k; |
| } |
| } |
| |
| Array<Scalar, Dynamic, Dynamic> m2(RowsAtCompileTime, ColsAtCompileTime); |
| k = 0; |
| for (int i = 0; i < RowsAtCompileTime; ++i) { |
| for (int j = 0; j < ColsAtCompileTime; ++j) { |
| m2(i, j) = raw[k]; |
| ++k; |
| } |
| } |
| VERIFY_IS_APPROX(m, m2); |
| } |
| } |
| |
| template <typename Scalar> |
| void dynamicVectorConstruction() { |
| const Index size = 4; |
| Scalar raw[size]; |
| for (int i = 0; i < size; ++i) { |
| raw[i] = internal::random<Scalar>(); |
| } |
| |
| typedef Matrix<Scalar, Dynamic, 1> VectorX; |
| |
| { |
| VectorX v{{raw[0], raw[1], raw[2], raw[3]}}; |
| for (int i = 0; i < size; ++i) { |
| VERIFY(v(i) == raw[i]); |
| } |
| VERIFY(v.rows() == size); |
| VERIFY(v.cols() == 1); |
| VERIFY_IS_EQUAL(v, (VectorX{{raw[0], raw[1], raw[2], raw[3]}})); |
| } |
| } |
| |
| EIGEN_DECLARE_TEST(initializer_list_construction) { |
| CALL_SUBTEST_1(initializerListVectorConstruction<unsigned char>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<float>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<double>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<int>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<long int>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<std::ptrdiff_t>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<std::complex<double>>()); |
| CALL_SUBTEST_1(initializerListVectorConstruction<std::complex<float>>()); |
| |
| CALL_SUBTEST_2(initializerListMatrixConstruction<unsigned char>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<float>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<double>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<int>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<long int>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<std::ptrdiff_t>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<std::complex<double>>()); |
| CALL_SUBTEST_2(initializerListMatrixConstruction<std::complex<float>>()); |
| |
| CALL_SUBTEST_3(initializerListArrayConstruction<unsigned char>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<float>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<double>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<int>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<long int>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<std::ptrdiff_t>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<std::complex<double>>()); |
| CALL_SUBTEST_3(initializerListArrayConstruction<std::complex<float>>()); |
| |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<unsigned char>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<float>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<double>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<int>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<long int>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<std::ptrdiff_t>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<std::complex<double>>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction<std::complex<float>>()); |
| CALL_SUBTEST_4(fixedsizeVariadicVectorConstruction3<0>()); |
| |
| CALL_SUBTEST_5(TestMethodDispatching<int>::run()); |
| CALL_SUBTEST_5(TestMethodDispatching<long int>::run()); |
| |
| CALL_SUBTEST_6(dynamicVectorConstruction<unsigned char>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<float>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<double>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<int>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<long int>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<std::ptrdiff_t>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<std::complex<double>>()); |
| CALL_SUBTEST_6(dynamicVectorConstruction<std::complex<float>>()); |
| } |