Commit ef378198 authored by Davis King's avatar Davis King

Relaxed the test for float matrices a little.

--HG--
extra : convert_revision : svn%3Afdd8eb12-d10e-0410-9acb-85c331704f74/trunk%403481
parent d7c051ea
......@@ -69,7 +69,7 @@ namespace
void test_eigenvalue ( const matrix_type& m)
{
typedef typename matrix_type::type type;
const type eps = max(abs(m))*sqrt(std::numeric_limits<type>::epsilon());
const type eps = 10*max(abs(m))*sqrt(std::numeric_limits<type>::epsilon());
dlog << LDEBUG << "test_eigenvalue(): " << m.nr() << " x " << m.nc() << " eps: " << eps;
print_spinner();
......@@ -106,7 +106,7 @@ namespace
DLIB_TEST(diagm(diag(CD)) == CD);
// verify that these things are actually eigenvalues and eigenvectors of m
DLIB_TEST(max(abs(m*V - V*D)) < eps);
DLIB_TEST_MSG(max(abs(m*V - V*D)) < eps, max(abs(m*V - V*D)) << " " << eps);
DLIB_TEST(max(norm(CM*CV - CV*CD)) < eps);
// if m is a symmetric matrix
......@@ -163,7 +163,7 @@ namespace
test_eigenvalue(10*randm<float>(3,3));
test_eigenvalue(10*randm<float>(4,4));
test_eigenvalue(10*randm<float>(15,15));
test_eigenvalue(10*randm<float>(150,150));
test_eigenvalue(10*randm<float>(50,50));
test_eigenvalue(10*randm<float,1,1>());
test_eigenvalue(10*randm<float,2,2>());
......@@ -174,7 +174,7 @@ namespace
test_eigenvalue(10*symm(randm<float>(3,3)));
test_eigenvalue(10*symm(randm<float>(4,4)));
test_eigenvalue(10*symm(randm<float>(15,15)));
test_eigenvalue(10*symm(randm<float>(150,150)));
test_eigenvalue(10*symm(randm<float>(50,50)));
}
// ----------------------------------------------------------------------------------------
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment