|  | // 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>>()); | 
|  | } |