$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
quat.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 <iostream>
27 
28 #include <mln/core/contract.hh>
29 
30 #include <mln/algebra/quat.hh>
31 #include <mln/norm/l2.hh>
32 
33 
34 int main()
35 {
36  using namespace mln;
37 
38  algebra::quat q1(3.f, 4.f, 1.6f, 0.5f);
39  algebra::quat q2(1.2f, make::vec(3, 6, 4));
40  algebra::quat q3(make::vec(1.3, 3., -6., 4.));
41 
42  std::cout << q1 << std::endl;
43  std::cout << q2 << std::endl;
44  std::cout << q3 << std::endl;
45 
46  std::cout << q1.s() << std::endl;
47 
48  q1.s() = 2.6f;
49  std::cout << q1 << std::endl;
50 
51  std::cout << q1.v() << std::endl;
52 
53  q2.v() = make::vec(1.4, 5.9, 3.1);
54  std::cout << q2 << std::endl;
55 
56  std::cout << q2 * q3 << std::endl;
57 
58  mln_assertion(! q3.is_unit());
59  q3.set_unit();
60  std::cout << q3 << std::endl;
61  mln_assertion(q3.is_unit());
62 
63  std::cout << q2.conj() << std::endl;
64  std::cout << q2.inv() << std::endl;
65  std::cout << norm::l2(q2.to_vec()) << ' '<< norm::l2(q2.inv().to_vec())
66  << std::endl;
67  std::cout << q2.inv().inv() << std::endl;
68 
69 }