Add static assertion for incompatible scalar types in decompositions
Submitted by seb..@..uni.lu
Assigned to Nobody
Link to original bugzilla bug (#949)
Version: 3.2
Description
Created attachment 537
Demo
For a regular square matrix with n >= 5 rows, the determinant computation returns completely wrong results. I attach a sample code.
For instance for the following matrix:
Eigen::Matrix<int,5,5> C0, C1, C2;
C0 <<
1, 2, 3, 4, 5,
3, 0, 4, 5, 6,
2, 1, 2, 3, 4,
0, 0, 0, 6, 5,
0, 0, 0, 5, 6;
C1 <<
0, 0, 3, 1, 1,
3, 1, 1, 0, 1,
1, 0, 2, 3, 3,
3, 1, 2, 3, 3,
0, 2, 1, 2, 1;
C2 <<
1, 1, 1, 1, 1,
1, 2, 3, 4, 5,
1, 3, 6, 10, 15,
1, 4, 10, 20, 35,
1, 5, 15, 35, 69;
Expected det(C0) = 99 (when C0.determinant() = -432
Expected det(C1) = -25 (when C1.determinant() = -54
Expected det(C2) = 0 (when C2.determinant() = 1296
Attachment 537, "Demo":
eigen_demo.cpp