cart-elc

Source code for CART-ELC
git clone git://git.laack.co/cart-elc.git
Log | Files | Refs | README | LICENSE

zerosized.cpp (3734B)

      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 #include "main.h"
     11 
     12 
     13 template<typename MatrixType> void zeroReduction(const MatrixType& m) {
     14   // Reductions that must hold for zero sized objects
     15   VERIFY(m.all());
     16   VERIFY(!m.any());
     17   VERIFY(m.prod()==1);
     18   VERIFY(m.sum()==0);
     19   VERIFY(m.norm()==0);
     20   VERIFY(m.squaredNorm()==0);
     21   VERIFY(m.count()==0);
     22   VERIFY(m.allFinite());
     23   VERIFY(!m.hasNaN());
     24   VERIFY_RAISES_ASSERT( m.minCoeff() );
     25   VERIFY_RAISES_ASSERT( m.maxCoeff() );
     26   Index i,j;
     27   VERIFY_RAISES_ASSERT( m.minCoeff(&i,&j) );
     28   VERIFY_RAISES_ASSERT( m.maxCoeff(&i,&j) );
     29   VERIFY_RAISES_ASSERT( m.reshaped().minCoeff(&i) );
     30   VERIFY_RAISES_ASSERT( m.reshaped().maxCoeff(&i) );
     31 }
     32 
     33 
     34 template<typename MatrixType> void zeroSizedMatrix()
     35 {
     36   MatrixType t1;
     37   typedef typename MatrixType::Scalar Scalar;
     38 
     39   if (MatrixType::SizeAtCompileTime == Dynamic || MatrixType::SizeAtCompileTime == 0)
     40   {
     41     zeroReduction(t1);
     42     if (MatrixType::RowsAtCompileTime == Dynamic)
     43       VERIFY(t1.rows() == 0);
     44     if (MatrixType::ColsAtCompileTime == Dynamic)
     45       VERIFY(t1.cols() == 0);
     46 
     47     if (MatrixType::RowsAtCompileTime == Dynamic && MatrixType::ColsAtCompileTime == Dynamic)
     48     {
     49 
     50       MatrixType t2(0, 0), t3(t1);
     51       VERIFY(t2.rows() == 0);
     52       VERIFY(t2.cols() == 0);
     53 
     54       zeroReduction(t2);
     55       VERIFY(t1==t2);
     56     }
     57   }
     58 
     59   if(MatrixType::MaxColsAtCompileTime!=0 && MatrixType::MaxRowsAtCompileTime!=0)
     60   {
     61     Index rows = MatrixType::RowsAtCompileTime==Dynamic ? internal::random<Index>(1,10) : Index(MatrixType::RowsAtCompileTime);
     62     Index cols = MatrixType::ColsAtCompileTime==Dynamic ? internal::random<Index>(1,10) : Index(MatrixType::ColsAtCompileTime);
     63     MatrixType m(rows,cols);
     64     zeroReduction(m.template block<0,MatrixType::ColsAtCompileTime>(0,0,0,cols));
     65     zeroReduction(m.template block<MatrixType::RowsAtCompileTime,0>(0,0,rows,0));
     66     zeroReduction(m.template block<0,1>(0,0));
     67     zeroReduction(m.template block<1,0>(0,0));
     68     Matrix<Scalar,Dynamic,Dynamic> prod = m.template block<MatrixType::RowsAtCompileTime,0>(0,0,rows,0) * m.template block<0,MatrixType::ColsAtCompileTime>(0,0,0,cols);
     69     VERIFY(prod.rows()==rows && prod.cols()==cols);
     70     VERIFY(prod.isZero());
     71     prod = m.template block<1,0>(0,0) * m.template block<0,1>(0,0);
     72     VERIFY(prod.size()==1);
     73     VERIFY(prod.isZero());
     74   }
     75 }
     76 
     77 template<typename VectorType> void zeroSizedVector()
     78 {
     79   VectorType t1;
     80 
     81   if (VectorType::SizeAtCompileTime == Dynamic || VectorType::SizeAtCompileTime==0)
     82   {
     83     zeroReduction(t1);
     84     VERIFY(t1.size() == 0);
     85     VectorType t2(DenseIndex(0)); // DenseIndex disambiguates with 0-the-null-pointer (error with gcc 4.4 and MSVC8)
     86     VERIFY(t2.size() == 0);
     87     zeroReduction(t2);
     88 
     89     VERIFY(t1==t2);
     90   }
     91 }
     92 
     93 EIGEN_DECLARE_TEST(zerosized)
     94 {
     95   zeroSizedMatrix<Matrix2d>();
     96   zeroSizedMatrix<Matrix3i>();
     97   zeroSizedMatrix<Matrix<float, 2, Dynamic> >();
     98   zeroSizedMatrix<MatrixXf>();
     99   zeroSizedMatrix<Matrix<float, 0, 0> >();
    100   zeroSizedMatrix<Matrix<float, Dynamic, 0, 0, 0, 0> >();
    101   zeroSizedMatrix<Matrix<float, 0, Dynamic, 0, 0, 0> >();
    102   zeroSizedMatrix<Matrix<float, Dynamic, Dynamic, 0, 0, 0> >();
    103   zeroSizedMatrix<Matrix<float, 0, 4> >();
    104   zeroSizedMatrix<Matrix<float, 4, 0> >();
    105 
    106   zeroSizedVector<Vector2d>();
    107   zeroSizedVector<Vector3i>();
    108   zeroSizedVector<VectorXf>();
    109   zeroSizedVector<Matrix<float, 0, 1> >();
    110   zeroSizedVector<Matrix<float, 1, 0> >();
    111 }