$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
make/h_mat.cc
1 // Copyright (C) 2009, 2012 EPITA Research and Development Laboratory
2 // (LRDE)
3 //
4 // This file is part of Olena.
5 //
6 // Olena is free software: you can redistribute it and/or modify it under
7 // the terms of the GNU General Public License as published by the Free
8 // Software Foundation, version 2 of the License.
9 //
10 // Olena is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // General Public License for more details.
14 //
15 // You should have received a copy of the GNU General Public License
16 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
17 //
18 // As a special exception, you may use this file as part of a free
19 // software project without restriction. Specifically, if other files
20 // instantiate templates or use macros or inline functions from this
21 // file, or you compile this file and link it with other files to produce
22 // an executable, this file does not by itself cause the resulting
23 // executable to be covered by the GNU General Public License. This
24 // exception does not however invalidate any other reasons why the
25 // executable file might be covered by the GNU General Public License.
26 
27 #include <mln/make/h_mat.hh>
28 #include <mln/math/abs.hh>
29 
30 
31 static const int result[4][4] = { { 2, 3, 4, 0 },
32  { 5, 6, 7, 0 },
33  { 8, 9, 10, 0 },
34  { 0, 0, 0, 1 } };
35 
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 },
39  { 0, 0, 0, 1 } };
40 
41 
42 template <typename T>
43 bool about_equal(const T& f, const T& q)
44 {
45  return mln::math::abs(q - f) <= 0.000001;
46 }
47 
48 
49 int main()
50 {
51  using namespace mln;
52 
53  int vals[9] = { 2, 3, 4,
54  5, 6, 7,
55  8, 9, 10 };
57 
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]);
61 
62 
63  // Checking creation from quaternions.
64  {
66  v[0] = 1;
67  v[1] = 2;
68  v[2] = 3;
69  v.normalize();
70 
71  algebra::quat q(0.92388, 0.186238, 0.310397, 0.124159);
72  q.set_unit();
73 
75 
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]));
79  }
80 }