GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-12 11:49:02 +08:00
Code Issues Packages Projects Releases Wiki Activity

ComplexSchur: fix bug introduced in my previous commit.

Browse Source 
The value of c is actually used a few lines later.
This commit is contained in:
Jitse Niesen 2010-02-26 09:47:17 +00:00
parent 90e4a605ef
commit d86f5339b2
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Eigen/src/Eigenvalues/ComplexSchur.h
View File

@ -218,10 +218,12 @@ void ComplexSchur<MatrixType>::compute(const MatrixType& matrix, bool skipU)
sf = t.cwiseAbs().sum(); sf = t.cwiseAbs().sum();
t /= sf; // the normalization by sf is to avoid under/overflow t /= sf; // the normalization by sf is to avoid under/overflow
b = t.coeff(0,0) + t.coeff(1,1); b = t.coeff(0,1) * t.coeff(1,0);
c = t.coeff(0,0) - t.coeff(1,1); c = t.coeff(0,0) - t.coeff(1,1);
disc = ei_sqrt(c*c + RealScalar(4)*t.coeff(0,1)*t.coeff(1,0)); disc = ei_sqrt(c*c + RealScalar(4)*b);
c = t.coeff(0,0) * t.coeff(1,1) - b;
b = t.coeff(0,0) + t.coeff(1,1);
r1 = (b+disc)/RealScalar(2); r1 = (b+disc)/RealScalar(2);
r2 = (b-disc)/RealScalar(2); r2 = (b-disc)/RealScalar(2);