DiagonalMatrix: release-quality documentation

BandMatrix: rename toDense() ---> toDenseMatrix() for consistency
This commit is contained in:
Benoit Jacob 2009-11-16 15:25:58 -05:00
parent e8d0dbf82e
commit eb6df28c6c
3 changed files with 35 additions and 18 deletions

View File

@ -171,7 +171,7 @@ class BandMatrix : public AnyMatrixBase<BandMatrix<_Scalar,Rows,Cols,Supers,Subs
return Block<DataType,1,Dynamic>(m_data, supers()-i, std::max(0,i), 1, diagonalLength(i)); return Block<DataType,1,Dynamic>(m_data, supers()-i, std::max(0,i), 1, diagonalLength(i));
} }
PlainMatrixType toDense() const PlainMatrixType toDenseMatrix() const
{ {
PlainMatrixType res(rows(),cols()); PlainMatrixType res(rows(),cols());
res.setZero(); res.setZero();

View File

@ -26,6 +26,7 @@
#ifndef EIGEN_DIAGONALMATRIX_H #ifndef EIGEN_DIAGONALMATRIX_H
#define EIGEN_DIAGONALMATRIX_H #define EIGEN_DIAGONALMATRIX_H
#ifndef EIGEN_PARSED_BY_DOXYGEN
template<typename Derived> template<typename Derived>
class DiagonalBase : public AnyMatrixBase<Derived> class DiagonalBase : public AnyMatrixBase<Derived>
{ {
@ -44,10 +45,8 @@ class DiagonalBase : public AnyMatrixBase<Derived>
typedef Matrix<Scalar, RowsAtCompileTime, ColsAtCompileTime, 0, MaxRowsAtCompileTime, MaxColsAtCompileTime> DenseMatrixType; typedef Matrix<Scalar, RowsAtCompileTime, ColsAtCompileTime, 0, MaxRowsAtCompileTime, MaxColsAtCompileTime> DenseMatrixType;
#ifndef EIGEN_PARSED_BY_DOXYGEN
inline const Derived& derived() const { return *static_cast<const Derived*>(this); } inline const Derived& derived() const { return *static_cast<const Derived*>(this); }
inline Derived& derived() { return *static_cast<Derived*>(this); } inline Derived& derived() { return *static_cast<Derived*>(this); }
#endif // not EIGEN_PARSED_BY_DOXYGEN
DenseMatrixType toDenseMatrix() const { return derived(); } DenseMatrixType toDenseMatrix() const { return derived(); }
template<typename DenseDerived> template<typename DenseDerived>
@ -77,16 +76,18 @@ void DiagonalBase<Derived>::evalTo(MatrixBase<DenseDerived> &other) const
other.setZero(); other.setZero();
other.diagonal() = diagonal(); other.diagonal() = diagonal();
} }
#endif
/** \class DiagonalMatrix /** \class DiagonalMatrix
* \nonstableyet
* *
* \brief Represents a diagonal matrix with its storage * \brief Represents a diagonal matrix with its storage
* *
* \param _Scalar the type of coefficients * \param _Scalar the type of coefficients
* \param _Size the dimension of the matrix, or Dynamic * \param SizeAtCompileTime the dimension of the matrix, or Dynamic
* \param MaxSizeAtCompileTime the dimension of the matrix, or Dynamic. This parameter is optional and defaults
* to SizeAtCompileTime. Most of the time, you do not need to specify it.
* *
* \sa class Matrix * \sa class DiagonalWrapper
*/ */
template<typename _Scalar, int SizeAtCompileTime, int MaxSizeAtCompileTime> template<typename _Scalar, int SizeAtCompileTime, int MaxSizeAtCompileTime>
struct ei_traits<DiagonalMatrix<_Scalar,SizeAtCompileTime,MaxSizeAtCompileTime> > struct ei_traits<DiagonalMatrix<_Scalar,SizeAtCompileTime,MaxSizeAtCompileTime> >
@ -100,10 +101,11 @@ class DiagonalMatrix
: public DiagonalBase<DiagonalMatrix<_Scalar,SizeAtCompileTime,MaxSizeAtCompileTime> > : public DiagonalBase<DiagonalMatrix<_Scalar,SizeAtCompileTime,MaxSizeAtCompileTime> >
{ {
public: public:
#ifndef EIGEN_PARSED_BY_DOXYGEN
typedef typename ei_traits<DiagonalMatrix>::DiagonalVectorType DiagonalVectorType; typedef typename ei_traits<DiagonalMatrix>::DiagonalVectorType DiagonalVectorType;
typedef const DiagonalMatrix& Nested; typedef const DiagonalMatrix& Nested;
typedef _Scalar Scalar; typedef _Scalar Scalar;
#endif
protected: protected:
@ -111,7 +113,9 @@ class DiagonalMatrix
public: public:
/** const version of diagonal(). */
inline const DiagonalVectorType& diagonal() const { return m_diagonal; } inline const DiagonalVectorType& diagonal() const { return m_diagonal; }
/** \returns a reference to the stored vector of diagonal coefficients. */
inline DiagonalVectorType& diagonal() { return m_diagonal; } inline DiagonalVectorType& diagonal() { return m_diagonal; }
/** Default constructor without initialization */ /** Default constructor without initialization */
@ -120,23 +124,27 @@ class DiagonalMatrix
/** Constructs a diagonal matrix with given dimension */ /** Constructs a diagonal matrix with given dimension */
inline DiagonalMatrix(int dim) : m_diagonal(dim) {} inline DiagonalMatrix(int dim) : m_diagonal(dim) {}
/** 2D only */ /** 2D constructor. */
inline DiagonalMatrix(const Scalar& x, const Scalar& y) : m_diagonal(x,y) {} inline DiagonalMatrix(const Scalar& x, const Scalar& y) : m_diagonal(x,y) {}
/** 3D only */ /** 3D constructor. */
inline DiagonalMatrix(const Scalar& x, const Scalar& y, const Scalar& z) : m_diagonal(x,y,z) {} inline DiagonalMatrix(const Scalar& x, const Scalar& y, const Scalar& z) : m_diagonal(x,y,z) {}
/** Copy constructor. */
template<typename OtherDerived> template<typename OtherDerived>
inline DiagonalMatrix(const DiagonalBase<OtherDerived>& other) : m_diagonal(other.diagonal()) {} inline DiagonalMatrix(const DiagonalBase<OtherDerived>& other) : m_diagonal(other.diagonal()) {}
#ifndef EIGEN_PARSED_BY_DOXYGEN
/** copy constructor. prevent a default copy constructor from hiding the other templated constructor */ /** copy constructor. prevent a default copy constructor from hiding the other templated constructor */
inline DiagonalMatrix(const DiagonalMatrix& other) : m_diagonal(other.diagonal()) {} inline DiagonalMatrix(const DiagonalMatrix& other) : m_diagonal(other.diagonal()) {}
#endif
/** generic constructor from expression of the diagonal coefficients */ /** generic constructor from expression of the diagonal coefficients */
template<typename OtherDerived> template<typename OtherDerived>
explicit inline DiagonalMatrix(const MatrixBase<OtherDerived>& other) : m_diagonal(other) explicit inline DiagonalMatrix(const MatrixBase<OtherDerived>& other) : m_diagonal(other)
{} {}
/** Copy operator. */
template<typename OtherDerived> template<typename OtherDerived>
DiagonalMatrix& operator=(const DiagonalBase<OtherDerived>& other) DiagonalMatrix& operator=(const DiagonalBase<OtherDerived>& other)
{ {
@ -144,6 +152,7 @@ class DiagonalMatrix
return *this; return *this;
} }
#ifndef EIGEN_PARSED_BY_DOXYGEN
/** This is a special case of the templated operator=. Its purpose is to /** This is a special case of the templated operator=. Its purpose is to
* prevent a default operator= from hiding the templated operator=. * prevent a default operator= from hiding the templated operator=.
*/ */
@ -152,23 +161,28 @@ class DiagonalMatrix
m_diagonal = other.m_diagonal(); m_diagonal = other.m_diagonal();
return *this; return *this;
} }
#endif
/** Resizes to given size. */
inline void resize(int size) { m_diagonal.resize(size); } inline void resize(int size) { m_diagonal.resize(size); }
/** Sets all coefficients to zero. */
inline void setZero() { m_diagonal.setZero(); } inline void setZero() { m_diagonal.setZero(); }
/** Resizes and sets all coefficients to zero. */
inline void setZero(int size) { m_diagonal.setZero(size); } inline void setZero(int size) { m_diagonal.setZero(size); }
/** Sets this matrix to be the identity matrix of the current size. */
inline void setIdentity() { m_diagonal.setOnes(); } inline void setIdentity() { m_diagonal.setOnes(); }
/** Sets this matrix to be the identity matrix of the given size. */
inline void setIdentity(int size) { m_diagonal.setOnes(size); } inline void setIdentity(int size) { m_diagonal.setOnes(size); }
}; };
/** \class DiagonalWrapper /** \class DiagonalWrapper
* \nonstableyet
* *
* \brief Expression of a diagonal matrix * \brief Expression of a diagonal matrix
* *
* \param _DiagonalVectorType the type of the vector of diagonal coefficients * \param _DiagonalVectorType the type of the vector of diagonal coefficients
* *
* This class is an expression of a diagonal matrix with given vector of diagonal * This class is an expression of a diagonal matrix, but not storing its own vector of diagonal coefficients,
* coefficients. It is the return type of MatrixBase::asDiagonal() * instead wrapping an existing vector expression. It is the return type of MatrixBase::asDiagonal()
* and most of the time this is the only way that it is used. * and most of the time this is the only way that it is used.
* *
* \sa class DiagonalMatrix, class DiagonalBase, MatrixBase::asDiagonal() * \sa class DiagonalMatrix, class DiagonalBase, MatrixBase::asDiagonal()
@ -192,18 +206,22 @@ class DiagonalWrapper
: public DiagonalBase<DiagonalWrapper<_DiagonalVectorType> >, ei_no_assignment_operator : public DiagonalBase<DiagonalWrapper<_DiagonalVectorType> >, ei_no_assignment_operator
{ {
public: public:
#ifndef EIGEN_PARSED_BY_DOXYGEN
typedef _DiagonalVectorType DiagonalVectorType; typedef _DiagonalVectorType DiagonalVectorType;
typedef DiagonalWrapper Nested; typedef DiagonalWrapper Nested;
#endif
/** Constructor from expression of diagonal coefficients to wrap. */
inline DiagonalWrapper(const DiagonalVectorType& diagonal) : m_diagonal(diagonal) {} inline DiagonalWrapper(const DiagonalVectorType& diagonal) : m_diagonal(diagonal) {}
/** \returns a const reference to the wrapped expression of diagonal coefficients. */
const DiagonalVectorType& diagonal() const { return m_diagonal; } const DiagonalVectorType& diagonal() const { return m_diagonal; }
protected: protected:
const typename DiagonalVectorType::Nested m_diagonal; const typename DiagonalVectorType::Nested m_diagonal;
}; };
/** \nonstableyet /** \returns a pseudo-expression of a diagonal matrix with *this as vector of diagonal coefficients
* \returns a pseudo-expression of a diagonal matrix with *this as vector of diagonal coefficients
* *
* \only_for_vectors * \only_for_vectors
* *
@ -219,8 +237,7 @@ MatrixBase<Derived>::asDiagonal() const
return derived(); return derived();
} }
/** \nonstableyet /** \returns true if *this is approximately equal to a diagonal matrix,
* \returns true if *this is approximately equal to a diagonal matrix,
* within the precision given by \a prec. * within the precision given by \a prec.
* *
* Example: \include MatrixBase_isDiagonal.cpp * Example: \include MatrixBase_isDiagonal.cpp

View File

@ -53,7 +53,7 @@ template<typename MatrixType> void bandmatrix(const MatrixType& _m)
dm1.diagonal(-i).setConstant(-static_cast<RealScalar>(i)); dm1.diagonal(-i).setConstant(-static_cast<RealScalar>(i));
} }
//std::cerr << m.m_data << "\n\n" << m.toDense() << "\n\n" << dm1 << "\n\n\n\n"; //std::cerr << m.m_data << "\n\n" << m.toDense() << "\n\n" << dm1 << "\n\n\n\n";
VERIFY_IS_APPROX(dm1,m.toDense()); VERIFY_IS_APPROX(dm1,m.toDenseMatrix());
for (int i=0; i<cols; ++i) for (int i=0; i<cols; ++i)
{ {
@ -68,7 +68,7 @@ template<typename MatrixType> void bandmatrix(const MatrixType& _m)
dm1.block(subs+1,0,rows-subs-1-b,rows-subs-1-b).template triangularView<LowerTriangular>().setZero(); dm1.block(subs+1,0,rows-subs-1-b,rows-subs-1-b).template triangularView<LowerTriangular>().setZero();
if(b>0) dm1.block(d+subs,0,b,cols).setZero(); if(b>0) dm1.block(d+subs,0,b,cols).setZero();
//std::cerr << m.m_data << "\n\n" << m.toDense() << "\n\n" << dm1 << "\n\n"; //std::cerr << m.m_data << "\n\n" << m.toDense() << "\n\n" << dm1 << "\n\n";
VERIFY_IS_APPROX(dm1,m.toDense()); VERIFY_IS_APPROX(dm1,m.toDenseMatrix());
} }