00001 // Copyright (C) 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #ifndef MLN_MAKE_H_MAT_HH 00027 # define MLN_MAKE_H_MAT_HH 00028 00032 00033 # include <mln/algebra/h_mat.hh> 00034 # include <mln/metal/math/sqrt.hh> 00035 00036 00037 namespace mln 00038 { 00039 00040 namespace make 00041 { 00042 00043 00045 /* 00046 * \param[in] tab C-array of values. 00047 * 00048 * \pre The array dimension N has to be square (N = n * n). 00049 */ 00050 template <typename T, unsigned N> 00051 algebra::h_mat<mlc_sqrt_int(N), T> h_mat(const T (&tab)[N]); 00052 00053 00054 00055 # ifndef MLN_INCLUDE_ONLY 00056 00057 template <typename T, unsigned N> 00058 inline 00059 algebra::h_mat<mlc_sqrt_int(N), T> 00060 h_mat(const T (&tab)[N]) 00061 { 00062 enum { n = mlc_sqrt_int(N) }; 00063 mlc_bool(N == n * n)::check(); 00064 algebra::h_mat<n,T> tmp; 00065 for (unsigned i = 0; i < N; ++i) 00066 tmp(i / n, i % n) = tab[i]; 00067 return tmp; 00068 } 00069 00070 # endif // ! MLN_INCLUDE_ONLY 00071 00072 } // end of namespace mln::make 00073 00074 } // end of namespace mln 00075 00076 00077 #endif // ! MLN_MAKE_H_MAT_HH