$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
x2x/translation.cc
1 // Copyright (C) 2007, 2008, 2009 EPITA Research and Development Laboratory (LRDE)
2 //
3 // This file is part of Olena.
4 //
5 // Olena is free software: you can redistribute it and/or modify it under
6 // the terms of the GNU General Public License as published by the Free
7 // Software Foundation, version 2 of the License.
8 //
9 // Olena is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // General Public License for more details.
13 //
14 // You should have received a copy of the GNU General Public License
15 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
16 //
17 // As a special exception, you may use this file as part of a free
18 // software project without restriction. Specifically, if other files
19 // instantiate templates or use macros or inline functions from this
20 // file, or you compile this file and link it with other files to produce
21 // an executable, this file does not by itself cause the resulting
22 // executable to be covered by the GNU General Public License. This
23 // exception does not however invalidate any other reasons why the
24 // executable file might be covered by the GNU General Public License.
25 
26 #include <cmath>
27 #include <iostream>
28 #include <mln/fun/x2x/translation.hh>
29 #include <mln/fun/i2v/all_to.hh>
30 
31 
32 
33 int main()
34 {
35  using namespace mln;
36 
37  float
38  a = 2.3,
39  b = 0,
40  c = 2.9;
41 
42  algebra::vec<3,float> vec1 = make::vec(a, b, c);
44 
45  std::cout << vec1 << std::endl;
46  std::cout << tr1(vec1) << std::endl;
47  std::cout << tr1.inv()(vec1) << std::endl;
48  mln_assertion(tr1(tr1.inv()(vec1)) == vec1);
49  mln_assertion(fabs(tr1(vec1)[0] - 3.9) <= 0.125 &&
50  fabs(tr1(vec1)[1] - 1.6) <= 0.125 &&
51  fabs(tr1(vec1)[2] - 4.5) <= 0.125);
52 }