1 // This file is part of Eigen, a lightweight C++ template library 2 // for linear algebra. 3 // 4 // Copyright (C) 2009 Gael Guennebaud <gael.guennebaud (at) inria.fr> 5 // 6 // This Source Code Form is subject to the terms of the Mozilla 7 // Public License v. 2.0. If a copy of the MPL was not distributed 8 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9 10 #define EIGEN2_SUPPORT 11 12 #include "main.h" 13 14 template<typename MatrixType> void eigen2support(const MatrixType& m) 15 { 16 typedef typename MatrixType::Index Index; 17 typedef typename MatrixType::Scalar Scalar; 18 19 Index rows = m.rows(); 20 Index cols = m.cols(); 21 22 MatrixType m1 = MatrixType::Random(rows, cols), 23 m3(rows, cols); 24 25 Scalar s1 = internal::random<Scalar>(), 26 s2 = internal::random<Scalar>(); 27 28 // scalar addition 29 VERIFY_IS_APPROX(m1.cwise() + s1, s1 + m1.cwise()); 30 VERIFY_IS_APPROX(m1.cwise() + s1, MatrixType::Constant(rows,cols,s1) + m1); 31 VERIFY_IS_APPROX((m1*Scalar(2)).cwise() - s2, (m1+m1) - MatrixType::Constant(rows,cols,s2) ); 32 m3 = m1; 33 m3.cwise() += s2; 34 VERIFY_IS_APPROX(m3, m1.cwise() + s2); 35 m3 = m1; 36 m3.cwise() -= s1; 37 VERIFY_IS_APPROX(m3, m1.cwise() - s1); 38 39 VERIFY_IS_EQUAL((m1.corner(TopLeft,1,1)), (m1.block(0,0,1,1))); 40 VERIFY_IS_EQUAL((m1.template corner<1,1>(TopLeft)), (m1.template block<1,1>(0,0))); 41 VERIFY_IS_EQUAL((m1.col(0).start(1)), (m1.col(0).segment(0,1))); 42 VERIFY_IS_EQUAL((m1.col(0).template start<1>()), (m1.col(0).segment(0,1))); 43 VERIFY_IS_EQUAL((m1.col(0).end(1)), (m1.col(0).segment(rows-1,1))); 44 VERIFY_IS_EQUAL((m1.col(0).template end<1>()), (m1.col(0).segment(rows-1,1))); 45 46 using namespace internal; 47 VERIFY_IS_EQUAL(ei_cos(s1), cos(s1)); 48 VERIFY_IS_EQUAL(ei_real(s1), real(s1)); 49 VERIFY_IS_EQUAL(ei_abs2(s1), abs2(s1)); 50 51 m1.minor(0,0); 52 } 53 54 void test_eigen2support() 55 { 56 for(int i = 0; i < g_repeat; i++) { 57 CALL_SUBTEST_1( eigen2support(Matrix<double,1,1>()) ); 58 CALL_SUBTEST_2( eigen2support(MatrixXd(1,1)) ); 59 CALL_SUBTEST_4( eigen2support(Matrix3f()) ); 60 CALL_SUBTEST_5( eigen2support(Matrix4d()) ); 61 CALL_SUBTEST_2( eigen2support(MatrixXf(200,200)) ); 62 CALL_SUBTEST_6( eigen2support(MatrixXcd(100,100)) ); 63 } 64 } 65