X-Git-Url: http://git.tdb.fi/?p=libs%2Fmath.git;a=blobdiff_plain;f=tests%2Fmatrix.cpp;fp=tests%2Fmatrix.cpp;h=1b01cbc94e5ab5cd940446d4a85c8748896765ae;hp=5b303dbe8cd99f706d4ce4889d766ee71f291657;hb=157764dbef39bf2f516fe98ef5f8efcbf2d3782a;hpb=9ab061034720c1c0db4e8bff5dae4191ccd2a26f

diff --git a/tests/matrix.cpp b/tests/matrix.cpp
index 5b303db..1b01cbc 100644
--- a/tests/matrix.cpp
+++ b/tests/matrix.cpp
@@ -1,7 +1,9 @@
+#include <limits>
 #include <msp/linal/matrix.h>
 #include <msp/linal/squarematrix.h>
 #include <msp/test/test.h>
 
+using namespace std;
 using namespace Msp;
 
 class MatrixTests: public Test::RegisteredTest<MatrixTests>
@@ -42,9 +44,10 @@ void MatrixTests::multiply()
 template<typename T, unsigned N>
 bool is_identity(const LinAl::Matrix<T, N, N> &m)
 {
+	static const T limit = numeric_limits<T>::epsilon()*4;
 	for(unsigned i=0; i<N; ++i)
 		for(unsigned j=0; j<N; ++j)
-			if(abs(T(i==j)-m(i, j))>1e-10)
+			if(abs(T(i==j)-m(i, j))>limit)
 				return false;
 	return true;
 }