1 // This file is part of Eigen, a lightweight C++ template library 2 // for linear algebra. 3 // 4 // Copyright (C) 2011 Benoit Jacob <jacob.benoit.1 (at) gmail.com> 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 #if defined EIGEN_TEST_PART_1 || defined EIGEN_TEST_PART_2 || defined EIGEN_TEST_PART_3 || defined EIGEN_TEST_PART_4 11 #define EIGEN_DONT_ALIGN 12 #elif defined EIGEN_TEST_PART_5 || defined EIGEN_TEST_PART_6 || defined EIGEN_TEST_PART_7 || defined EIGEN_TEST_PART_8 13 #define EIGEN_DONT_ALIGN_STATICALLY 14 #endif 15 16 #include "main.h" 17 #include <Eigen/Dense> 18 19 template<typename MatrixType> 20 void dontalign(const MatrixType& m) 21 { 22 typedef typename MatrixType::Index Index; 23 typedef typename MatrixType::Scalar Scalar; 24 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType; 25 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> SquareMatrixType; 26 27 Index rows = m.rows(); 28 Index cols = m.cols(); 29 30 MatrixType a = MatrixType::Random(rows,cols); 31 SquareMatrixType square = SquareMatrixType::Random(rows,rows); 32 VectorType v = VectorType::Random(rows); 33 34 VERIFY_IS_APPROX(v, square * square.colPivHouseholderQr().solve(v)); 35 square = square.inverse().eval(); 36 a = square * a; 37 square = square*square; 38 v = square * v; 39 v = a.adjoint() * v; 40 VERIFY(square.determinant() != Scalar(0)); 41 42 // bug 219: MapAligned() was giving an assert with EIGEN_DONT_ALIGN, because Map Flags were miscomputed 43 Scalar* array = internal::aligned_new<Scalar>(rows); 44 v = VectorType::MapAligned(array, rows); 45 internal::aligned_delete(array, rows); 46 } 47 48 void test_dontalign() 49 { 50 #if defined EIGEN_TEST_PART_1 || defined EIGEN_TEST_PART_5 51 dontalign(Matrix3d()); 52 dontalign(Matrix4f()); 53 #elif defined EIGEN_TEST_PART_2 || defined EIGEN_TEST_PART_6 54 dontalign(Matrix3cd()); 55 dontalign(Matrix4cf()); 56 #elif defined EIGEN_TEST_PART_3 || defined EIGEN_TEST_PART_7 57 dontalign(Matrix<float, 32, 32>()); 58 dontalign(Matrix<std::complex<float>, 32, 32>()); 59 #elif defined EIGEN_TEST_PART_4 || defined EIGEN_TEST_PART_8 60 dontalign(MatrixXd(32, 32)); 61 dontalign(MatrixXcf(32, 32)); 62 #endif 63 } 64
