Division by Zero
Submitted by Kacper Pluta
Assigned to Nobody
Link to original bugzilla bug (#1679)
Version: 3.3 (current stable)
typename ComplexSchur<MatrixType>::ComplexScalar ComplexSchur<MatrixType>::computeShift(Index iu, Index iter)
if(numext::norm1(eival1) > numext::norm1(eival2))
eival2 = det / eival1;
eival1 = det / eival2;
but it can happen that eival1 and eival2 are zeros and this leads to division by zero.
I found the bug while computing eigenvalues and placing in my code:
feenableexcept ( FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW );
Do not hesitate to ask me for more information.