1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra. Eigen itself is part of the KDE project. 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2008 Gael Guennebaud <g.gael@free.fr> 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2008 Benoit Jacob <jacob.benoit.1@gmail.com> 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define EIGEN_WORK_AROUND_QT_BUG_CALLING_WRONG_OPERATOR_NEW_FIXED_IN_QT_4_5 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "main.h" 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <Eigen/Geometry> 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <Eigen/QtAlignedMalloc> 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <QtCore/QVector> 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename MatrixType> 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid check_qtvector_matrix(const MatrixType& m) 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int rows = m.rows(); 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int cols = m.cols(); 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType x = MatrixType::Random(rows,cols), y = MatrixType::Random(rows,cols); 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QVector<MatrixType> v(10, MatrixType(rows,cols)), w(20, y); 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < 20; i++) 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[i], y); 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[5] = x; 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath w[6] = v[5]; 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[6], v[5]); 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v = w; 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < 20; i++) 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[i], v[i]); 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.resize(21); 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[20] = x; 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[20], x); 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.fill(y,22); 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[21], y); 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(x); 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[22], x); 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY((size_t)&(v[22]) == (size_t)&(v[21]) + sizeof(MatrixType)); 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // do a lot of push_back such that the vector gets internally resized 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // (with memory reallocation) 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType* ref = &w[0]; 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i=0; i<30 || ((ref==&w[0]) && i<300); ++i) 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(w[i%w.size()]); 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i=23; i<v.size(); ++i) 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY(v[i]==w[(i-23)%w.size()]); 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename TransformType> 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid check_qtvector_transform(const TransformType&) 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename TransformType::MatrixType MatrixType; 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath TransformType x(MatrixType::Random()), y(MatrixType::Random()); 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QVector<TransformType> v(10), w(20, y); 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[5] = x; 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath w[6] = v[5]; 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[6], v[5]); 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v = w; 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < 20; i++) 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[i], v[i]); 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.resize(21); 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[20] = x; 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[20], x); 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.fill(y,22); 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[21], y); 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(x); 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[22], x); 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY((size_t)&(v[22]) == (size_t)&(v[21]) + sizeof(TransformType)); 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // do a lot of push_back such that the vector gets internally resized 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // (with memory reallocation) 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath TransformType* ref = &w[0]; 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i=0; i<30 || ((ref==&w[0]) && i<300); ++i) 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(w[i%w.size()]); 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(unsigned int i=23; int(i)<v.size(); ++i) 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY(v[i].matrix()==w[(i-23)%w.size()].matrix()); 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename QuaternionType> 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid check_qtvector_quaternion(const QuaternionType&) 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef typename QuaternionType::Coefficients Coefficients; 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QuaternionType x(Coefficients::Random()), y(Coefficients::Random()); 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QVector<QuaternionType> v(10), w(20, y); 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[5] = x; 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath w[6] = v[5]; 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[6], v[5]); 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v = w; 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < 20; i++) 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(w[i], v[i]); 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.resize(21); 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v[20] = x; 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[20], x); 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.fill(y,22); 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[21], y); 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(x); 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY_IS_APPROX(v[22], x); 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY((size_t)&(v[22]) == (size_t)&(v[21]) + sizeof(QuaternionType)); 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // do a lot of push_back such that the vector gets internally resized 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // (with memory reallocation) 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QuaternionType* ref = &w[0]; 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i=0; i<30 || ((ref==&w[0]) && i<300); ++i) 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath v.push_back(w[i%w.size()]); 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(unsigned int i=23; int(i)<v.size(); ++i) 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VERIFY(v[i].coeffs()==w[(i-23)%w.size()].coeffs()); 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid test_eigen2_qtvector() 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // some non vectorizable fixed sizes 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_1(check_qtvector_matrix(Vector2f())); 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_1(check_qtvector_matrix(Matrix3f())); 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_1(check_qtvector_matrix(Matrix3d())); 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // some vectorizable fixed sizes 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_2(check_qtvector_matrix(Matrix2f())); 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_2(check_qtvector_matrix(Vector4f())); 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_2(check_qtvector_matrix(Matrix4f())); 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_2(check_qtvector_matrix(Matrix4d())); 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // some dynamic sizes 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_3(check_qtvector_matrix(MatrixXd(1,1))); 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_3(check_qtvector_matrix(VectorXd(20))); 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_3(check_qtvector_matrix(RowVectorXf(20))); 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_3(check_qtvector_matrix(MatrixXcf(10,10))); 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // some Transform 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_4(check_qtvector_transform(Transform2f())); 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_4(check_qtvector_transform(Transform3f())); 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_4(check_qtvector_transform(Transform3d())); 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath //CALL_SUBTEST_4(check_qtvector_transform(Transform4d())); 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // some Quaternion 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_5(check_qtvector_quaternion(Quaternionf())); 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath CALL_SUBTEST_5(check_qtvector_quaternion(Quaternionf())); 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 159