1 // This file is part of Eigen, a lightweight C++ template library 2 // for linear algebra. 3 // 4 // Copyright (C) 2009 Gael Guennebaud <g.gael (at) free.fr> 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 #include "main.h" 11 #include <unsupported/Eigen/AlignedVector3> 12 13 template<typename Scalar> 14 void alignedvector3() 15 { 16 Scalar s1 = internal::random<Scalar>(); 17 Scalar s2 = internal::random<Scalar>(); 18 typedef Matrix<Scalar,3,1> RefType; 19 typedef Matrix<Scalar,3,3> Mat33; 20 typedef AlignedVector3<Scalar> FastType; 21 RefType r1(RefType::Random()), r2(RefType::Random()), r3(RefType::Random()), 22 r4(RefType::Random()), r5(RefType::Random()), r6(RefType::Random()); 23 FastType f1(r1), f2(r2), f3(r3), f4(r4), f5(r5), f6(r6); 24 Mat33 m1(Mat33::Random()); 25 26 VERIFY_IS_APPROX(f1,r1); 27 VERIFY_IS_APPROX(f4,r4); 28 29 VERIFY_IS_APPROX(f4+f1,r4+r1); 30 VERIFY_IS_APPROX(f4-f1,r4-r1); 31 VERIFY_IS_APPROX(f4+f1-f2,r4+r1-r2); 32 VERIFY_IS_APPROX(f4+=f3,r4+=r3); 33 VERIFY_IS_APPROX(f4-=f5,r4-=r5); 34 VERIFY_IS_APPROX(f4-=f5+f1,r4-=r5+r1); 35 VERIFY_IS_APPROX(f5+f1-s1*f2,r5+r1-s1*r2); 36 VERIFY_IS_APPROX(f5+f1/s2-s1*f2,r5+r1/s2-s1*r2); 37 38 VERIFY_IS_APPROX(m1*f4,m1*r4); 39 VERIFY_IS_APPROX(f4.transpose()*m1,r4.transpose()*m1); 40 41 VERIFY_IS_APPROX(f2.dot(f3),r2.dot(r3)); 42 VERIFY_IS_APPROX(f2.cross(f3),r2.cross(r3)); 43 VERIFY_IS_APPROX(f2.norm(),r2.norm()); 44 45 VERIFY_IS_APPROX(f2.normalized(),r2.normalized()); 46 47 VERIFY_IS_APPROX((f2+f1).normalized(),(r2+r1).normalized()); 48 49 f2.normalize(); 50 r2.normalize(); 51 VERIFY_IS_APPROX(f2,r2); 52 } 53 54 void test_alignedvector3() 55 { 56 for(int i = 0; i < g_repeat; i++) { 57 CALL_SUBTEST( alignedvector3<float>() ); 58 } 59 } 60
