Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 2008, 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_W_WINDOW3D_HH 00027 # define MLN_MAKE_W_WINDOW3D_HH 00028 00034 # include <cmath> 00035 00036 # include <mln/core/w_window.hh> 00037 # include <mln/core/alias/dpoint3d.hh> 00038 00039 #include <iostream> 00040 00041 namespace mln 00042 { 00043 00044 namespace make 00045 { 00046 00055 template <typename W, unsigned M> 00056 mln::w_window<mln::dpoint3d, W> w_window3d(W (&weights)[M]); 00057 00058 00059 # ifndef MLN_INCLUDE_ONLY 00060 00061 template <typename W, unsigned M> 00062 inline 00063 mln::w_window<mln::dpoint3d, W> 00064 w_window3d(W (&weights)[M]) 00065 { 00066 int h = unsigned(std::pow(float(M), float(1. / 3))) / 2; 00067 mln_precondition((2 * h + 1) * (2 * h + 1) * (2 * h + 1) == M); 00068 mln::w_window<mln::dpoint3d, W> tmp; 00069 unsigned i = 0; 00070 for (int sli = - h; sli <= h; ++sli) 00071 for (int row = - h; row <= h; ++row) 00072 for (int col = - h; col <= h; ++col) 00073 { 00074 const W& cur = weights[i]; 00075 00076 if (cur != 0) 00077 tmp.insert(cur, dpoint3d(sli, row, col)); 00078 i++; 00079 } 00080 return tmp; 00081 } 00082 00083 # endif // ! MLN_INCLUDE_ONLY 00084 00085 } // end of namespace mln::make 00086 00087 } // end of namespace mln 00088 00089 00090 #endif // ! MLN_MAKE_W_WINDOW3D_HH