ground/gcs/src/libs/eigen/test/zerosized.cpp

   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.count()==0);
  20   VERIFY(m.allFinite());
  21   VERIFY(!m.hasNaN());
  22 }
  23
  24
  25 template<typename MatrixType> void zeroSizedMatrix()
  26 {
  27   MatrixType t1;
  28   typedef typename MatrixType::Scalar Scalar;
  29
  30   if (MatrixType::SizeAtCompileTime == Dynamic || MatrixType::SizeAtCompileTime == 0)
  31   {
  32     zeroReduction(t1);
  33     if (MatrixType::RowsAtCompileTime == Dynamic)
  34       VERIFY(t1.rows() == 0);
  35     if (MatrixType::ColsAtCompileTime == Dynamic)
  36       VERIFY(t1.cols() == 0);
  37
  38     if (MatrixType::RowsAtCompileTime == Dynamic && MatrixType::ColsAtCompileTime == Dynamic)
  39     {
  40
  41       MatrixType t2(0, 0), t3(t1);
  42       VERIFY(t2.rows() == 0);
  43       VERIFY(t2.cols() == 0);
  44
  45       zeroReduction(t2);
  46       VERIFY(t1==t2);
  47     }
  48   }
  49
  50   if(MatrixType::MaxColsAtCompileTime!=0 && MatrixType::MaxRowsAtCompileTime!=0)
  51   {
  52     Index rows = MatrixType::RowsAtCompileTime==Dynamic ? internal::random<Index>(1,10) : Index(MatrixType::RowsAtCompileTime);
  53     Index cols = MatrixType::ColsAtCompileTime==Dynamic ? internal::random<Index>(1,10) : Index(MatrixType::ColsAtCompileTime);
  54     MatrixType m(rows,cols);
  55     zeroReduction(m.template block<0,MatrixType::ColsAtCompileTime>(0,0,0,cols));
  56     zeroReduction(m.template block<MatrixType::RowsAtCompileTime,0>(0,0,rows,0));
  57     zeroReduction(m.template block<0,1>(0,0));
  58     zeroReduction(m.template block<1,0>(0,0));
  59     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);
  60     VERIFY(prod.rows()==rows && prod.cols()==cols);
  61     VERIFY(prod.isZero());
  62     prod = m.template block<1,0>(0,0) * m.template block<0,1>(0,0);
  63     VERIFY(prod.size()==1);
  64     VERIFY(prod.isZero());
  65   }
  66 }
  67
  68 template<typename VectorType> void zeroSizedVector()
  69 {
  70   VectorType t1;
  71
  72   if (VectorType::SizeAtCompileTime == Dynamic || VectorType::SizeAtCompileTime==0)
  73   {
  74     zeroReduction(t1);
  75     VERIFY(t1.size() == 0);
  76     VectorType t2(DenseIndex(0)); // DenseIndex disambiguates with 0-the-null-pointer (error with gcc 4.4 and MSVC8)
  77     VERIFY(t2.size() == 0);
  78     zeroReduction(t2);
  79
  80     VERIFY(t1==t2);
  81   }
  82 }
  83
  84 void test_zerosized()
  85 {
  86   zeroSizedMatrix<Matrix2d>();
  87   zeroSizedMatrix<Matrix3i>();
  88   zeroSizedMatrix<Matrix<float, 2, Dynamic> >();
  89   zeroSizedMatrix<MatrixXf>();
  90   zeroSizedMatrix<Matrix<float, 0, 0> >();
  91   zeroSizedMatrix<Matrix<float, Dynamic, 0, 0, 0, 0> >();
  92   zeroSizedMatrix<Matrix<float, 0, Dynamic, 0, 0, 0> >();
  93   zeroSizedMatrix<Matrix<float, Dynamic, Dynamic, 0, 0, 0> >();
  94   zeroSizedMatrix<Matrix<float, 0, 4> >();
  95   zeroSizedMatrix<Matrix<float, 4, 0> >();
  96
  97   zeroSizedVector<Vector2d>();
  98   zeroSizedVector<Vector3i>();
  99   zeroSizedVector<VectorXf>();
 100   zeroSizedVector<Matrix<float, 0, 1> >();
 101   zeroSizedVector<Matrix<float, 1, 0> >();
 102 }