GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 04:35:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix stable_norm unit test for complexes

Browse Source 
(grafted from 0715d49908e59487e040d49e03319b0a9093b969
)
This commit is contained in:
Gael Guennebaud 2014-02-13 15:49:54 +01:00
parent 52dc1d7ffd
commit 9a09b75df3
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
test/stable_norm.cpp
View File

@ -57,12 +57,12 @@ template<typename MatrixType> void stable_norm(const MatrixType& m)
// get a non-zero random factor
Scalar factor = internal::random<Scalar>();
while(factor<RealScalar(1e-3))
while(numext::abs2(factor)<RealScalar(1e-4))
factor = internal::random<Scalar>();
Scalar big = factor * ((std::numeric_limits<RealScalar>::max)() * RealScalar(1e-4));
factor = internal::random<Scalar>();
while(factor<RealScalar(1e-3))
while(numext::abs2(factor)<RealScalar(1e-4))
factor = internal::random<Scalar>();
Scalar small = factor * ((std::numeric_limits<RealScalar>::min)() * RealScalar(1e4));