27 #include <mln/make/h_mat.hh>
28 #include <mln/math/abs.hh>
31 static const int result[4][4] = { { 2, 3, 4, 0 },
36 static const double result_q[4][4] = { { 0.776477, -0.113801, 0.619785, 0 },
37 { 0.345031, 0.8998, -0.267046, 0 },
38 { -0.527293, 0.4212, 0.737938, 0 },
53 int vals[9] = { 2, 3, 4,
58 for (
unsigned i = 0; i < 4; ++i)
59 for (
unsigned j = 0; j < 4; ++j)
60 mln_assertion(m(i,j) == result[i][j]);
76 for (
unsigned i = 0; i < 4; ++i)
77 for (
unsigned j = 0; j < 4; ++j)
78 mln_assertion(
about_equal(mat(i,j), result_q[i][j]));