Fixed division by zero corner case in array unit test.

This commit is contained in:
Hauke Heibel 2012-03-09 14:04:13 +01:00
parent d7da6f63a8
commit dd9365e089

View File

@ -220,8 +220,9 @@ template<typename ArrayType> void array_real(const ArrayType& m)
VERIFY_IS_APPROX(std::pow(m3,RealScalar(0.5)), m3.sqrt());
// scalar by array division
ArrayType t1 = s1/m1;
ArrayType t2 = s1 * m1.inverse();
const RealScalar tiny = std::sqrt(std::numeric_limits<RealScalar>::epsilon());
s1 += Scalar(tiny);
m1 += ArrayType::Constant(rows,cols,Scalar(tiny));
VERIFY_IS_APPROX(s1/m1, s1 * m1.inverse());
}