From 335d3bcf054cd0d1328cc12fbc72bb8463aef859 Mon Sep 17 00:00:00 2001 From: Benoit Jacob Date: Fri, 9 Jan 2009 23:26:45 +0000 Subject: Based on code + help from Alex Stapleton: *Add Eigen/StdVector header. Including it #includes and "Core" and generates a partial specialization of std::vector for T=Eigen::Matrix<...> that will work even with vectorizable fixed-size Eigen types (working around a design issue in the c++ STL) *Add unit-test CCMAIL: alex.stapleton@gmail.com --- test/stdvector.cpp | 64 ++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 64 insertions(+) create mode 100644 test/stdvector.cpp (limited to 'test/stdvector.cpp') diff --git a/test/stdvector.cpp b/test/stdvector.cpp new file mode 100644 index 000000000..4e6aa9e3a --- /dev/null +++ b/test/stdvector.cpp @@ -0,0 +1,64 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. Eigen itself is part of the KDE project. +// +// Copyright (C) 2008 Gael Guennebaud +// +// Eigen is free software; you can redistribute it and/or +// modify it under the terms of the GNU Lesser General Public +// License as published by the Free Software Foundation; either +// version 3 of the License, or (at your option) any later version. +// +// Alternatively, you can redistribute it and/or +// modify it under the terms of the GNU General Public License as +// published by the Free Software Foundation; either version 2 of +// the License, or (at your option) any later version. +// +// Eigen is distributed in the hope that it will be useful, but WITHOUT ANY +// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS +// FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License or the +// GNU General Public License for more details. +// +// You should have received a copy of the GNU Lesser General Public +// License and a copy of the GNU General Public License along with +// Eigen. If not, see . + +#include "main.h" +#include + +template +void check_stdvector_fixedsize() +{ + MatrixType x = MatrixType::Random(), y = MatrixType::Random(); + std::vector v(10), w(20, y); + v[5] = x; + w[6] = v[5]; + VERIFY_IS_APPROX(w[6], v[5]); + v = w; + for(int i = 0; i < 20; i++) + { + VERIFY_IS_APPROX(w[i], v[i]); + } + v.resize(21); + v[20] = x; + VERIFY_IS_APPROX(v[20], x); + v.resize(22,y); + VERIFY_IS_APPROX(v[21], y); + v.push_back(x); + VERIFY_IS_APPROX(v[22], x); +} + + +void test_stdvector() +{ + // some non vectorizable fixed sizes + CALL_SUBTEST(check_stdvector_fixedsize()); + CALL_SUBTEST(check_stdvector_fixedsize()); + CALL_SUBTEST(check_stdvector_fixedsize()); + + // some vectorizable fixed sizes + CALL_SUBTEST(check_stdvector_fixedsize()); + CALL_SUBTEST(check_stdvector_fixedsize()); + CALL_SUBTEST(check_stdvector_fixedsize()); + CALL_SUBTEST(check_stdvector_fixedsize()); + +} -- cgit v1.2.3