diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-04-22 20:49:01 -0400 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-04-22 20:49:01 -0400 |
commit | a4f9ca44ab3de019cbe9bb1247eab47ac5e22ff3 (patch) | |
tree | 10186ed9061e26f74b74c6937c3777eba55d5d8d | |
parent | 2362eadcdd710038b43ba8bde17dc124fe0d00bb (diff) |
add minor to Eigen2Support
-rw-r--r-- | Eigen/Core | 4 | ||||
-rw-r--r-- | Eigen/Eigen2Support | 1 | ||||
-rw-r--r-- | Eigen/src/Core/MatrixBase.h | 3 | ||||
-rw-r--r-- | Eigen/src/Core/util/ForwardDeclarations.h | 1 | ||||
-rw-r--r-- | Eigen/src/Eigen2Support/Minor.h | 124 | ||||
-rw-r--r-- | test/eigen2support.cpp | 2 |
6 files changed, 135 insertions, 0 deletions
diff --git a/Eigen/Core b/Eigen/Core index f06b10064..e92ea5700 100644 --- a/Eigen/Core +++ b/Eigen/Core @@ -191,6 +191,10 @@ inline static const char *SimdInstructionSetsInUse(void) { #endif } +#ifdef EIGEN2_SUPPORT +#undef minor +#endif + // we use size_t frequently and we'll never remember to prepend it with std:: everytime just to // ensure QNX/QCC support using std::size_t; diff --git a/Eigen/Eigen2Support b/Eigen/Eigen2Support index b876cc4c3..44e862226 100644 --- a/Eigen/Eigen2Support +++ b/Eigen/Eigen2Support @@ -49,6 +49,7 @@ namespace Eigen { #include "src/Eigen2Support/TriangularSolver.h" #include "src/Eigen2Support/Block.h" #include "src/Eigen2Support/VectorBlock.h" +#include "src/Eigen2Support/Minor.h" } // namespace Eigen diff --git a/Eigen/src/Core/MatrixBase.h b/Eigen/src/Core/MatrixBase.h index 4b0ee1187..4496eccf1 100644 --- a/Eigen/src/Core/MatrixBase.h +++ b/Eigen/src/Core/MatrixBase.h @@ -409,6 +409,9 @@ template<typename Derived> class MatrixBase template<int Size> const VectorBlock<Derived,Size> start() const; template<int Size> VectorBlock<Derived,Size> end(); template<int Size> const VectorBlock<Derived,Size> end() const; + + Minor<Derived> minor(int row, int col); + const Minor<Derived> minor(int row, int col) const; #endif protected: diff --git a/Eigen/src/Core/util/ForwardDeclarations.h b/Eigen/src/Core/util/ForwardDeclarations.h index 25f808377..faee3be8e 100644 --- a/Eigen/src/Core/util/ForwardDeclarations.h +++ b/Eigen/src/Core/util/ForwardDeclarations.h @@ -190,6 +190,7 @@ struct ei_stem_function #ifdef EIGEN2_SUPPORT template<typename ExpressionType> class Cwise; +template<typename MatrixType> class Minor; #endif #endif // EIGEN_FORWARDDECLARATIONS_H diff --git a/Eigen/src/Eigen2Support/Minor.h b/Eigen/src/Eigen2Support/Minor.h new file mode 100644 index 000000000..e7e164a16 --- /dev/null +++ b/Eigen/src/Eigen2Support/Minor.h @@ -0,0 +1,124 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. +// +// Copyright (C) 2006-2009 Benoit Jacob <jacob.benoit.1@gmail.com> +// +// 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/>. + +#ifndef EIGEN_MINOR_H +#define EIGEN_MINOR_H + +/** \nonstableyet + * \class Minor + * + * \brief Expression of a minor + * + * \param MatrixType the type of the object in which we are taking a minor + * + * This class represents an expression of a minor. It is the return + * type of MatrixBase::minor() and most of the time this is the only way it + * is used. + * + * \sa MatrixBase::minor() + */ +template<typename MatrixType> +struct ei_traits<Minor<MatrixType> > + : ei_traits<MatrixType> +{ + typedef typename ei_nested<MatrixType>::type MatrixTypeNested; + typedef typename ei_unref<MatrixTypeNested>::type _MatrixTypeNested; + enum { + RowsAtCompileTime = (MatrixType::RowsAtCompileTime != Dynamic) ? + int(MatrixType::RowsAtCompileTime) - 1 : Dynamic, + ColsAtCompileTime = (MatrixType::ColsAtCompileTime != Dynamic) ? + int(MatrixType::ColsAtCompileTime) - 1 : Dynamic, + MaxRowsAtCompileTime = (MatrixType::MaxRowsAtCompileTime != Dynamic) ? + int(MatrixType::MaxRowsAtCompileTime) - 1 : Dynamic, + MaxColsAtCompileTime = (MatrixType::MaxColsAtCompileTime != Dynamic) ? + int(MatrixType::MaxColsAtCompileTime) - 1 : Dynamic, + Flags = _MatrixTypeNested::Flags & HereditaryBits, + CoeffReadCost = _MatrixTypeNested::CoeffReadCost // minor is used typically on tiny matrices, + // where loops are unrolled and the 'if' evaluates at compile time + }; +}; + +template<typename MatrixType> class Minor + : public MatrixBase<Minor<MatrixType> > +{ + public: + + typedef MatrixBase<Minor> Base; + EIGEN_DENSE_PUBLIC_INTERFACE(Minor) + + inline Minor(const MatrixType& matrix, + int row, int col) + : m_matrix(matrix), m_row(row), m_col(col) + { + ei_assert(row >= 0 && row < matrix.rows() + && col >= 0 && col < matrix.cols()); + } + + EIGEN_INHERIT_ASSIGNMENT_OPERATORS(Minor) + + inline int rows() const { return m_matrix.rows() - 1; } + inline int cols() const { return m_matrix.cols() - 1; } + + inline Scalar& coeffRef(int row, int col) + { + return m_matrix.const_cast_derived().coeffRef(row + (row >= m_row), col + (col >= m_col)); + } + + inline const Scalar coeff(int row, int col) const + { + return m_matrix.coeff(row + (row >= m_row), col + (col >= m_col)); + } + + protected: + const typename MatrixType::Nested m_matrix; + const int m_row, m_col; +}; + +/** \nonstableyet + * \return an expression of the (\a row, \a col)-minor of *this, + * i.e. an expression constructed from *this by removing the specified + * row and column. + * + * Example: \include MatrixBase_minor.cpp + * Output: \verbinclude MatrixBase_minor.out + * + * \sa class Minor + */ +template<typename Derived> +inline Minor<Derived> +MatrixBase<Derived>::minor(int row, int col) +{ + return Minor<Derived>(derived(), row, col); +} + +/** \nonstableyet + * This is the const version of minor(). */ +template<typename Derived> +inline const Minor<Derived> +MatrixBase<Derived>::minor(int row, int col) const +{ + return Minor<Derived>(derived(), row, col); +} + +#endif // EIGEN_MINOR_H diff --git a/test/eigen2support.cpp b/test/eigen2support.cpp index f8d601984..a18e7ed8a 100644 --- a/test/eigen2support.cpp +++ b/test/eigen2support.cpp @@ -57,6 +57,8 @@ template<typename MatrixType> void eigen2support(const MatrixType& m) VERIFY_IS_EQUAL((m1.col(0).template start<1>()), (m1.col(0).segment(0,1))); VERIFY_IS_EQUAL((m1.col(0).end(1)), (m1.col(0).segment(rows-1,1))); VERIFY_IS_EQUAL((m1.col(0).template end<1>()), (m1.col(0).segment(rows-1,1))); + + m1.minor(0,0); } void test_eigen2support() |