Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
algebra/h_mat.cc
1 // Copyright (C) 2007, 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 <mln/fun/i2v/all_to.hh>
27 #include <mln/algebra/mat.hh>
28 #include <mln/algebra/h_mat.hh>
29 
30 
31 
32 int main()
33 {
34  using namespace mln;
35 
36  algebra::mat<1,3,float> m1;
37  m1.set_all(4);
38  algebra::mat<2,2,float> m2 = literal::identity;
39 
42  algebra::h_mat<3,float> hm3(all_to(1.5));
43 
44  algebra::mat<4,4,float> m4 = hm3;
45 
46  std::cout << "m1 = " << m1 << ";" << std::endl;
47  std::cout << "m2 = " << m2 << ";" << std::endl;
48  std::cout << "m4 = " << m4 << ";" << std::endl;
49  std::cout << "hm1 = " << hm1 << ";" << std::endl;
50  std::cout << "hm2 = " << hm2 << ";" << std::endl;
51  std::cout << "hm3 = " << hm3 << ";" << std::endl;
52 
53  {
55  m = m2;
56  // FIXME: Test *many* => runs ok...
57  }
58 
59 }