devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

bug #1509: fix computeInverseWithCheck for complexes

Browse Source 
(grafted from 2f833b1c64
)
This commit is contained in:
Gael Guennebaud
2018-04-04 15:47:46 +02:00
parent 90cd199d4b
commit a2a2c3c865
2 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
Eigen/src/LU/InverseImpl.h
View File

@@ -404,7 +404,7 @@ inline void MatrixBase<Derived>::computeInverseWithCheck(
const RealScalar& absDeterminantThreshold const RealScalar& absDeterminantThreshold
) const ) const
{ {
RealScalar determinant; Scalar determinant;
// i'd love to put some static assertions there, but SFINAE means that they have no effect... // i'd love to put some static assertions there, but SFINAE means that they have no effect...
eigen_assert(rows() == cols()); eigen_assert(rows() == cols());
computeInverseAndDetWithCheck(inverse,determinant,invertible,absDeterminantThreshold); computeInverseAndDetWithCheck(inverse,determinant,invertible,absDeterminantThreshold);

4
test/inverse.cpp
View File

@@ -47,7 +47,7 @@ template<typename MatrixType> void inverse(const MatrixType& m)
//computeInverseAndDetWithCheck tests //computeInverseAndDetWithCheck tests
//First: an invertible matrix //First: an invertible matrix
bool invertible; bool invertible;
RealScalar det; Scalar det;
m2.setZero(); m2.setZero();
m1.computeInverseAndDetWithCheck(m2, det, invertible); m1.computeInverseAndDetWithCheck(m2, det, invertible);
@@ -113,5 +113,7 @@ void test_inverse()
CALL_SUBTEST_7( inverse(Matrix4d()) ); CALL_SUBTEST_7( inverse(Matrix4d()) );
CALL_SUBTEST_7( inverse(Matrix<double,4,4,DontAlign>()) ); CALL_SUBTEST_7( inverse(Matrix<double,4,4,DontAlign>()) );
CALL_SUBTEST_8( inverse(Matrix4cd()) );
} }
} }