From 0715d49908e59487e040d49e03319b0a9093b969 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 13 Feb 2014 15:49:54 +0100
Subject: [PATCH] Fix stable_norm unit test for complexes

---
 test/stable_norm.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/test/stable_norm.cpp b/test/stable_norm.cpp
index 364170acd..549f91fbf 100644
--- a/test/stable_norm.cpp
+++ b/test/stable_norm.cpp
@@ -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));