mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-03-31 19:00:35 +08:00
Code Issues Packages Projects Releases Wiki Activity

restrict identity() to square matrices; small change helping g++ optimize.

Browse Source
This commit is contained in:
Benoit Jacob 2007-10-14 10:01:25 +00:00
parent 6c8f159635
commit 3f97918760
3 changed files with 10 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/Core/Identity.h
View File

@ -36,30 +36,31 @@ template<typename MatrixType> class Identity
static const int RowsAtCompileTime = MatrixType::RowsAtCompileTime,
ColsAtCompileTime = MatrixType::ColsAtCompileTime;
Identity(int rows, int cols) : m_rows(rows), m_cols(cols)
Identity(int rows) : m_rows(rows)
{
assert(rows > 0 && cols > 0);
assert(rows > 0);
assert(RowsAtCompileTime == ColsAtCompileTime);
}
private:
Identity& _ref() { return *this; }
const Identity& _constRef() const { return *this; }
int _rows() const { return m_rows; }
int _cols() const { return m_cols; }
int _cols() const { return m_rows; }
Scalar _read(int row, int col) const
{
return static_cast<Scalar>(row == col);
return row == col ? static_cast<Scalar>(1) : static_cast<Scalar>(0);
}
protected:
int m_rows, m_cols;
int m_rows;
};
template<typename Scalar, typename Derived>
Identity<Derived> Object<Scalar, Derived>::identity(int rows, int cols)
Identity<Derived> Object<Scalar, Derived>::identity(int rows)
{
return Identity<Derived>(rows, cols);
return Identity<Derived>(rows);
}
#endif // EI_IDENTITY_H

2
src/Core/Object.h
View File

@ -116,7 +116,7 @@ template<typename Scalar, typename Derived> class Object
static Zero<Derived>
zero(int rows = RowsAtCompileTime, int cols = ColsAtCompileTime);
static Identity<Derived>
identity(int rows = RowsAtCompileTime, int cols = ColsAtCompileTime);
identity(int rows = RowsAtCompileTime);
template<typename OtherDerived>
bool isApprox(

2
test/basicstuff.cpp
View File

@ -37,7 +37,7 @@ template<typename MatrixType> void basicStuff(const MatrixType& m)
m3,
mzero = MatrixType::zero(rows, cols),
identity = Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime>
::identity(rows, rows),
::identity(rows),
square = Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime>
::random(rows, rows);
VectorType v1 = VectorType::random(rows),