diff options
author | Gael Guennebaud <g.gael@free.fr> | 2009-05-03 17:19:19 +0000 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2009-05-03 17:19:19 +0000 |
commit | c25fc50d540501f5b937792c6f758d14d6317530 (patch) | |
tree | 197b0ad1b33b6bd74422e8e1af762b80b3237de6 /unsupported/test/alignedvector3.cpp | |
parent | 469b8aa380323f89c59cfd4cde4f9849aa5b9cc7 (diff) |
add an AlignedVector3 module suitable for vectorization of 3D vectors
Diffstat (limited to 'unsupported/test/alignedvector3.cpp')
-rw-r--r-- | unsupported/test/alignedvector3.cpp | 71 |
1 files changed, 71 insertions, 0 deletions
diff --git a/unsupported/test/alignedvector3.cpp b/unsupported/test/alignedvector3.cpp new file mode 100644 index 000000000..6a5a700ab --- /dev/null +++ b/unsupported/test/alignedvector3.cpp @@ -0,0 +1,71 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. Eigen itself is part of the KDE project. +// +// Copyright (C) 2009 Gael Guennebaud <g.gael@free.fr> +// +// 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 <http://www.gnu.org/licenses/>. + +#include "main.h" +#include <unsupported/Eigen/AlignedVector3> + +template<typename Scalar> +void alignedvector3() +{ + Scalar s1 = ei_random<Scalar>(); + Scalar s2 = ei_random<Scalar>(); + Scalar s3 = ei_random<Scalar>(); + typedef Matrix<Scalar,3,1> RefType; + typedef Matrix<Scalar,3,3> Mat33; + typedef AlignedVector3<Scalar> FastType; + RefType r1(RefType::Random()), r2(RefType::Random()), r3(RefType::Random()), + r4(RefType::Random()), r5(RefType::Random()), r6(RefType::Random()); + FastType f1(r1), f2(r2), f3(r3), f4(r4), f5(r5), f6(r6); + Mat33 m1(Mat33::Random()); + + VERIFY_IS_APPROX(f1,r1); + VERIFY_IS_APPROX(f4,r4); + + VERIFY_IS_APPROX(f4+f1,r4+r1); + VERIFY_IS_APPROX(f4-f1,r4-r1); + VERIFY_IS_APPROX(f4+f1-f2,r4+r1-r2); + VERIFY_IS_APPROX(f4+=f3,r4+=r3); + VERIFY_IS_APPROX(f4-=f5,r4-=r5); + VERIFY_IS_APPROX(f4-=f5+f1,r4-=r5+r1); + VERIFY_IS_APPROX(f5+f1-s1*f2,r5+r1-s1*r2); + VERIFY_IS_APPROX(f5+f1/s2-s1*f2,r5+r1/s2-s1*r2); + + VERIFY_IS_APPROX(m1*f4,m1*r4); + VERIFY_IS_APPROX(f4.transpose()*m1,r4.transpose()*m1); + + VERIFY_IS_APPROX(f2.dot(f3),r2.dot(r3)); + VERIFY_IS_APPROX(f2.cross(f3),r2.cross(r3)); + VERIFY_IS_APPROX(f2.norm(),r2.norm()); + + f2.normalize(); + r2.normalize(); + VERIFY_IS_APPROX(f2,r2); +} + +void test_alignedvector3() +{ + for(int i = 0; i < g_repeat; i++) { + CALL_SUBTEST(( alignedvector3<float>() )); + } +} |