Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
make/h_mat.hh
1 // Copyright (C) 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 #ifndef MLN_MAKE_H_MAT_HH
27 # define MLN_MAKE_H_MAT_HH
28 
32 
33 # include <mln/algebra/h_mat.hh>
34 # include <mln/metal/math/sqrt.hh>
35 
36 
37 namespace mln
38 {
39 
40  namespace make
41  {
42 
43 
45  /*
46  * \param[in] tab C-array of values.
47  *
48  * \pre The array dimension N has to be square (N = n * n).
49  */
50  template <typename T, unsigned N>
51  algebra::h_mat<mlc_sqrt_int(N), T> h_mat(const T (&tab)[N]);
52 
53 
54 
55 # ifndef MLN_INCLUDE_ONLY
56 
57  template <typename T, unsigned N>
58  inline
59  algebra::h_mat<mlc_sqrt_int(N), T>
60  h_mat(const T (&tab)[N])
61  {
62  enum { n = mlc_sqrt_int(N) };
63  mlc_bool(N == n * n)::check();
65  for (unsigned i = 0; i < N; ++i)
66  tmp(i / n, i % n) = tab[i];
67  return tmp;
68  }
69 
70 # endif // ! MLN_INCLUDE_ONLY
71 
72  } // end of namespace mln::make
73 
74 } // end of namespace mln
75 
76 
77 #endif // ! MLN_MAKE_H_MAT_HH