Design #1279
Tidy up code for matrix determinant
Description
I recently encountered an unexpected error while trying to compute the determinant of a 2x2 matrix over ZZ/(841)
where 841 = 29^2
.
DetDirect
has an incorrect assertion that the ring is a PolyRingDetDirect
probably should not have been called anyway- why wasn't
det2x2
called?
Review the code, and tidy it up!
Related issues
History
#1 Updated by John Abbott almost 5 years ago
- Related to Bug #1280: Determinant & Inverse of matrix over non integral domain added
#2 Updated by John Abbott over 4 years ago
- Target version changed from CoCoALib-0.99700 to CoCoALib-0.99800
#3 Updated by John Abbott over 3 years ago
- Status changed from New to Resolved
- Assignee set to John Abbott
- % Done changed from 0 to 70
Main dispatch fn is ConstMatrixViewBase::myDet
around matrix.C:221
The code has already been considerably cleaned up, but still performs poorly for some inputs.
#4 Updated by John Abbott over 3 years ago
- Description updated (diff)
- Status changed from Resolved to Closed
- % Done changed from 70 to 100
- Estimated time set to 2.66 h
The current version is acceptable; it could be made better, but there are plenty of other more important issues (at the moment).
Closing.