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

Fix ForceNonZeroDiag for complexes

Browse Source
This commit is contained in:
Gael Guennebaud 2015-08-04 16:12:16 +02:00
parent cbce0e3b12
commit b986c147cd
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
test/sparse.h
View File

@ -77,8 +77,11 @@ initSparse(double density,
Scalar v = (internal::random<double>(0,1) < density) ? internal::random<Scalar>() : Scalar(0); Scalar v = (internal::random<double>(0,1) < density) ? internal::random<Scalar>() : Scalar(0);
if ((flags&ForceNonZeroDiag) && (i==j)) if ((flags&ForceNonZeroDiag) && (i==j))
{ {
// FIXME: the following is too conservative
v = internal::random<Scalar>()*Scalar(3.); v = internal::random<Scalar>()*Scalar(3.);
v = v*v + Scalar(5.); v = v*v;
if(numext::real(v)>0) v += Scalar(5);
else v -= Scalar(5);
} }
if ((flags & MakeLowerTriangular) && aj>ai) if ((flags & MakeLowerTriangular) && aj>ai)
v = Scalar(0); v = Scalar(0);