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@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
19template<typename MatrixType>
20void 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
48void 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